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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273233,9310GA35.0% 169.9 / 159.0 152S37N (AGT→AAT) yafTlipoprotein
Reads supporting (aligned to +/- strand):  ref base G (68/30);  new base A (30/23);  total (98/54)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

Only 100 of 157 total aligned reads displayed.

AATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTT  >  CP009273/233782‑234079
                                                                                                                                                     |                                                                                                                                                    
aaTATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAg                                                                                                                                                      >  1:2877724/1‑150 (MQ=255)
   atTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTaa                                                                                                                                                   <  1:153352/150‑1 (MQ=255)
     tCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAataata                                                                                                                                                 <  1:930381/150‑1 (MQ=255)
       tAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATAtt                                                                                                                                               <  1:4631616/150‑1 (MQ=255)
       tAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATAtt                                                                                                                                               <  1:3927360/150‑1 (MQ=255)
        aaTCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATta                                                                                                                                              <  1:2859214/150‑1 (MQ=255)
          tCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATtaat                                                                                                                                            >  1:4238833/1‑150 (MQ=255)
           cATTATGACAGGAGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAATAGAAGAAAAAAGCTAAAGTCCCACAAAttaatattaata                                                                                                                                           >  2:2281278/1‑150 (MQ=255)
             ttATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATaaa                                                                                                                                         >  1:4459584/1‑150 (MQ=255)
                 gACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAaca                                                                                                                                     >  2:255992/1‑150 (MQ=255)
                    aGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTc                                                                                                                                  <  1:1949239/150‑1 (MQ=255)
                        gAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCa                                                                                                                              >  1:1470391/1‑150 (MQ=255)
                             aGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAAc                                                                                                                         >  1:3926623/1‑150 (MQ=255)
                                 tCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTgac                                                                                                                     >  1:154020/1‑150 (MQ=255)
                                 tCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTgac                                                                                                                     <  2:2667354/150‑1 (MQ=255)
                                   cAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTgacag                                                                                                                   >  1:2154897/1‑150 (MQ=255)
                                   cAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTgacag                                                                                                                   >  1:548333/1‑150 (MQ=255)
                                    aaTTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCAGCAACTgacagac                                                                                                                 <  1:4844700/150‑1 (MQ=255)
                                     aTTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTgacagac                                                                                                                 <  2:1142761/150‑1 (MQ=255)
                                      ttATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTgacagaca                                                                                                                >  1:1628420/1‑150 (MQ=255)
                                       tATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACaa                                                                                                               <  1:1324042/150‑1 (MQ=255)
                                       tATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACaa                                                                                                               <  1:1324237/150‑1 (MQ=255)
                                        aTGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACaaa                                                                                                              >  1:519970/1‑150 (MQ=255)
                                        aTGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACaaa                                                                                                              >  1:519824/1‑150 (MQ=255)
                                          gAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGa                                                                                                            >  1:1404168/1‑150 (MQ=255)
                                            aTTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGAtt                                                                                                          >  1:4080072/1‑150 (MQ=255)
                                            aTTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGAtt                                                                                                          >  2:3381872/1‑150 (MQ=255)
                                                aaaaaaGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTAtt                                                                                                      >  2:150658/1‑150 (MQ=255)
                                                aaaaaaGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTAtt                                                                                                      >  2:150181/1‑150 (MQ=255)
                                                    aaGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGt                                                                                                  >  1:2092447/1‑150 (MQ=255)
                                                        tttGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATg                                                                                              >  2:87051/1‑150 (MQ=255)
                                                          tgttgCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGaa                                                                                            >  2:3460456/1‑150 (MQ=255)
                                                          tgttgCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGaa                                                                                            >  2:1008826/1‑150 (MQ=255)
                                                               cATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGATGAAAAATGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAAGCaa                                                                                       >  2:2750127/1‑149 (MQ=255)
                                                                atatGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCAct                                                                                      >  1:3542264/1‑150 (MQ=255)
                                                                atatGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCAct                                                                                      <  2:1292527/150‑1 (MQ=255)
                                                                  atGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCActct                                                                                    >  1:1566548/1‑150 (MQ=255)
                                                                  atGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCActct                                                                                    >  2:1766106/1‑150 (MQ=255)
                                                                    gtgtgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAGGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTgg                                                                                  >  2:1946505/1‑150 (MQ=255)
                                                                    gtgtgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTgg                                                                                  >  2:4285646/1‑150 (MQ=255)
                                                                    gtgtgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTgg                                                                                  >  2:1649823/1‑150 (MQ=255)
                                                                    gtgtgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTgg                                                                                  >  2:4706494/1‑150 (MQ=255)
                                                                     tgtgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGc                                                                                 >  2:879739/1‑150 (MQ=255)
                                                                      gtgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGcc                                                                                >  2:924791/1‑150 (MQ=255)
                                                                       tgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCg                                                                               >  1:436331/1‑150 (MQ=255)
                                                                       tgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCg                                                                               >  1:436324/1‑150 (MQ=255)
                                                                         ttATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTa                                                                             <  2:4459584/150‑1 (MQ=255)
                                                                          tATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTATATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTTTTCGGTAATGAAACCACTCTGGCCGTAt                                                                            <  1:1957311/150‑1 (MQ=255)
                                                                             tctcGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACCGACAAAGATTTAGTCGGTAATGAAACCACTCTGGCCGTATCCg                                                                         >  2:1752485/1‑150 (MQ=255)
                                                                             tctcGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCg                                                                         >  2:1751414/1‑150 (MQ=255)
                                                                               tcGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGaa                                                                       >  2:1173728/1‑150 (MQ=255)
                                                                                      ttGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGata                                                                <  2:2877724/150‑1 (MQ=255)
                                                                                      ttGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGata                                                                <  2:713221/150‑1 (MQ=255)
                                                                                       tGCAGGATGTGCCTCTGAATCTTCTATTGATTAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTGGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGatat                                                               <  1:2234823/150‑1 (MQ=255)
                                                                                         cAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTc                                                             >  1:3657521/1‑150 (MQ=255)
                                                                                              tgtgCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCt                                                        <  1:4810360/150‑1 (MQ=255)
                                                                                              tgtgCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCt                                                        >  2:427537/1‑150 (MQ=255)
                                                                                              tgtgCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTGATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCt                                                        >  2:2859606/1‑150 (MQ=255)
                                                                                               gtgCCTCTTAATATTATATTGATTAAAACAATAAAAAAGCTCAAGTCACAAAAATTAATATTAATAAAAACACTCCCCAGCAACTGACAGAAAAAGATTTTTTCGTTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTg                                                       <  1:4180416/150‑1 (MQ=255)
                                                                                                 gCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTgcg                                                     >  1:130538/1‑150 (MQ=255)
                                                                                                    tctGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCtt                                                  <  1:4835993/150‑1 (MQ=255)
                                                                                                           cttctATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGgag                                           >  2:3770581/1‑150 (MQ=255)
                                                                                                           cttctATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGgag                                           >  2:2898331/1‑150 (MQ=255)
                                                                                                            ttctATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATACGAAGAGGATATTCAAGCTGCGCTTGATGgaga                                          >  2:2884409/1‑150 (MQ=255)
                                                                                                                  tgatgaAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGtt                                    >  2:3917394/1‑150 (MQ=255)
                                                                                                                    atgaAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTcc                                  <  2:2125118/150‑1 (MQ=255)
                                                                                                                          agaagaAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTc                            >  1:3219843/1‑150 (MQ=255)
                                                                                                                          agaagaAAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTc                            >  2:1026760/1‑150 (MQ=255)
                                                                                                                            aagaaAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTccc                          <  2:4134471/150‑1 (MQ=255)
                                                                                                                            aagaaAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTccc                          >  2:1279093/1‑150 (MQ=255)
                                                                                                                             agaaAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTcccc                         >  1:3313882/1‑150 (MQ=255)
                                                                                                                              gaaAAAAGCTCAAGTCACACAAAATAATATTAATAAAAACACGCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCt                        <  1:3909245/150‑1 (MQ=255)
                                                                                                                               aaaaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCtt                       >  2:2466793/1‑150 (MQ=255)
                                                                                                                               aaaaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCtt                       >  1:3291691/1‑150 (MQ=255)
                                                                                                                               aaaaaaGCTCAAGTCACAAAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCtt                       <  1:3376900/150‑1 (MQ=255)
                                                                                                                                aaaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCtaa                      >  2:4102112/1‑148 (MQ=255)
                                                                                                                                aaaaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTa                      >  2:4102006/1‑150 (MQ=255)
                                                                                                                                 aaaaGCTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTaa                     >  2:3020233/1‑150 (MQ=255)
                                                                                                                                  aaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGACGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAAt                    >  2:4532023/1‑150 (MQ=255)
                                                                                                                                  aaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAAt                    >  2:2079552/1‑150 (MQ=255)
                                                                                                                                  aaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAAt                    <  1:255992/150‑1 (MQ=255)
                                                                                                                                  aaaGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAAt                    <  1:255984/150‑1 (MQ=255)
                                                                                                                                   aaccTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTTGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAAtt                   <  2:3260816/147‑1 (MQ=255)
                                                                                                                                   aaGCTCAAGTCACAAAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGTTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAAtt                   <  1:150658/150‑1 (MQ=255)
                                                                                                                                    atctaaaGACACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATtc                  <  2:4088692/145‑1 (MQ=255)
                                                                                                                                    aGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATtc                  >  2:3026497/1‑150 (MQ=255)
                                                                                                                                    aGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATtc                  <  2:1404168/150‑1 (MQ=255)
                                                                                                                                      cTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATtctc                >  2:3911382/1‑150 (MQ=255)
                                                                                                                                      cTCAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATtctc                >  2:4597275/1‑150 (MQ=255)
                                                                                                                                        cAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCt              <  2:587795/150‑1 (MQ=255)
                                                                                                                                        cAAGTCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCt              <  2:588196/150‑1 (MQ=255)
                                                                                                                                         aaGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTg             >  2:2334130/1‑150 (MQ=255)
                                                                                                                                          aGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGt            >  2:2521643/1‑150 (MQ=255)
                                                                                                                                            tCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTaa          >  2:571888/1‑150 (MQ=255)
                                                                                                                                            tCACACAAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTaa          >  2:3929268/1‑150 (MQ=255)
                                                                                                                                             cacacaAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAAt         <  2:2770085/150‑1 (MQ=255)
                                                                                                                                              acacaAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATa        >  1:100366/1‑150 (MQ=255)
                                                                                                                                               cacaAAATAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATAc       >  2:4402448/1‑150 (MQ=255)
                                                                                                                                                  aaaGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTg    >  2:336718/1‑150 (MQ=255)
                                                                                                                                                    aGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTTTTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGtt  <  1:4601196/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTTATTCGGTAATGAAACCACTCTGGCCGTATCCGAAGAGGATATTCAAGCTGCGCTTGATGGAGATGAGTTCCGCGTTCCCCTTAATTCTCCTGTAATACTTGTT  >  CP009273/233782‑234079

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: