Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 233,931 G→A 25.7% S37N (AGT→AAT)  yafT → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273233,9310GA25.7% 435.6 / 173.1 231S37N (AGT→AAT) yafTlipoprotein
Reads supporting (aligned to +/- strand):  ref base G (111/59);  new base A (33/26);  total (144/87)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

Only 100 of 234 total aligned reads displayed.

AATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTC  >  CP009273/233782‑234080
                                                                                                                                                     |                                                                                                                                                     
aaTATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAg                                                                                                                                                       <  2:1815864/150‑1 (MQ=255)
aaTATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAg                                                                                                                                                       <  1:4219399/150‑1 (MQ=255)
   atTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTaa                                                                                                                                                    >  1:2010371/1‑150 (MQ=255)
    ttCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAat                                                                                                                                                   >  1:3039859/1‑150 (MQ=255)
      cTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAatat                                                                                                                                                 >  1:4058290/1‑150 (MQ=255)
       tAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATAtt                                                                                                                                                >  2:4311450/1‑150 (MQ=255)
       tAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATAtt                                                                                                                                                >  1:4222281/1‑150 (MQ=255)
         aTCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATtaa                                                                                                                                              <  1:2170757/150‑1 (MQ=255)
            aTTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTCTTGATGAAAAGAAGAAAAAAGCTCAAGTCACACCAAGTAATATtaataa                                                                                                                                           >  2:4199242/1‑150 (MQ=255)
            aTTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATtaataa                                                                                                                                           >  2:1114158/1‑150 (MQ=255)
               aTGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGAGGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATaaaaa                                                                                                                                        >  2:2520898/1‑150 (MQ=255)
                 gACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAaca                                                                                                                                      >  2:1839402/1‑150 (MQ=255)
                 gACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAaca                                                                                                                                      <  2:4243555/150‑1 (MQ=255)
                  aCAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAacac                                                                                                                                     <  1:786922/150‑1 (MQ=255)
                    aGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTc                                                                                                                                   >  2:2062458/1‑150 (MQ=255)
                    aGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTc                                                                                                                                   >  1:1510744/1‑150 (MQ=255)
                    aGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTc                                                                                                                                   >  2:3420435/1‑150 (MQ=255)
                     ggCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTcc                                                                                                                                  <  1:309294/150‑1 (MQ=255)
                          gggAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGc                                                                                                                             <  2:589105/150‑1 (MQ=255)
                               tgtCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTg                                                                                                                        >  2:3217126/1‑150 (MQ=255)
                               tgtCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTg                                                                                                                        >  2:943576/1‑150 (MQ=255)
                                gtCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTga                                                                                                                       >  1:1772372/1‑150 (MQ=255)
                                    aaTTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTgacaga                                                                                                                   >  1:4241399/1‑150 (MQ=255)
                                     aTTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTgacagac                                                                                                                  <  2:88353/150‑1 (MQ=255)
                                         tGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAg                                                                                                              <  1:1137327/150‑1 (MQ=255)
                                           aaTTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGAt                                                                                                            >  1:617918/1‑150 (MQ=255)
                                           aaTTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGAt                                                                                                            >  1:1720606/1‑150 (MQ=255)
                                           aaTTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGAt                                                                                                            >  1:1064077/1‑150 (MQ=255)
                                             ttCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGAttt                                                                                                          <  1:2894513/150‑1 (MQ=255)
                                             ttCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGAttt                                                                                                          <  1:2218734/150‑1 (MQ=255)
                                             ttCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGAttt                                                                                                          <  1:2894612/150‑1 (MQ=255)
                                               cAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTAt                                                                                                        >  1:2431745/1‑150 (MQ=255)
                                                aaaaaaGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTAtt                                                                                                       >  2:3060820/1‑150 (MQ=255)
                                                       cTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAAt                                                                                                <  2:4339831/150‑1 (MQ=255)
                                                       cTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAAt                                                                                                <  1:957742/150‑1 (MQ=255)
                                                       cTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAAt                                                                                                >  1:851171/1‑150 (MQ=255)
                                                       cTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCACAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAAt                                                                                                <  2:468144/150‑1 (MQ=255)
                                                          tgttgCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGaa                                                                                             >  1:2203103/1‑150 (MQ=255)
                                                           gttgCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTCATGaaa                                                                                            >  2:3843929/1‑150 (MQ=255)
                                                           gttgCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGaaa                                                                                            >  2:2098049/1‑150 (MQ=255)
                                                           gttgCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGaaa                                                                                            >  2:3155928/1‑150 (MQ=255)
                                                                 tatGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCActc                                                                                      >  1:4493346/1‑150 (MQ=255)
                                                                  atGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCActct                                                                                     >  2:4496731/1‑150 (MQ=255)
                                                                   tgtgtgTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTg                                                                                    >  2:685893/1‑150 (MQ=255)
                                                                    gtgtgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTgg                                                                                   >  2:1508844/1‑150 (MQ=255)
                                                                    gtgtgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTgg                                                                                   >  2:1508677/1‑150 (MQ=255)
                                                                      gtgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGcc                                                                                 >  1:3222752/1‑150 (MQ=255)
                                                                      gtgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGcc                                                                                 >  2:1875452/1‑150 (MQ=255)
                                                                       tgtTTTTCTCGCTGCTTGCAGGATGTGCCTCTTAATCTTCTATTGATGAAAAGAATAAAAAAGCTCAAGTCAAAAAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTTTTCGGTAATGAAACCACTCTGGCCg                                                                                <  1:4092730/150‑1 (MQ=255)
                                                                          tATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTAt                                                                             >  1:1120772/1‑150 (MQ=255)
                                                                          tATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTAt                                                                             >  1:668830/1‑150 (MQ=255)
                                                                            ttCTCGCTGCTTGCAGGTTGTGCCTGTGAATCTTCTTTTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAATGATTTATTCGGTAATGAAACCACTCTGGCCGTATcc                                                                           <  2:216969/150‑1 (MQ=255)
                                                                            ttCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATcc                                                                           >  1:2584827/1‑150 (MQ=255)
                                                                                cGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGacg                                                                       >  2:3530845/1‑148 (MQ=255)
                                                                                    gctTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGa                                                                   >  2:2978206/1‑150 (MQ=255)
                                                                                    gctTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGa                                                                   <  2:1543702/150‑1 (MQ=255)
                                                                                    gctTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGa                                                                   >  1:2051429/1‑150 (MQ=255)
                                                                                      ttGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGata                                                                 <  1:4129076/150‑1 (MQ=255)
                                                                                      ttGCAGGATGTGCATCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGata                                                                 <  1:3336750/150‑1 (MQ=255)
                                                                                       tGCAGGATGTGCCTCTGAATCTTCTATTGATTAAAATAATAAAAAAGCTCAAGTCAAACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTTTTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGatat                                                                <  1:4389191/150‑1 (MQ=255)
                                                                                               gtgCCTCTTAATCTTCTATTTATGAAAAGAAGAAAAAAGCTAAAGACACAAAAATTAATATTAATAAAAACCCTCCCCCGCAACTGACAGACAAATATTTATTCGGTAATGACACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTg                                                        <  1:2302622/150‑1 (MQ=255)
                                                                                                  ccTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCAATCTGGCCGTATCCGAAGAGGATATTCAAGCTgcgc                                                     >  2:3577891/1‑150 (MQ=255)
                                                                                                   ctctGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCt                                                    >  1:1693899/1‑150 (MQ=255)
                                                                                                     ctGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTg                                                  >  2:485355/1‑150 (MQ=255)
                                                                                                     ctGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTg                                                  <  2:4241389/150‑1 (MQ=255)
                                                                                                      ttaaTATTCTATTGATGAAAAGAAAAAAAAAGCTCAAGACACACAAAGTAATATTAATAAAAACACACCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGa                                                 <  1:1918885/148‑1 (MQ=255)
                                                                                                      tGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGa                                                 >  2:1797215/1‑150 (MQ=255)
                                                                                                       gAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGAt                                                >  2:554880/1‑150 (MQ=255)
                                                                                                        aaTCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATg                                               >  2:794392/1‑150 (MQ=255)
                                                                                                         aTCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATgg                                              >  2:1402128/1‑150 (MQ=255)
                                                                                                            ttctATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGgaga                                           >  1:3613735/1‑150 (MQ=255)
                                                                                                                aTTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCGTGATGGAGATGAg                                       >  2:3982737/1‑150 (MQ=255)
                                                                                                                 ttGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGt                                      >  2:2287158/1‑150 (MQ=255)
                                                                                                                    atgaAAATAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTcc                                   <  1:2468497/150‑1 (MQ=255)
                                                                                                                    atgaAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGGTATTCAAGCTGCGCTTGATGGAGATGAGTTcc                                   <  2:851171/150‑1 (MQ=255)
                                                                                                                        aaaGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGt                               >  1:3525232/1‑150 (MQ=255)
                                                                                                                        aaaGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGt                               >  1:2383394/1‑150 (MQ=255)
                                                                                                                         aagaagAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGtt                              >  2:2462810/1‑150 (MQ=255)
                                                                                                                         aagaagAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGtt                              >  2:1980689/1‑150 (MQ=255)
                                                                                                                         aagaagAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGtt                              <  1:2541829/150‑1 (MQ=255)
                                                                                                                          agaagaAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTc                             >  2:3079423/1‑150 (MQ=255)
                                                                                                                          agaagaAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTc                             >  2:1102916/1‑150 (MQ=255)
                                                                                                                          agaagaAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTc                             >  2:3388232/1‑150 (MQ=255)
                                                                                                                            aataaaaaaTCAAAATTCACACAAAATAATATTAATAAAAACACTCCCCACCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTccc                           <  2:1273823/147‑1 (MQ=255)
                                                                                                                            aagaaAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGTTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTccc                           <  2:1273658/150‑1 (MQ=255)
                                                                                                                               aaaaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAATATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCTCCtt                        <  1:4065552/150‑1 (MQ=255)
                                                                                                                               aaaaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCtt                        >  1:4535123/1‑150 (MQ=255)
                                                                                                                                aaaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTa                       <  2:3707225/150‑1 (MQ=255)
                                                                                                                                  aaaGCTCAAGTCAAAAAAAATAATATTAATAAAAACACTCCCCAGAAACTGACAGACAAAGATTTATTCGTTAATGACACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAAt                     <  1:3784988/150‑1 (MQ=255)
                                                                                                                                  aaaGCTCAAGACACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACACACAAAGATTTTTTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAAt                     <  1:3783099/150‑1 (MQ=255)
                                                                                                                                    aGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATtc                   >  1:1712328/1‑150 (MQ=255)
                                                                                                                                      cTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTCATGAAACCCCTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCGTTAATtctc                 >  2:2954608/1‑150 (MQ=255)
                                                                                                                                       tCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTcc                <  1:2053974/150‑1 (MQ=255)
                                                                                                                                       tCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTcc                <  1:2054975/150‑1 (MQ=255)
                                                                                                                                        cAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCt               >  2:1134172/1‑150 (MQ=255)
                                                                                                                                         aaGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTg              <  1:1857523/150‑1 (MQ=255)
                                                                                                                                         aaGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTg              <  1:2644937/150‑1 (MQ=255)
                                                                                                                                                  aaaGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTg     >  1:1220177/1‑150 (MQ=255)
                                                                                                                                                   aaGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGt    >  2:1587615/1‑150 (MQ=255)
                                                                                                                                                    aGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGtt   >  2:2814855/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTTC  >  CP009273/233782‑234080

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: