Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,904,977 A→T 7.0% intergenic (‑92/+222) pstS ← / ← glmS periplasmic phosphate binding protein, high‑affinity/L‑glutamine:D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,904,9770AT 7.0% 831.8 / 41.4 271intergenic (‑92/+222)pstS/glmSperiplasmic phosphate binding protein, high‑affinity/L‑glutamine:D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (99/152);  new base T (6/13);  total (105/166)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 272 total aligned reads displayed.

CTCATCGATAAGGTCGCGGCGACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTC  >  CP009273/3904830‑3905122
                                                                                                                                                   |                                                                                                                                                 
cTCATCGATAAGGTCGCGGCGACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGa                                                                                                                                                   <  2:738348/148‑1 (MQ=255)
 tCATCGATAAGGTCGCGGCGACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGAc                                                                                                                                                  <  2:350314/148‑1 (MQ=255)
 tCATCGATAAGGTCGCGGCGACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGAc                                                                                                                                                  <  2:3018464/148‑1 (MQ=255)
  cATCGATAAGGTCGCGGCGACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACa                                                                                                                                                 <  1:163019/148‑1 (MQ=255)
    tCGAGAAGGTCGCGGCGACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGt                                                                                                                                               <  1:2476956/148‑1 (MQ=255)
     cGATAAGGTCGCGGCGACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGtt                                                                                                                                              >  2:1287583/1‑148 (MQ=255)
      gATAAGGTCGCGGCGACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGttt                                                                                                                                             <  2:3163874/148‑1 (MQ=255)
       aTAAGGTCGCGGCGACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGTCAGtttt                                                                                                                                            <  2:661465/148‑1 (MQ=255)
       aTAAGGTCGCGGCGACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGtttt                                                                                                                                            <  1:2398410/148‑1 (MQ=255)
        tAAGGTCGCGGCGACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGTCAGTTTTa                                                                                                                                           >  2:3419673/1‑148 (MQ=255)
        tAAGGTCGCGGCGACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTa                                                                                                                                           <  1:804540/148‑1 (MQ=255)
        tAAGGTCGCGGCGACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTa                                                                                                                                           <  1:2051159/148‑1 (MQ=255)
            gTCGCGGCGACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGAc                                                                                                                                       >  1:214158/1‑148 (MQ=255)
            gTCGCGGCGACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGAc                                                                                                                                       <  1:2303002/148‑1 (MQ=255)
                 ggcgACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACagaga                                                                                                                                  >  1:1935339/1‑148 (MQ=255)
                  gcgACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACagagag                                                                                                                                 <  1:711134/148‑1 (MQ=255)
                   cgACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAgagaga                                                                                                                                <  1:1327009/148‑1 (MQ=255)
                     acaacaGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATa                                                                                                                              >  2:447560/1‑148 (MQ=255)
                      caacaGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATaa                                                                                                                             >  1:3342964/1‑148 (MQ=255)
                       aacaGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATaaa                                                                                                                            <  2:1085070/148‑1 (MQ=255)
                            ttGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCtt                                                                                                                       >  2:2221085/1‑148 (MQ=255)
                                gACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGt                                                                                                                   <  1:1595194/148‑1 (MQ=255)
                                gACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGt                                                                                                                   <  2:2626211/148‑1 (MQ=255)
                                 aCGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTc                                                                                                                  >  2:3125625/1‑148 (MQ=255)
                                  cGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCt                                                                                                                 >  2:3516496/1‑148 (MQ=255)
                                    gtggtACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGa                                                                                                               <  1:1580318/148‑1 (MQ=255)
                                       gtACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGAttt                                                                                                            >  2:2436159/1‑148 (MQ=255)
                                         aCGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTaa                                                                                                          <  1:2554837/148‑1 (MQ=255)
                                           gCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAAt                                                                                                        >  1:3344827/1‑148 (MQ=255)
                                           gCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAAt                                                                                                        >  1:3181272/1‑148 (MQ=255)
                                            cATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATa                                                                                                       <  1:740276/148‑1 (MQ=255)
                                              tAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGTCAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAg                                                                                                     <  2:374850/148‑1 (MQ=255)
                                                aCTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCg                                                                                                   <  1:3302825/148‑1 (MQ=255)
                                                aCTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCg                                                                                                   <  2:3447911/148‑1 (MQ=255)
                                                aCTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCg                                                                                                   <  2:1825026/148‑1 (MQ=255)
                                                 cTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGt                                                                                                  <  1:833844/148‑1 (MQ=255)
                                                   ttCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTg                                                                                                >  2:2919597/1‑148 (MQ=255)
                                                      aTAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGata                                                                                             >  1:1265897/1‑148 (MQ=255)
                                                       tAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGatat                                                                                            <  2:793244/148‑1 (MQ=255)
                                                       tAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGatat                                                                                            <  2:2607635/148‑1 (MQ=255)
                                                        aaTGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATAtt                                                                                           >  2:1782686/1‑148 (MQ=255)
                                                             ctcCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAATTTGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTc                                                                                      <  2:295501/148‑1 (MQ=255)
                                                               ccTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCaa                                                                                    <  1:732520/148‑1 (MQ=255)
                                                                cTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAAt                                                                                   <  2:1952070/148‑1 (MQ=255)
                                                                     aggaTTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACaa                                                                              >  2:1174986/1‑148 (MQ=255)
                                                                       gaTTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAAc                                                                            >  1:994284/1‑148 (MQ=255)
                                                                       gaTTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAAc                                                                            >  2:2456838/1‑148 (MQ=255)
                                                                        aTTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTATGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACa                                                                           >  2:1948472/1‑148 (MQ=255)
                                                                         ttCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACAt                                                                          <  2:1626015/148‑1 (MQ=255)
                                                                          tCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACAtt                                                                         >  1:606692/1‑148 (MQ=255)
                                                                          tCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACAtt                                                                         >  2:520424/1‑148 (MQ=255)
                                                                            aTAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATtaa                                                                       <  2:2227958/148‑1 (MQ=255)
                                                                            aTAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATtaa                                                                       >  1:3422932/1‑148 (MQ=255)
                                                                             tAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTAAATTACAAACATtaat                                                                      <  1:1020296/148‑1 (MQ=255)
                                                                             tAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTAAATTACAAACATtaat                                                                      <  1:2864350/148‑1 (MQ=255)
                                                                              aaaGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATtaata                                                                     <  1:579445/148‑1 (MQ=255)
                                                                              aaaGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATtaata                                                                     >  1:1327934/1‑148 (MQ=255)
                                                                              aaaGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATtaata                                                                     <  2:1982281/148‑1 (MQ=255)
                                                                               aaGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATtaataa                                                                    >  1:1058875/1‑148 (MQ=255)
                                                                                 gCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACg                                                                  <  2:3087723/148‑1 (MQ=255)
                                                                                   aTTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGaa                                                                >  1:3376760/1‑148 (MQ=255)
                                                                                      gttttttGTATAAGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATTTAATGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAagag                                                             <  2:3181272/148‑1 (MQ=255)
                                                                                      gtttgttGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAagag                                                             <  2:3344827/148‑1 (MQ=255)
                                                                                       tttgttGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAgaga                                                            <  1:2039078/148‑1 (MQ=255)
                                                                                           ttgGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGAc                                                        >  2:2541914/1‑148 (MQ=255)
                                                                                            tgGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACa                                                       >  1:1498601/1‑148 (MQ=255)
                                                                                             ggCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAg                                                      <  2:144066/148‑1 (MQ=255)
                                                                                                tACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGTCAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGaaa                                                   <  1:2345063/148‑1 (MQ=255)
                                                                                                  cGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGaaaaa                                                 <  2:3342835/148‑1 (MQ=255)
                                                                                                  cGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGaaaaa                                                 <  2:86532/148‑1 (MQ=255)
                                                                                                      aaGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGTCAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAAtttt                                             >  2:1848907/1‑148 (MQ=255)
                                                                                                      aaGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAAtttt                                             <  2:109465/148‑1 (MQ=255)
                                                                                                        gCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCa                                           >  2:3435171/1‑148 (MQ=255)
                                                                                                         cAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCAt                                          <  2:1434042/148‑1 (MQ=255)
                                                                                                          aaaaTAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCAtt                                         <  2:2319830/148‑1 (MQ=255)
                                                                                                          aaaaTAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCAtt                                         <  1:3030700/148‑1 (MQ=255)
                                                                                                             aTAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCtg                                      >  2:2674239/1‑148 (MQ=255)
                                                                                                                 gACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGAc                                  <  1:183394/148‑1 (MQ=255)
                                                                                                                 gACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGAc                                  >  1:186415/1‑148 (MQ=255)
                                                                                                                    aaaCAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACaga                               <  1:1780041/148‑1 (MQ=255)
                                                                                                                    aaaCAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACaga                               <  1:2726150/148‑1 (MQ=255)
                                                                                                                     aaCAGGTGACAGTTATATGTAAGGAATATGTCAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACagag                              >  1:865666/1‑148 (MQ=255)
                                                                                                                     aaCAGGTGACAGTTATATGTAAGGAATATGTCAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACagag                              <  2:642376/148‑1 (MQ=255)
                                                                                                                         ggTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGaaaa                          <  1:1705099/148‑1 (MQ=255)
                                                                                                                         ggTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGaaaa                          <  1:294796/148‑1 (MQ=255)
                                                                                                                          gTGACAGTTATATGTAAGGAATATGTCAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGaaaaa                         <  1:3051963/148‑1 (MQ=255)
                                                                                                                          gTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGaaaaa                         <  2:1561868/148‑1 (MQ=255)
                                                                                                                            gACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGt                       <  2:3462956/148‑1 (MQ=255)
                                                                                                                               aGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGc                    <  2:2042598/148‑1 (MQ=255)
                                                                                                                                 ttATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCg                  <  1:3217112/148‑1 (MQ=255)
                                                                                                                                 ttATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCg                  >  1:521457/1‑148 (MQ=255)
                                                                                                                                  tataTGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGa                 >  1:391001/1‑148 (MQ=255)
                                                                                                                                     atGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGa              <  1:2436159/148‑1 (MQ=255)
                                                                                                                                       gTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATTACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATg            <  2:796752/148‑1 (MQ=255)
                                                                                                                                       gTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATTACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATg            <  2:797019/148‑1 (MQ=255)
                                                                                                                                          aGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACg         >  1:3493891/1‑148 (MQ=255)
                                                                                                                                          aGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACg         >  1:3508435/1‑148 (MQ=255)
                                                                                                                                           ggAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACgg        >  2:532720/1‑148 (MQ=255)
                                                                                                                                              atatGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGttt     >  1:877145/1‑148 (MQ=255)
                                                                                                                                                 tGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAATGTAGCCGAAGATGACGGTTTGTc  >  2:2402330/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTCATCGATAAGGTCGCGGCGACAACAGTTGCGACGGTGGTACGCATAACTTTCATAATGTCTCCTGGGAGGATTCATAAAGCATTGTTTGTTGGCTACGAGAAGCAAAATAGGACAAACAGGTGACAGTTATATGTAAGGAATATGACAGTTTTATGACAGAGAGATAAAGTCTTCAGTCTGATTTAAATAAGCGTTGATATTCAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTC  >  CP009273/3904830‑3905122

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