Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 667,017 T→G 18.1% E15D (GAA→GAC holA ← DNA polymerase III, delta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273667,0170TG18.1% 417.2 / 12.5 189E15D (GAA→GACholADNA polymerase III, delta subunit
Reads supporting (aligned to +/- strand):  ref base T (70/75);  new base G (36/2);  total (112/77)
Fisher's exact test for biased strand distribution p-value = 3.05e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 192 total aligned reads displayed.

AGATCGCATTCCAGTCAGTGTTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATGCTTGG  >  CP009273/666869‑667165
                                                                                                                                                    |                                                                                                                                                    
aGATCGCATTCCAGTCAGTGTTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCt                                                                                                                                                      <  1:2633444/149‑1 (MQ=255)
  aTCGCATTCCAGTCAGTGTTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTc                                                                                                                                                    <  1:3238494/149‑1 (MQ=255)
  aTCGCATTCCAGTCAGTGTTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTc                                                                                                                                                    >  1:4362956/1‑149 (MQ=255)
    cGCCTTCCAGTCAGTGTTGGGATCAATGGAAAAAGTGTGGTGGTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCAt                                                                                                                                                  <  1:3306737/149‑1 (MQ=255)
       aTTCCAGTCAGTGTTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGa                                                                                                                                               >  1:3552063/1‑149 (MQ=255)
       aTTCCAGTCAGTGTTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCGTCAGTGa                                                                                                                                               >  1:4137433/1‑149 (MQ=255)
        ttCCAGTCAGTGTTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCGTCATTGAg                                                                                                                                              >  2:1994261/1‑149 (MQ=255)
        ttCCAGTCAGTGTTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCATCATTGAg                                                                                                                                              >  1:2101267/1‑149 (MQ=255)
          ccAGTCAGTGTTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCt                                                                                                                                            <  1:1492705/149‑1 (MQ=255)
              tcagtGTTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTGCTTGTGCTgcgc                                                                                                                                        >  1:1191472/1‑149 (MQ=255)
                 gtgtTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGcgcgcg                                                                                                                                     >  1:2808098/1‑149 (MQ=255)
                  tgtTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCgg                                                                                                                                    >  2:1717135/1‑149 (MQ=255)
                  tgtTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCgg                                                                                                                                    >  2:1926272/1‑149 (MQ=255)
                  tgtTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCATCATTGAGCTGCGCGCgg                                                                                                                                    >  1:676924/1‑149 (MQ=255)
                        gATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAgttgt                                                                                                                              >  1:1325471/1‑149 (MQ=255)
                             aTGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCgg                                                                                                                         >  2:1870930/1‑149 (MQ=255)
                             aTGGAAAAAGTGTGGTGTGCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCGGCAGTGAGCTGCGCGCGGAGTTGTTCCgg                                                                                                                         >  2:2353006/1‑149 (MQ=255)
                              tGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCggg                                                                                                                        <  1:1882085/149‑1 (MQ=255)
                               ggAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGt                                                                                                                       <  1:167709/149‑1 (MQ=255)
                                   aaaGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACaa                                                                                                                   >  2:3052037/1‑149 (MQ=255)
                                   aaaGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACaa                                                                                                                   <  1:2362198/149‑1 (MQ=255)
                                    aaGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAAc                                                                                                                  <  2:2260547/149‑1 (MQ=255)
                                     aGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCGGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAAcc                                                                                                                 <  1:3022654/149‑1 (MQ=255)
                                         tggtgTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCGTCATTGAGCTGCGCGCGGAGTTGGTCCGGGTACAACCCaaa                                                                                                             >  1:538281/1‑148 (MQ=255)
                                            tgtTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCGTCCTTGAGCGGCGCGCGGAGTTGGTCCGGGTACAACCGAatcat                                                                                                          >  1:878102/1‑149 (MQ=255)
                                             gtTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAatcatc                                                                                                         <  1:209727/149‑1 (MQ=255)
                                             gtTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCAGTGAGCTGCGCGCGGAGTTTTTCCGGGTGCAACCCAatcatc                                                                                                         >  1:2991388/1‑149 (MQ=255)
                                              ttcttcGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCATCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAAtcatca                                                                                                        >  1:1091641/1‑149 (MQ=255)
                                              gtcttcGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAAtcatca                                                                                                        <  2:2753302/148‑1 (MQ=255)
                                               tcttcGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAg                                                                                                       <  2:86799/149‑1 (MQ=255)
                                               tcttcGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAg                                                                                                       <  2:4111509/149‑1 (MQ=255)
                                                cttcGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGt                                                                                                      <  2:1981188/149‑1 (MQ=255)
                                                cttcGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGt                                                                                                      <  2:4202396/149‑1 (MQ=255)
                                                cttcGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCGGCCGTTAGCGGCGCGCGGAGTTGTTCCGGGTACAACCGAATAATCAgg                                                                                                      >  1:532671/1‑148 (MQ=255)
                                                 ttcGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGGGCCGCCCGTCATTGATCTGCGCGCGGAGTTGTTACGGGTACAAACCAATCAAAAGtt                                                                                                     >  1:1460066/1‑149 (MQ=255)
                                                     aaTCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTTAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTAccc                                                                                                 >  1:258487/1‑149 (MQ=255)
                                                     aaTCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTAccc                                                                                                 <  1:8273/149‑1 (MQ=255)
                                                       tCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGTGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAg                                                                                               >  2:4278418/1‑149 (MQ=255)
                                                       tCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAg                                                                                               >  2:2481077/1‑149 (MQ=255)
                                                                  cAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCCAGTAAAAGATACGCCGCGCGCAGCCCGTCAGTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGgaga                                                                                    <  2:3713529/149‑1 (MQ=255)
                                                                  cAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGgaga                                                                                    <  2:180677/149‑1 (MQ=255)
                                                                       aCCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGTTACGCCGCGCGCTGCCCTTCATTGAGCTGGGCGCGGAGGTGTTCCGGGTACAACCGTTTCATCAGTTTCCCAGCGTGGTGGAGACGCGt                                                                               >  1:344413/1‑149 (MQ=255)
                                                                       aCCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCGTCATTTAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGt                                                                               >  1:561903/1‑149 (MQ=255)
                                                                       aCCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCGGCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGGGGTGGAGACGCGt                                                                               >  2:3371577/1‑149 (MQ=255)
                                                                          tGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGTGCTGCGCGCGGAGTTGTTACGGGTACAACCGAATCATCAGTTACCCCGCGGGGGGGAGACGCGTGCa                                                                            >  2:1542082/1‑149 (MQ=255)
                                                                           gACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCGTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAg                                                                           >  2:1205414/1‑149 (MQ=255)
                                                                               aaCAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCACTTAACCAGCGTGGTGGAGACGCGTGCAGGCGt                                                                       >  1:247242/1‑149 (MQ=255)
                                                                                 cAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTg                                                                     >  2:3910755/1‑149 (MQ=255)
                                                                                   gCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGCATCATCAGTTACCCAGCGTGGGGGAGACGCGTGCAGGCGTTGcc                                                                   >  2:3663358/1‑149 (MQ=255)
                                                                                   gCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGcc                                                                   >  2:3663450/1‑149 (MQ=255)
                                                                                   gCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGcc                                                                   >  2:1438941/1‑149 (MQ=255)
                                                                                      tCCTGGCTGTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGa                                                                <  2:1469151/149‑1 (MQ=255)
                                                                                       ccgggCTTTCCTGCACTCACAGAGGATCGTTACCAAGTAACAGATACGCCACGCGCAGCCCGTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAg                                                               <  1:3878466/146‑1 (MQ=255)
                                                                                       ccTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAg                                                               <  2:2879901/149‑1 (MQ=255)
                                                                                         tGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTa                                                             <  1:4012793/149‑1 (MQ=255)
                                                                                          ggCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCGTCATTTAGCTGCGCGCGGGGTTGTTCCGGGTACAACCGCATCATAAGTTAACCAGCCGGGTGGAGAAGCGTGCAGGCGTTGCCGGCGTAt                                                            >  1:873886/1‑149 (MQ=255)
                                                                                           gCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTGATTGTGGTGCGCGCGGTGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGTCGCGTGCCGGCGTTGCCGGAGTATc                                                           >  2:2675480/1‑149 (MQ=255)
                                                                                                 cTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTtgt                                                     <  1:4113447/149‑1 (MQ=255)
                                                                                                    cAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTCGTGGt                                                  <  1:3118168/149‑1 (MQ=255)
                                                                                                     aaTAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTc                                                 <  1:3392399/149‑1 (MQ=255)
                                                                                                        aaCAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGAc                                              <  2:1673626/149‑1 (MQ=255)
                                                                                                        aaCAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGAc                                              >  2:3062126/1‑149 (MQ=255)
                                                                                                         aCAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCGTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACg                                             >  1:2164618/1‑149 (MQ=255)
                                                                                                           agagGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTc                                           <  1:3579828/149‑1 (MQ=255)
                                                                                                            gagGATCGTTACCAAGTAAAAGATACGCCGCGCGCACCCCTTCCATGAGCTGCGCGCGGAGTTGGTCCGGGTACAACCCAATCATCAGTTACCCAGCGTGGGGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCt                                          >  1:2911278/1‑149 (MQ=255)
                                                                                                              ggATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGt                                        <  2:2146874/149‑1 (MQ=255)
                                                                                                              ggATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCGGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTgg                                        <  1:3165688/149‑2 (MQ=255)
                                                                                                                aTCGTTACCAAGTAAAAGATACGCCGCGCGCCGCCCTTCATTGAGCGGCGCGCGGAGGTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGGGGTGGAGACGCGTGCAGGCGGTGGCGGAGTATCCGTTGTGGTCGACGTCTGttc                                      >  2:2588123/1‑149 (MQ=255)
                                                                                                                aTCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGttc                                      >  2:3245117/1‑149 (MQ=255)
                                                                                                                 tCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGttct                                     <  1:283383/149‑1 (MQ=255)
                                                                                                                 tCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTACGGGTACACCCGAATCATCAGGTACCCAGCGTGGTGGAGACGCGGGCAGGCGTTGGCGGAGTATCCGTTGTGGTCGACGTCTGttct                                     >  2:250122/1‑149 (MQ=255)
                                                                                                                   gTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCGGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGttcttc                                   >  2:110921/1‑149 (MQ=255)
                                                                                                                    ttACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTtcg                                  <  1:364518/149‑1 (MQ=255)
                                                                                                                     tACCAAGTAAAAGATACGCCGCGCGCAGCCCGTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTtcgt                                 >  2:1831637/1‑149 (MQ=255)
                                                                                                                      aCCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTtcgtc                                <  2:2086819/149‑1 (MQ=255)
                                                                                                                      aCCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTtcgtc                                <  1:3740975/149‑1 (MQ=255)
                                                                                                                      aCCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTtcgtc                                <  1:4228094/149‑1 (MQ=255)
                                                                                                                       ccAAGTAAAAGATACGCCGCGCGCAGCCCGTCATTGTGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCACCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGGCGACGTCTGTTCTtcgtcg                               >  1:3983313/1‑149 (MQ=255)
                                                                                                                        cAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCgg                              <  1:3944886/149‑1 (MQ=255)
                                                                                                                        cAAGTAAAAGATACGCCGCGCGCAGCCCGTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGCATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCgg                              >  2:2660319/1‑149 (MQ=255)
                                                                                                                        cAAGTAAAAGATACGCCGCGCGCAGCCCGTCATTGAGCGGCGCGCGGAGTTGTTCCGGGGACAACCGAATCATCACTTACCCAGCGGGGTGGCGACGCGTGCAGGCGGTGCCGGCGTATCCGGTGTGGGCGACGGCTGTTCTTCGTCgg                              >  2:458787/1‑149 (MQ=255)
                                                                                                                          aGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGaa                            <  1:3052040/149‑1 (MQ=255)
                                                                                                                               aaGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAat                       <  1:1101160/149‑1 (MQ=255)
                                                                                                                                aGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAata                      <  1:293918/149‑1 (MQ=255)
                                                                                                                                 gATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAatat                     <  1:1480272/149‑1 (MQ=255)
                                                                                                                                 gATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAAACATAAGTTAACCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTATTCGTCGGAACGAatat                     >  2:3527326/1‑149 (MQ=255)
                                                                                                                                  aTACGCCGCGCGGAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTAAAAACGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTATTCGTCGGAACGAatata                    >  2:2741986/1‑148 (MQ=255)
                                                                                                                                   tACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATcc                   >  2:4021746/1‑149 (MQ=255)
                                                                                                                                    aCGCCGCGCGCAGCCCTTCATTTAGCTGCGCGCGGAGTTGTTCCGGGTACAAACGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTATTAGTAGGAACGAATATACg                  >  1:4251545/1‑149 (MQ=255)
                                                                                                                                         gcgcgcAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCgcagca             <  2:676924/149‑1 (MQ=255)
                                                                                                                                         gcgcgcAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCgcagca             <  1:2128671/149‑1 (MQ=255)
                                                                                                                                           gcgcAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACg           >  2:590960/1‑149 (MQ=255)
                                                                                                                                             gcAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGa         >  2:3107723/1‑149 (MQ=255)
                                                                                                                                             gcAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGa         >  1:1331107/1‑149 (MQ=255)
                                                                                                                                               aGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATg       <  2:3187444/149‑1 (MQ=255)
                                                                                                                                               aGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATg       <  2:3423397/149‑1 (MQ=255)
                                                                                                                                                 cccTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATGCt     <  1:1419629/149‑1 (MQ=255)
                                                                                                                                                  ccTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATGCtt    <  2:2944882/149‑1 (MQ=255)
                                                                                                                                                  ccTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATGCtt    >  2:3833096/1‑149 (MQ=255)
                                                                                                                                                  ccTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATGCtt    >  2:1612861/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGATCGCATTCCAGTCAGTGTTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACAGCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATGCTTGG  >  CP009273/666869‑667165

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: