Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,986,289 T→A 83.7% N4Y (AAT→TAT)  pgsA ← phosphatidylglycerophosphate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,986,2890TA83.7% 740.4 / 125.2 258N4Y (AAT→TAT) pgsAphosphatidylglycerophosphate synthetase
Reads supporting (aligned to +/- strand):  ref base T (19/23);  new base A (108/107);  total (128/130)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 259 total aligned reads displayed.

CCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGC  >  CP009273/1986142‑1986435
                                                                                                                                                   |                                                                                                                                                  
ccAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATaaa                                                                                                                                                  >  1:665451/1‑150 (MQ=255)
  aGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAAt                                                                                                                                                >  1:2040059/1‑150 (MQ=255)
    aaaaCCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTg                                                                                                                                              <  1:3217372/150‑1 (MQ=255)
     aaaCCATCGCACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGc                                                                                                                                             <  1:1143827/150‑1 (MQ=255)
     aaaCCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATCAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCCACGGAGGGATCTCAAATTGc                                                                                                                                             >  2:2363985/1‑150 (MQ=255)
       aCCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCAt                                                                                                                                           >  1:4513722/1‑150 (MQ=255)
       aCCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCAt                                                                                                                                           <  1:1370468/150‑1 (MQ=255)
       aCCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCAt                                                                                                                                           <  2:3547093/150‑1 (MQ=255)
       aCCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCAt                                                                                                                                           <  2:1261909/150‑1 (MQ=255)
         cATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATaa                                                                                                                                         >  1:2466779/1‑150 (MQ=255)
          aTCGAACCAGTCAGTCACCGACTCTACGCAGAAAATGAGCGCGGAGTCAAACGTCTACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTATGGATATAAAATTGCATAAt                                                                                                                                        <  2:2766472/150‑1 (MQ=255)
            cGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGa                                                                                                                                      >  1:4539162/1‑150 (MQ=255)
               aCCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACgg                                                                                                                                   <  2:502020/150‑1 (MQ=255)
               aCCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGGAGGGATATAAAATTGCATAATGACgg                                                                                                                                   >  2:4128628/1‑150 (MQ=255)
                ccAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACggg                                                                                                                                  >  1:1089383/1‑150 (MQ=255)
                 cagtcagtCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGt                                                                                                                                 >  2:3266537/1‑150 (MQ=255)
                  agtcagtcACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGAGATAAAATTGCATAATGACGGGTa                                                                                                                                >  2:3097394/1‑150 (MQ=255)
                   gtcagtcaCCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTaa                                                                                                                               <  2:1706984/150‑1 (MQ=255)
                     cagtcaCCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACt                                                                                                                             <  1:2931916/150‑1 (MQ=255)
                       gtcaCCGCCGCGACGCAGAAAAGGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTAt                                                                                                                           <  1:861083/150‑1 (MQ=255)
                             gccgcGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCtgttg                                                                                                                     <  1:402054/150‑1 (MQ=255)
                                   aCGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTa                                                                                                               >  2:3873952/1‑150 (MQ=255)
                                     gCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAg                                                                                                             >  2:896999/1‑150 (MQ=255)
                                      cAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGa                                                                                                            <  1:2020176/150‑1 (MQ=255)
                                       aGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGAt                                                                                                           >  1:2665998/1‑150 (MQ=255)
                                        gAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGAtt                                                                                                          <  2:1419042/150‑1 (MQ=255)
                                          aaaTGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTAc                                                                                                        <  1:4456223/150‑1 (MQ=255)
                                             tGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTAcccc                                                                                                     <  1:539256/150‑1 (MQ=255)
                                             tGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTAcccc                                                                                                     <  1:3639017/150‑1 (MQ=255)
                                              gAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCt                                                                                                    >  1:63102/1‑150 (MQ=255)
                                              gAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCt                                                                                                    <  2:3901586/150‑1 (MQ=255)
                                              gAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCt                                                                                                    >  1:62955/1‑150 (MQ=255)
                                              gAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCt                                                                                                    >  1:4186861/1‑150 (MQ=255)
                                                  gcggcgGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGt                                                                                                <  2:2515945/150‑1 (MQ=255)
                                                  gcggcgGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGt                                                                                                <  2:3790405/150‑1 (MQ=255)
                                                       ggcAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTa                                                                                           >  1:174145/1‑150 (MQ=255)
                                                       ggcAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTa                                                                                           >  1:2052951/1‑150 (MQ=255)
                                                        gcAAACGGTGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTAc                                                                                          <  2:4303918/150‑1 (MQ=255)
                                                          aaaCGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACag                                                                                        >  1:2100415/1‑150 (MQ=255)
                                                          aaaCGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACag                                                                                        <  1:196979/150‑1 (MQ=255)
                                                                  aCCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCa                                                                                >  1:1566045/1‑150 (MQ=255)
                                                                  aCCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCa                                                                                >  1:3559812/1‑150 (MQ=255)
                                                                       gTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGtt                                                                           <  2:2001055/150‑1 (MQ=255)
                                                                        tGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGttt                                                                          >  1:3453613/1‑150 (MQ=255)
                                                                          aCAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCa                                                                        <  2:3696695/150‑1 (MQ=255)
                                                                           cAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCaa                                                                       <  2:1932849/150‑1 (MQ=255)
                                                                            aGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAAc                                                                      >  2:2938300/1‑150 (MQ=255)
                                                                             ggCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACg                                                                     >  1:2208850/1‑150 (MQ=255)
                                                                               cAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGAc                                                                   >  1:2430103/1‑150 (MQ=255)
                                                                                 gATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCa                                                                 <  1:2764678/150‑1 (MQ=255)
                                                                                 gATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCa                                                                 <  1:1005713/150‑1 (MQ=255)
                                                                                   tAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCaga                                                               >  1:4093867/1‑150 (MQ=255)
                                                                                     aaaGACCAATACAAAGAATGGGATAATGATGACAAGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAAATATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTAAACGAACagaag                                                             >  2:4391112/1‑150 (MQ=255)
                                                                                      aaGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCagaaga                                                            >  2:109435/1‑150 (MQ=255)
                                                                                       aGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGAt                                                           <  1:1834101/150‑1 (MQ=255)
                                                                                       aGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGAt                                                           <  1:4111379/150‑1 (MQ=255)
                                                                                        gACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAATATTAACCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATc                                                          >  1:970513/1‑150 (MQ=255)
                                                                                        gACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATc                                                          >  1:4469483/1‑150 (MQ=255)
                                                                                         aCCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAACGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCt                                                         <  1:294055/150‑1 (MQ=255)
                                                                                          ccAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCtt                                                        <  1:4358988/150‑1 (MQ=255)
                                                                                           cAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCttt                                                       >  1:3012752/1‑150 (MQ=255)
                                                                                             aTACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTc                                                     <  2:786496/150‑1 (MQ=255)
                                                                                              tACAAAGAATGTGATAAGGATGACACGGAACAGTTTAAGCACAGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCt                                                    <  2:1815898/150‑1 (MQ=255)
                                                                                                  aaGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCa                                                >  1:3272661/1‑150 (MQ=255)
                                                                                                    gAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGa                                              >  2:3001845/1‑150 (MQ=255)
                                                                                                      aTGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGAcc                                            >  2:431304/1‑150 (MQ=255)
                                                                                                        gggATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCtt                                          <  1:1674477/150‑1 (MQ=255)
                                                                                                        gggATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCtt                                          <  1:896999/150‑1 (MQ=255)
                                                                                                        gggATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCtt                                          <  2:2111265/150‑1 (MQ=255)
                                                                                                         ggATAAGGATGACACGGAACAGTGTAAGCAACGTTGGGATATAAATTTGCATAATGAGGGGTAACTATCTGTTGTCTGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTg                                         <  1:2286179/150‑1 (MQ=255)
                                                                                                         ggATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTg                                         >  1:525720/1‑150 (MQ=255)
                                                                                                           aTAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCg                                       >  1:2721991/1‑150 (MQ=255)
                                                                                                           aTAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCg                                       >  1:457756/1‑150 (MQ=255)
                                                                                                             aaGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGaa                                     >  1:2556968/1‑150 (MQ=255)
                                                                                                             aaGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGaa                                     <  2:547053/150‑1 (MQ=255)
                                                                                                             aaGGATGACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGaa                                     >  2:1140985/1‑150 (MQ=255)
                                                                                                                   gACACGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATAccc                               <  2:1392906/150‑1 (MQ=255)
                                                                                                                     cacGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCgg                             <  2:1718633/150‑1 (MQ=255)
                                                                                                                       cGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCa                           >  2:4065343/1‑150 (MQ=255)
                                                                                                                       cGGAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCa                           >  1:1511221/1‑150 (MQ=255)
                                                                                                                         gAACAGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACt                         <  1:3332493/150‑1 (MQ=255)
                                                                                                                           aCAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACttt                       <  1:2423089/150‑1 (MQ=255)
                                                                                                                             aGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTg                     <  2:3688541/150‑1 (MQ=255)
                                                                                                                             aGTGTAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTg                     <  1:1277057/150‑1 (MQ=255)
                                                                                                                               tgtAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCa                   >  1:3580657/1‑150 (MQ=255)
                                                                                                                                gtAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCaa                  <  2:2238610/150‑1 (MQ=255)
                                                                                                                                 tAAGCAACGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAAt                 <  2:665451/150‑1 (MQ=255)
                                                                                                                                  aaGCAACGTATGGATATAAAATTTCTGAATTACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTAAACGACCAGTAGATCTTTTCTGAGATACCTTGGGAAATACCCTGCACTTGTGCAAtt                <  2:2068298/150‑1 (MQ=255)
                                                                                                                                      aaCGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCt            >  1:2200052/1‑150 (MQ=255)
                                                                                                                                        cGTAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCg          >  2:855351/1‑150 (MQ=255)
                                                                                                                                          tAGGGATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGAc        <  2:4175128/150‑1 (MQ=255)
                                                                                                                                           aGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACg       >  2:3656123/1‑150 (MQ=255)
                                                                                                                                             ggATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCt     <  2:1468436/150‑1 (MQ=255)
                                                                                                                                             ggATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCt     >  2:4409665/1‑150 (MQ=255)
                                                                                                                                             ggATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCt     <  1:3520256/150‑1 (MQ=255)
                                                                                                                                              gATATAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTg    >  1:2262232/1‑150 (MQ=255)
                                                                                                                                               atatAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTgg   >  1:4079843/1‑150 (MQ=255)
                                                                                                                                               atatAAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTgg   >  1:1216349/1‑150 (MQ=255)
                                                                                                                                                tatTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGc  <  1:1611676/150‑1 (MQ=255)
                                                                                                                                                tatTAAATTGCATAATGACGGGTAACTATCTGTGGTCAGTAATATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGc  <  1:4269360/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCAGAAAACCATCGAACCAGTCAGTCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGC  >  CP009273/1986142‑1986435

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