Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,681,873 T→A 89.7% Q97L (CAG→CTG)  glrR ← response regulator regulating glmY sRNA in two‑component system with sensor protein GlrK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,681,8730TA89.7% 849.8 / 69.8 253Q97L (CAG→CTG) glrRresponse regulator regulating glmY sRNA in two‑component system with sensor protein GlrK
Reads supporting (aligned to +/- strand):  ref base T (11/15);  new base A (90/137);  total (101/152)
Fisher's exact test for biased strand distribution p-value = 8.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 257 total aligned reads displayed.

AACGCAGCATCAGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGACTAAATCTACTTT  >  CP009273/2681727‑2682018
                                                                                                                                                  |                                                                                                                                                 
     agcaTCCGCGGGCTGCGGGTGACACTTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCACTTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAACACGCCCTGCTGTGTTgca                                                                                                                                           <  2:3825095/150‑1 (MQ=255)
     agcaTCAGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTgca                                                                                                                                           <  2:3384815/150‑1 (MQ=255)
       caTCAGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTgcagc                                                                                                                                         >  2:409026/1‑150 (MQ=255)
         tCAGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGAGTTGCAGCaa                                                                                                                                       <  2:2253464/150‑1 (MQ=255)
          cAGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAAc                                                                                                                                      >  1:840717/1‑150 (MQ=255)
          cAGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAAc                                                                                                                                      <  1:3587336/150‑1 (MQ=255)
            gCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGTAAACTAAAAACGCCCTGCAGTGTTGCAGCAACgg                                                                                                                                    <  1:419260/150‑1 (MQ=255)
               ggCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCAt                                                                                                                                 >  2:1723776/1‑150 (MQ=255)
               ggCTGCGGGTGACAAGGGCCTCGCGCCAGCGTTCTTCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGCAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCAt                                                                                                                                 <  1:3455723/150‑1 (MQ=255)
                gCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATc                                                                                                                                <  1:560255/150‑1 (MQ=255)
                   gCGGGTGAAAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATGGCCTTTTATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCggg                                                                                                                             <  1:926616/150‑1 (MQ=255)
                    cGGGTGACAAGTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGa                                                                                                                            <  2:438606/150‑1 (MQ=255)
                     gggTGACAAGTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGaa                                                                                                                           <  2:4174631/150‑1 (MQ=255)
                        tGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAg                                                                                                                        <  2:4093393/150‑1 (MQ=255)
                         gACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGa                                                                                                                       <  2:2941336/150‑1 (MQ=255)
                         gACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGa                                                                                                                       <  2:679269/150‑1 (MQ=255)
                           cAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAAc                                                                                                                     <  2:4108898/150‑1 (MQ=255)
                           cAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAAc                                                                                                                     <  2:1830331/150‑1 (MQ=255)
                             aTTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCa                                                                                                                   >  1:3893372/1‑150 (MQ=255)
                              ttGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCAt                                                                                                                  >  1:4102311/1‑150 (MQ=255)
                               tGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTAAATTGCCTGATATAGCGCGTCTTTTTCGACGGTCTTGGTGAGTAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATg                                                                                                                 <  1:1528588/150‑1 (MQ=255)
                                  cTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATgcgc                                                                                                              >  2:2331169/1‑150 (MQ=255)
                                       gcCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGtga                                                                                                         <  1:3741196/150‑3 (MQ=255)
                                       gcCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTaa                                                                                                         <  2:2208627/150‑1 (MQ=255)
                                       gcCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTaa                                                                                                         >  2:1215091/1‑150 (MQ=255)
                                       gcCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTaa                                                                                                         >  1:1211010/1‑150 (MQ=255)
                                          aGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGaa                                                                                                      >  2:463892/1‑150 (MQ=255)
                                          aGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGac                                                                                                      >  1:310120/1‑149 (MQ=255)
                                          aGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGaa                                                                                                      >  1:1992590/1‑150 (MQ=255)
                                              ttCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAACACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGaataat                                                                                                  <  2:3931190/150‑1 (MQ=255)
                                              ttCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGaataat                                                                                                  <  2:820617/150‑1 (MQ=255)
                                                cATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTa                                                                                                <  2:3725578/150‑1 (MQ=255)
                                                cATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTa                                                                                                >  2:1669802/1‑150 (MQ=255)
                                                 aTCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTAc                                                                                               >  2:3635212/1‑150 (MQ=255)
                                                  tCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTAcc                                                                                              <  2:1776635/150‑1 (MQ=255)
                                                  tCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTAcc                                                                                              >  1:493262/1‑150 (MQ=255)
                                                    ggTGGCTGGCGCGGATTGCTCCAGCGCATCGTAAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCgg                                                                                            <  1:3469564/150‑1 (MQ=255)
                                                        gctggcGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCAtt                                                                                        <  1:2522220/150‑1 (MQ=255)
                                                         ctggcGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTc                                                                                       <  1:4172254/150‑1 (MQ=255)
                                                         ctggcGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTc                                                                                       <  2:4314606/150‑1 (MQ=255)
                                                            gcgcGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCTGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCg                                                                                    <  1:147556/150‑1 (MQ=255)
                                                            gcgcGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCg                                                                                    <  1:2897786/150‑1 (MQ=255)
                                                                ggATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTg                                                                                <  1:2033185/150‑1 (MQ=255)
                                                                ggATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTg                                                                                <  2:959611/150‑1 (MQ=255)
                                                                  aTTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCa                                                                              >  2:1467655/1‑150 (MQ=255)
                                                                    tGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACt                                                                            >  2:4202305/1‑150 (MQ=255)
                                                                       tCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTc                                                                         >  1:3614312/1‑150 (MQ=255)
                                                                       tCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTc                                                                         <  1:2302053/150‑1 (MQ=255)
                                                                        ccAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCt                                                                        <  2:2783344/150‑1 (MQ=255)
                                                                          aGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTgg                                                                      >  1:2733388/1‑150 (MQ=255)
                                                                           gctcATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGATGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTTAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGa                                                                     <  2:253887/147‑1 (MQ=255)
                                                                              cATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGAttt                                                                  <  2:3563313/150‑1 (MQ=255)
                                                                                  gTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGc                                                              >  2:3962566/1‑150 (MQ=255)
                                                                                      aTTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGACTAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAAc                                                          >  1:4479182/1‑150 (MQ=255)
                                                                                      aTTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAAc                                                          <  1:4054885/150‑1 (MQ=255)
                                                                                           cTGATATAGCGCGTCTTTGTCGACCGGCTTGGTGAGGAAACTAAACACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTg                                                     <  1:1143398/150‑1 (MQ=255)
                                                                                           cTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTg                                                     >  2:2117441/1‑150 (MQ=255)
                                                                                              atatAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCAt                                                  <  2:542927/150‑1 (MQ=255)
                                                                                               tataGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTGTCTGGATTTCAGCAAACAGCTGCATa                                                 >  1:672185/1‑150 (MQ=255)
                                                                                                ataGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATAc                                                >  2:253278/1‑150 (MQ=255)
                                                                                                ataGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATAc                                                <  2:2510613/150‑1 (MQ=255)
                                                                                                ataGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATAc                                                >  1:756787/1‑150 (MQ=255)
                                                                                                 taGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATAcc                                               >  2:4249321/1‑150 (MQ=255)
                                                                                                  aGCGCGTCTTTGTCGACCGGCTTGGTGAGGAAAATAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCg                                              <  1:1341308/150‑1 (MQ=255)
                                                                                                  aGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCg                                              <  2:1852021/150‑1 (MQ=255)
                                                                                                    cgcgTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTc                                            >  1:1096944/1‑150 (MQ=255)
                                                                                                     gcgTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTcc                                           <  2:2641217/150‑1 (MQ=255)
                                                                                                      cgTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCa                                          <  2:732430/150‑1 (MQ=255)
                                                                                                      cgTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCa                                          >  2:3874961/1‑150 (MQ=255)
                                                                                                       gTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCAt                                         <  2:1644822/150‑1 (MQ=255)
                                                                                                        tCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCAtt                                        <  2:782993/150‑1 (MQ=255)
                                                                                                         cTTTGTCGACAGGCTTGGTGATGAAACTAAAATCGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCAttt                                       <  2:840604/150‑1 (MQ=255)
                                                                                                          tttGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTc                                      >  1:1028808/1‑150 (MQ=255)
                                                                                                           ttGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTca                                     <  2:2035818/150‑1 (MQ=255)
                                                                                                                gACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTcatccat                                <  2:1923093/150‑1 (MQ=255)
                                                                                                                 aCAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTcatccatc                               >  1:3929239/1‑150 (MQ=255)
                                                                                                                   aGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCg                             >  1:1015897/1‑150 (MQ=255)
                                                                                                                   aGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCg                             >  1:3373518/1‑150 (MQ=255)
                                                                                                                   aGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCg                             >  2:1376510/1‑150 (MQ=255)
                                                                                                                    ggCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGc                            >  1:3182374/1‑150 (MQ=255)
                                                                                                                       ttGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAgg                         <  1:1810152/150‑1 (MQ=255)
                                                                                                                       ttGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAgg                         <  1:2410428/150‑1 (MQ=255)
                                                                                                                         ggTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTc                       <  1:167069/150‑1 (MQ=255)
                                                                                                                         ggTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTc                       <  1:1576545/150‑1 (MQ=255)
                                                                                                                         ggTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTc                       <  2:1296112/150‑1 (MQ=255)
                                                                                                                          gTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCg                      >  1:1702991/1‑150 (MQ=255)
                                                                                                                            gAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCt                    <  1:4146783/150‑1 (MQ=255)
                                                                                                                            gAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCt                    >  1:6622/1‑150 (MQ=255)
                                                                                                                             aGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCtg                   <  1:2421649/150‑1 (MQ=255)
                                                                                                                             aGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCtg                   <  1:2100334/150‑1 (MQ=255)
                                                                                                                               taaaCTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCtgat                 <  1:1225523/149‑1 (MQ=255)
                                                                                                                               gAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCtgat                 <  1:1166855/150‑1 (MQ=255)
                                                                                                                                aaaCTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCtgatg                <  1:3762524/150‑1 (MQ=255)
                                                                                                                                aaaCTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAAGTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCtgatg                <  1:2169644/150‑1 (MQ=255)
                                                                                                                                 aaCTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACGTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCtgatga               >  2:34314/1‑150 (MQ=255)
                                                                                                                                    tAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGACTa            >  2:1058915/1‑150 (MQ=255)
                                                                                                                                     aaaaaCGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGACTaa           >  2:4360819/1‑150 (MQ=255)
                                                                                                                                      aaaaCGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGACTaaa          >  2:1173855/1‑150 (MQ=255)
                                                                                                                                          cGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGACTAAATCTa      >  1:4291907/1‑150 (MQ=255)
                                                                                                                                             ccTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGACTAAATCTACtt   >  1:2869532/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AACGCAGCATCAGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGACTAAATCTACTTT  >  CP009273/2681727‑2682018

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: