Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,432,682 G→T 9.0% E162* (GAA→TAA)  tamA → translocation and assembly module for autotransporter export, outer membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,432,6820GT 9.0% 907.1 / 48.5 245E162* (GAA→TAA) tamAtranslocation and assembly module for autotransporter export, outer membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (86/137);  new base T (9/13);  total (95/150)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 246 total aligned reads displayed.

TATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCC  >  CP009273/4432533‑4432830
                                                                                                                                                     |                                                                                                                                                    
 aTTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGa                                                                                                                                                     <  1:3742727/150‑1 (MQ=255)
     cgcgGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAAttt                                                                                                                                                 <  1:419266/150‑1 (MQ=255)
       cggcggcgCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTAc                                                                                                                                               <  1:379626/150‑1 (MQ=255)
         gcggcgCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCa                                                                                                                                             <  2:665272/150‑1 (MQ=255)
              gcgcgGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAgcg                                                                                                                                        <  1:4716529/150‑1 (MQ=255)
                gcgGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCa                                                                                                                                      <  1:33342/150‑1 (MQ=255)
                   gACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAATGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCt                                                                                                                                   >  2:303674/1‑150 (MQ=255)
                     ccGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTgg                                                                                                                                 >  2:2217422/1‑150 (MQ=255)
                      cGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTggg                                                                                                                                >  1:1636969/1‑150 (MQ=255)
                      cGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTggg                                                                                                                                >  1:989645/1‑150 (MQ=255)
                       gATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGc                                                                                                                               <  1:3104326/150‑1 (MQ=255)
                         tAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCAt                                                                                                                             >  2:504483/1‑150 (MQ=255)
                          aaaGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCAtt                                                                                                                            >  1:865434/1‑150 (MQ=255)
                             gACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTgcg                                                                                                                         <  2:4562949/150‑1 (MQ=255)
                               cTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCt                                                                                                                       >  2:308839/1‑150 (MQ=255)
                                tATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTc                                                                                                                      <  1:570031/150‑1 (MQ=255)
                                tATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTc                                                                                                                      <  1:2662049/150‑1 (MQ=255)
                                tATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTc                                                                                                                      <  1:289918/150‑1 (MQ=255)
                                tATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTc                                                                                                                      <  2:4302859/150‑1 (MQ=255)
                                     gAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCt                                                                                                                 <  1:3844035/150‑1 (MQ=255)
                                     gAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCt                                                                                                                 <  1:3418209/150‑1 (MQ=255)
                                     gAAAGTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCt                                                                                                                 <  2:1096893/150‑1 (MQ=255)
                                      aacttGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTg                                                                                                                <  2:3718729/147‑1 (MQ=255)
                                       aaTTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGc                                                                                                               <  1:39741/150‑1 (MQ=255)
                                          tGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATa                                                                                                            <  1:3683259/150‑1 (MQ=255)
                                          tGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATa                                                                                                            <  1:4295691/150‑1 (MQ=255)
                                           gCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATaa                                                                                                           <  2:2252975/150‑1 (MQ=255)
                                            cTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCTAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATaaa                                                                                                          >  2:952934/1‑150 (MQ=255)
                                             tCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAg                                                                                                         <  1:2283396/150‑1 (MQ=255)
                                             tCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAg                                                                                                         <  2:2857921/150‑1 (MQ=255)
                                               gATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGcc                                                                                                       >  2:1718983/1‑150 (MQ=255)
                                                aTACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCTAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCt                                                                                                      <  1:941921/150‑1 (MQ=255)
                                                 tACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCtt                                                                                                     >  1:1276411/1‑150 (MQ=255)
                                                  aCTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTc                                                                                                    <  1:176615/150‑1 (MQ=255)
                                                    tCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTg                                                                                                  <  1:1108606/150‑1 (MQ=255)
                                                      gCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTggg                                                                                                <  2:205082/150‑1 (MQ=255)
                                                         cGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGata                                                                                             <  1:2020198/150‑1 (MQ=255)
                                                          ggCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCTAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGatat                                                                                            <  2:1807354/150‑1 (MQ=255)
                                                          ggCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGatat                                                                                            <  1:4752208/150‑1 (MQ=255)
                                                            cTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCTAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTg                                                                                          >  2:1346714/1‑150 (MQ=255)
                                                             tATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGa                                                                                         >  2:46735/1‑150 (MQ=255)
                                                              aTTGGCACGGTACTGAACCAGGGCGATTATGAAACTTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGAt                                                                                        <  2:3689457/150‑1 (MQ=255)
                                                               ttGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGAtt                                                                                       <  1:131951/150‑1 (MQ=255)
                                                                tGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATta                                                                                      >  1:2765086/1‑150 (MQ=255)
                                                                 ggCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATtat                                                                                     <  1:4465730/150‑1 (MQ=255)
                                                                 ggCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATtat                                                                                     <  2:656801/150‑1 (MQ=255)
                                                                    aCGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAAc                                                                                  >  2:3131055/1‑150 (MQ=255)
                                                                    aCGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAAc                                                                                  >  1:4170208/1‑150 (MQ=255)
                                                                    aCGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAAc                                                                                  >  2:761919/1‑150 (MQ=255)
                                                                    aCGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAAc                                                                                  >  1:3392417/1‑150 (MQ=255)
                                                                       gTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCTAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGt                                                                               <  1:4566647/150‑1 (MQ=255)
                                                                        tACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTg                                                                              <  2:4597192/150‑1 (MQ=255)
                                                                        tACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTg                                                                              <  2:2892058/150‑1 (MQ=255)
                                                                        tACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTg                                                                              <  2:2468273/150‑1 (MQ=255)
                                                                                cAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGt                                                                      >  1:1936134/1‑150 (MQ=255)
                                                                                 aGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGtt                                                                     <  2:1449263/150‑1 (MQ=255)
                                                                                    gCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTAcc                                                                  <  1:3531345/150‑1 (MQ=255)
                                                                                         tATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCttt                                                             >  2:2346982/1‑150 (MQ=255)
                                                                                         tATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATAGCGCTCGGCATGCATAAAGACTTCTGGGATATTGATTATAACAGTGGAGAACGTAAACGAttt                                                             >  2:4281274/1‑150 (MQ=255)
                                                                                          aTGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTACTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTg                                                            >  2:1972346/1‑150 (MQ=255)
                                                                                            gAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTggg                                                          >  2:1690822/1‑150 (MQ=255)
                                                                                            gAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTggg                                                          >  2:4076207/1‑150 (MQ=255)
                                                                                               aaTTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCAt                                                       <  1:2724219/150‑1 (MQ=255)
                                                                                               aaTTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCAt                                                       <  1:3582622/150‑1 (MQ=255)
                                                                                                aTTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCAtg                                                      >  1:2805116/1‑150 (MQ=255)
                                                                                                 tttCAAAAAGTCCTTAACCAGCAGTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCAtgt                                                     <  2:4436948/150‑1 (MQ=255)
                                                                                                     aaaaaGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGAcc                                                 >  2:3892838/1‑150 (MQ=255)
                                                                                                       aaaGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCtt                                               >  2:3263008/1‑150 (MQ=255)
                                                                                                       aaaGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCtt                                               >  1:924414/1‑150 (MQ=255)
                                                                                                        aaGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCttt                                              <  2:1842917/150‑1 (MQ=255)
                                                                                                           tCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGaa                                           <  1:2500267/150‑1 (MQ=255)
                                                                                                             cTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAgg                                         <  1:3926151/150‑1 (MQ=255)
                                                                                                               tAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGAt                                       <  1:1078937/150‑1 (MQ=255)
                                                                                                                aaCCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATc                                      <  1:2217422/150‑1 (MQ=255)
                                                                                                                 aCCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATca                                     <  2:3410798/150‑1 (MQ=255)
                                                                                                                  ccAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATcac                                    >  1:1665880/1‑150 (MQ=255)
                                                                                                                   cAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATcaca                                   <  1:3196681/150‑1 (MQ=255)
                                                                                                                     gCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACaaa                                 <  1:1617857/150‑1 (MQ=255)
                                                                                                                      taTTGCGTTGCGTAAAGGTTATTTCGATAGCTAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAAt                                >  2:2640461/2‑150 (MQ=255)
                                                                                                                           cgttgcgtAAAGGTTATTTCGATAGCTAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCcgcg                           >  1:731657/1‑150 (MQ=255)
                                                                                                                           cgttgcgtAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCcgcg                           <  1:2345651/150‑1 (MQ=255)
                                                                                                                            gttgcgtAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGa                          <  1:3179093/150‑1 (MQ=255)
                                                                                                                               gcgtAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGa                       <  2:3862686/150‑1 (MQ=255)
                                                                                                                               gcgtAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGa                       >  1:797461/1‑150 (MQ=255)
                                                                                                                                cgtAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGaa                      >  2:1917685/1‑150 (MQ=255)
                                                                                                                                  tAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATa                    >  2:2220375/1‑150 (MQ=255)
                                                                                                                                  tAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATa                    >  1:1385273/1‑150 (MQ=255)
                                                                                                                                    aaGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATAcc                  <  2:871743/150‑1 (MQ=255)
                                                                                                                                       gTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGc               >  2:1062590/1‑150 (MQ=255)
                                                                                                                                        ttATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAACGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCa              <  2:3967244/150‑1 (MQ=255)
                                                                                                                                            ttCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAAt          <  2:747287/150‑1 (MQ=255)
                                                                                                                                              cGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCt        <  1:1490128/150‑1 (MQ=255)
                                                                                                                                              cGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCt        >  1:2589829/1‑150 (MQ=255)
                                                                                                                                               gATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTg       >  1:1136528/1‑150 (MQ=255)
                                                                                                                                               gATAGCGAATTTACCAAAGCGCAGCTGGGCACTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTg       >  2:4128593/1‑150 (MQ=255)
                                                                                                                                                aTAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTgg      >  2:2438832/1‑150 (MQ=255)
                                                                                                                                                  aGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTg    <  2:1534445/150‑1 (MQ=255)
                                                                                                                                                   gCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGc   >  2:4355715/1‑150 (MQ=255)
                                                                                                                                                   gCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGc   <  1:572828/150‑1 (MQ=255)
                                                                                                                                                    cGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGcc  >  2:2884402/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCC  >  CP009273/4432533‑4432830

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