Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,165,868 +C 99.0% intergenic (‑173/‑211) RCAP_rcc02007 ← / → tig phage integrase/trigger factor

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,165,8681.→C99.0% 2.9 0.0 ‑0.8383intergenic (‑173/‑211)RCAP_rcc02007/tigphage integrase/trigger factor
Reads supporting (aligned to +/- strand):  new base (214/164):  ref base (2/2):  total (217/166)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.76e-01

Only 100 of 421 total aligned reads displayed.

TCTTTGCGGCGTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTT.CTCTG.CCTG.C.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGAAAACGGCCCGGAAAA  >  NC_014034/2165770‑2165955
                                                                                                    |                                                                                         
tCTTTGCGGCGTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACGCGCACCATCTGACCTTTCCTctg.c                                                                                          <  2:6096161/101‑1 (MQ=255)
tCTTTGCGGCGTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTctg.c                                                                                          <  1:7745028/101‑1 (MQ=255)
 cTTTGCGGCGTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTctg.cc                                                                                         <  2:3180534/101‑1 (MQ=255)
 cTTTGCGGCGTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTctg.cc                                                                                         >  1:2454832/1‑101 (MQ=255)
 cTTTGCGGCGTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTctg.cc                                                                                         >  2:6880830/1‑101 (MQ=255)
 cTTTGCGGCGTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTctg.cc                                                                                         >  1:8499755/1‑101 (MQ=255)
 cTTTGCGGCGTCGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTctg.cc                                                                                         >  1:2087277/1‑101 (MQ=255)
  tttGCGGCGTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTctg.ccc                                                                                        <  1:9165603/101‑2 (MQ=255)
    tGCGGCGTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTctg.ccctg                                                                                      >  2:3039059/1‑98 (MQ=255)
     gcggcgTAGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTctg.ccctgc                                                                                     >  1:8681572/1‑97 (MQ=255)
      cggcgTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTctg.ccctgcc                                                                                    <  2:8045636/101‑6 (MQ=255)
        gcgTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTG.c.ccc                                                                                 <  1:5045124/101‑2 (MQ=255)
        gcgTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTG.c.ccc                                                                                 <  1:4849527/101‑2 (MQ=255)
        gcgTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTG.c.ccc                                                                                 <  1:3242908/101‑2 (MQ=255)
        gcgTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTG.c.ccc                                                                                 >  2:1353732/1‑100 (MQ=255)
        gcgTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTG.c.ccc                                                                                 <  1:1250922/101‑2 (MQ=255)
         cgTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTG.c.cccg                                                                                >  2:3605308/1‑99 (MQ=255)
          gTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTG.c.cccgc                                                                               >  2:4992819/1‑98 (MQ=255)
          gTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTG.c.cccgc                                                                               >  1:8769204/1‑98 (MQ=255)
           ttGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTG.c.cccgcg                                                                              <  2:1628655/101‑5 (MQ=255)
             gagaTTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCt                                                                             >  1:5897448/1‑101 (MQ=255)
             gagaTTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCt                                                                             >  2:3121862/1‑101 (MQ=255)
             gagaTTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCt                                                                             <  1:2608299/101‑1 (MQ=255)
             gagaTTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCt                                                                             >  2:7494227/1‑101 (MQ=255)
               gaTTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGc                                                                           >  1:7074813/1‑101 (MQ=255)
                 tttCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCg                                                                         >  2:6499927/1‑101 (MQ=255)
                 tttCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCg                                                                         >  1:4356711/1‑101 (MQ=255)
                 tttCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCg                                                                         >  2:2376015/1‑101 (MQ=255)
                     ccGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGtctc                                                                     <  2:5065328/101‑1 (MQ=255)
                      cGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCg                                                                    <  2:2740820/101‑1 (MQ=255)
                       ggggCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCACCTGACCTTTCCCCTGCCCTG.CTCCGCGCTCCCGtccggg                                                                   >  2:4276711/1‑97 (MQ=255)
                         ggCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCgggg                                                                 <  1:1353764/101‑1 (MQ=255)
                         ggCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCgggg                                                                 <  1:2874530/101‑1 (MQ=255)
                         ggCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCgggg                                                                 >  2:447999/1‑101 (MQ=255)
                         ggCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCgggg                                                                 >  2:3735991/1‑101 (MQ=255)
                             cAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCgtc                                                             >  2:6038983/1‑101 (MQ=255)
                             cAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCgtc                                                             >  1:1568420/1‑101 (MQ=255)
                              agagTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCgtct                                                            >  1:3118479/1‑101 (MQ=255)
                               gagTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCgtctg                                                           <  1:6674036/101‑1 (MQ=255)
                               gagTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCgtctg                                                           >  1:5305773/1‑101 (MQ=255)
                               gagTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCgtctg                                                           >  1:1878386/1‑101 (MQ=255)
                                agTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCgtctgt                                                          >  1:71361/1‑101 (MQ=255)
                                agTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCgtctgt                                                          >  2:7116582/1‑101 (MQ=255)
                                     aCCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTggg                                                     >  1:5826822/1‑101 (MQ=255)
                                      ccTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGc                                                    <  2:4356585/101‑1 (MQ=255)
                                       cTGAGACGGCATGCAACACGGCACTGGTGCAGCCCCCCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGCCTCGGGGGGTCTGTCCGGCCg                                                   >  1:3918998/1‑101 (MQ=255)
                                        tGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGa                                                  >  1:8880287/1‑101 (MQ=255)
                                         gagaCGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGaa                                                 >  2:838021/1‑101 (MQ=255)
                                         gagaCGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGaa                                                 >  1:170361/1‑101 (MQ=255)
                                          agaCGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAAt                                                <  1:4507180/101‑1 (MQ=255)
                                          agaCGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAAt                                                <  2:4112025/101‑1 (MQ=255)
                                           gaCGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATc                                               <  2:7939942/101‑1 (MQ=255)
                                           gaCGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATc                                               <  2:7815199/101‑1 (MQ=255)
                                            aCGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCt                                              <  1:1047373/101‑1 (MQ=255)
                                              ggCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCttt                                            <  1:838034/101‑1 (MQ=255)
                                                cATGCAACACGGCACTGGTGCAGACCACCCGCGCCATCTGGCCTTTCCTCTGCCCTGCC.CCGCGCTGCCGGCTGGGGGCGTCTGGCTGGGCGAATCTTTTc                                          <  2:3918883/101‑1 (MQ=255)
                                                cATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTc                                          <  2:9096230/101‑1 (MQ=255)
                                                 aTGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCg                                         <  1:3441945/101‑1 (MQ=255)
                                                   gCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGAt                                       >  1:7498758/1‑101 (MQ=255)
                                                    cAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATc                                      <  1:4599204/101‑1 (MQ=255)
                                                     aaCACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCt                                     <  1:7646427/101‑1 (MQ=255)
                                                      acacGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCtt                                    <  2:5537651/101‑1 (MQ=255)
                                                       cacGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCttt                                   >  1:4171033/1‑101 (MQ=255)
                                                       cacGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCttt                                   >  2:8235395/1‑101 (MQ=255)
                                                        acGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTc                                  <  2:568951/101‑1 (MQ=255)
                                                         cGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCt                                 >  2:620585/1‑101 (MQ=255)
                                                         cGGCACTGGTGCAGACCACCCGCACCATCTGACCTTTACTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCt                                 <  1:6845175/101‑1 (MQ=255)
                                                            cACTGGTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGcc                              <  1:6375580/101‑1 (MQ=255)
                                                                 gTGCATCCCACCCGCGCCATCTGCCCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGcc                         <  1:551155/101‑1 (MQ=255)
                                                                 gTGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGcc                         >  1:6026773/1‑101 (MQ=255)
                                                                  tGCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGGTCTTTCTGCCgcgctg                        >  2:7339039/1‑99 (MQ=255)
                                                                   ggaGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCgg                       <  2:3308251/99‑1 (MQ=255)
                                                                   gCAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCgg                       >  2:2510597/1‑101 (MQ=255)
                                                                    cAGACCACCCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCggg                      >  2:8601671/1‑101 (MQ=255)
                                                                       accaccCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGAcc                   <  1:5414869/101‑1 (MQ=255)
                                                                       accaccCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGAcc                   <  1:371372/101‑1 (MQ=255)
                                                                        ccaccCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCg                  <  2:1111855/101‑1 (MQ=255)
                                                                        ccaccCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCg                  >  1:6063449/1‑101 (MQ=255)
                                                                        ccaccCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCg                  >  1:8265681/1‑101 (MQ=255)
                                                                         caccCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGa                 >  2:5471108/1‑101 (MQ=255)
                                                                         caccCGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGa                 <  2:7211301/101‑1 (MQ=255)
                                                                           cccGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGaaa               >  2:1956420/1‑101 (MQ=255)
                                                                            ccGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGaaaa              >  2:1765408/1‑101 (MQ=255)
                                                                            ccGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGaaaa              >  2:5243509/1‑101 (MQ=255)
                                                                             cGCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGAAAAc             <  1:4287655/101‑1 (MQ=255)
                                                                              gCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGACAACg            >  1:7300872/1‑101 (MQ=255)
                                                                              gCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGAAAACg            >  2:8476928/1‑101 (MQ=255)
                                                                              gCACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGAAAACg            >  1:6924896/1‑101 (MQ=255)
                                                                               cACCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGAAAACgg           >  1:4628264/1‑101 (MQ=255)
                                                                               cACCATCTGACCGTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGAAAACgg           >  1:8072222/1‑101 (MQ=255)
                                                                                aCCATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGAAAACGGc          <  2:6340040/101‑1 (MQ=255)
                                                                                 ccATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGAAAACGGcc         >  1:7974755/1‑101 (MQ=255)
                                                                                 ccATCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGAAAACGGcc         <  2:1873439/101‑1 (MQ=255)
                                                                                   aTCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGCCCGAAAACGGCCCg       >  2:5825118/1‑101 (MQ=255)
                                                                                   aTCTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGAAAACGGCCCg       <  2:4679018/101‑1 (MQ=255)
                                                                                     cTGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGAAAACGGCCCGGa     <  1:6542296/101‑1 (MQ=255)
                                                                                      tGACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGAAAACGGCCCGGaa    >  2:2209229/1‑101 (MQ=255)
                                                                                       gACCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGAAAACGGCCCGGaaa   <  2:6977155/101‑1 (MQ=255)
                                                                                        aCCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGAAAACGGCCCGGaaaa  >  1:6291585/1‑101 (MQ=255)
                                                                                        aCCTTTCCTCTGCCCTGCC.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGAAAACGGCCCGGaaaa  <  1:2361612/101‑1 (MQ=255)
                                                                                                    |                                                                                         
TCTTTGCGGCGTTGAGATTTCCCGGGGCGCAGAGTCAACCTGAGACGGCATGCAACACGGCACTGGTGCAGACCACCCGCACCATCTGACCTTT.CTCTG.CCTG.C.CCGCGCTGCCGTCTCGGGGCGTCTGTCTGGGCGAATCTTTTCGATCTTTCTGCCGCGCCGGGACCGAAAACGGCCCGGAAAA  >  NC_014034/2165770‑2165955

Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 31 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG