Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 634,078 T→A 10.7% intergenic (‑49/‑323) dsbG ← / → ahpC thiol:disulfide interchange protein, periplasmic/alkyl hydroperoxide reductase, C22 subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273634,0780TA10.7% 743.8 / 64.5 252intergenic (‑49/‑323)dsbG/ahpCthiol:disulfide interchange protein, periplasmic/alkyl hydroperoxide reductase, C22 subunit
Reads supporting (aligned to +/- strand):  ref base T (132/93);  new base A (13/14);  total (145/107)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 255 total aligned reads displayed.

TGTAATGCCCTGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAGTTCAGATTACAC  >  CP009273/633931‑634224
                                                                                                                                                   |                                                                                                                                                  
tGTAATGCCCTGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCt                                                                                                                                                    <  1:2602022/148‑1 (MQ=255)
  tAATGCCCTGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTAt                                                                                                                                                  <  2:14997/148‑1 (MQ=255)
   aaTGCCCTGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAATGGCATAACCTATc                                                                                                                                                 >  1:2858153/1‑148 (MQ=255)
     tGCCCTGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCAc                                                                                                                                               >  2:932612/1‑148 (MQ=255)
        ccTGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGt                                                                                                                                            <  2:1583123/148‑1 (MQ=255)
          tGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCa                                                                                                                                          <  2:976499/148‑1 (MQ=255)
          tGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCa                                                                                                                                          <  2:2222366/148‑1 (MQ=255)
           gTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCAt                                                                                                                                         >  1:1530344/1‑148 (MQ=255)
            tttttCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATa                                                                                                                                        <  2:1187205/148‑1 (MQ=255)
                  aaTCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAag                                                                                                                                  <  1:597676/148‑1 (MQ=255)
                   aTCGCTTGTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAaga                                                                                                                                 <  2:225576/148‑1 (MQ=255)
                     cGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGc                                                                                                                               >  2:16098/1‑148 (MQ=255)
                     cGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGc                                                                                                                               >  1:323081/1‑148 (MQ=255)
                      gCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCt                                                                                                                              >  2:3023554/1‑148 (MQ=255)
                      gCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCt                                                                                                                              <  1:931885/148‑1 (MQ=255)
                        ttttACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTa                                                                                                                            <  2:3297923/148‑1 (MQ=255)
                        ttttACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTa                                                                                                                            >  1:2190688/1‑148 (MQ=255)
                         tttACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAg                                                                                                                           >  2:2055421/1‑148 (MQ=255)
                          ttACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGa                                                                                                                          <  1:1661965/148‑1 (MQ=255)
                             cTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCa                                                                                                                       <  1:1700007/148‑1 (MQ=255)
                               ggAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAgg                                                                                                                     <  2:459779/148‑1 (MQ=255)
                                gAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGt                                                                                                                    >  2:2194623/1‑148 (MQ=255)
                                gAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGt                                                                                                                    >  2:609328/1‑148 (MQ=255)
                                 aGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTg                                                                                                                   <  1:819103/148‑1 (MQ=255)
                                   cAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGAt                                                                                                                 >  1:3030434/1‑148 (MQ=255)
                                   cAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGAt                                                                                                                 >  1:35104/1‑148 (MQ=255)
                                   cAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGAt                                                                                                                 <  1:3334706/148‑1 (MQ=255)
                                     ggAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTg                                                                                                               <  1:1133896/148‑1 (MQ=255)
                                      gAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGc                                                                                                              >  2:342701/1‑148 (MQ=255)
                                             ctctGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGt                                                                                                       <  1:751013/148‑1 (MQ=255)
                                             ctctGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGt                                                                                                       <  1:1213956/148‑1 (MQ=255)
                                             ctctGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTATTAACATTAAATTGTACTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATAACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGt                                                                                                       <  1:1001878/148‑1 (MQ=255)
                                              tctGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCgttgtt                                                                                                      >  1:1911083/1‑148 (MQ=255)
                                               ctGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGttt                                                                                                     <  2:2866996/148‑1 (MQ=255)
                                                tGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTa                                                                                                    <  2:1729848/148‑1 (MQ=255)
                                                 gCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTAt                                                                                                   <  2:1219880/148‑1 (MQ=255)
                                                  cGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATg                                                                                                  >  2:299067/1‑148 (MQ=255)
                                                   gAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGa                                                                                                 <  2:686105/148‑1 (MQ=255)
                                                     aGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAgg                                                                                               >  2:1548441/1‑148 (MQ=255)
                                                     aGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAgg                                                                                               >  2:406361/1‑148 (MQ=255)
                                                      ggCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAggg                                                                                              >  1:835756/1‑148 (MQ=255)
                                                           ttGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGtgttg                                                                                         >  2:1284254/1‑148 (MQ=255)
                                                             gCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTa                                                                                       <  1:2496385/148‑1 (MQ=255)
                                                                ggAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATc                                                                                    <  1:3075243/148‑1 (MQ=255)
                                                                ggAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATc                                                                                    <  1:3074627/148‑1 (MQ=255)
                                                                ggAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATc                                                                                    <  1:495214/148‑1 (MQ=255)
                                                                ggAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATc                                                                                    <  2:2929836/148‑1 (MQ=255)
                                                                    gCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATg                                                                                >  1:1994191/1‑148 (MQ=255)
                                                                      agagCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCAATGTc                                                                              >  1:1368153/1‑148 (MQ=255)
                                                                         gCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCgtt                                                                           >  2:144646/1‑148 (MQ=255)
                                                                           cAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCgttgt                                                                         >  1:3399370/1‑148 (MQ=255)
                                                                                aaaGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCAt                                                                    >  2:793773/1‑148 (MQ=255)
                                                                                 aaGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCAtt                                                                   <  2:3539157/148‑1 (MQ=255)
                                                                                 aaGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCAtt                                                                   >  1:739355/1‑148 (MQ=255)
                                                                                 aaGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCAtt                                                                   >  1:29651/1‑148 (MQ=255)
                                                                                 aaGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCAtt                                                                   >  1:3249358/1‑148 (MQ=255)
                                                                                   gTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTg                                                                 >  1:1575598/1‑148 (MQ=255)
                                                                                    tATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGt                                                                >  2:2404393/1‑148 (MQ=255)
                                                                                    tATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGAAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGACCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGt                                                                <  1:3553831/148‑1 (MQ=255)
                                                                                      tCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTTTTGCATTTGTaa                                                              >  1:1828310/1‑148 (MQ=255)
                                                                                      tCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTaa                                                              >  1:3254605/1‑148 (MQ=255)
                                                                                        ttttttAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAgg                                                            >  1:1778289/1‑148 (MQ=255)
                                                                                         tttttAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAggg                                                           >  2:552026/1‑148 (MQ=255)
                                                                                         tttttAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAggg                                                           >  2:1772515/1‑148 (MQ=255)
                                                                                          ttttAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGc                                                          <  1:1457251/148‑1 (MQ=255)
                                                                                             tAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAac                                                       >  1:1285932/1‑148 (MQ=255)
                                                                                               aCATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAacac                                                     <  1:1603497/148‑1 (MQ=255)
                                                                                               aCATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAacac                                                     >  1:2214131/1‑148 (MQ=255)
                                                                                               aCATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAacac                                                     <  1:3364698/148‑1 (MQ=255)
                                                                                               aCATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAacac                                                     >  1:45959/1‑148 (MQ=255)
                                                                                               aCATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTGTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAacac                                                     >  1:45990/1‑148 (MQ=255)
                                                                                                 aTTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCt                                                   <  1:1706981/148‑1 (MQ=255)
                                                                                                 aTTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCt                                                   >  1:1982962/1‑148 (MQ=255)
                                                                                                 aTTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCt                                                   >  2:2603510/1‑148 (MQ=255)
                                                                                                  ttAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTc                                                  <  1:1248807/148‑1 (MQ=255)
                                                                                                  ttAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTc                                                  <  2:3117543/148‑1 (MQ=255)
                                                                                                   tAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCa                                                 >  2:2819034/1‑148 (MQ=255)
                                                                                                   tAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCa                                                 >  2:1866368/1‑148 (MQ=255)
                                                                                                     aTTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGc                                               <  1:1962088/148‑1 (MQ=255)
                                                                                                     aTTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGc                                               >  1:2007815/1‑148 (MQ=255)
                                                                                                      tttGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGcc                                              <  2:2337867/148‑1 (MQ=255)
                                                                                                             ttttCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTgcaggc                                       >  1:2852617/1‑148 (MQ=255)
                                                                                                                  agtcagtGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAggcaggca                                  >  2:2796575/1‑148 (MQ=255)
                                                                                                                  agtcagtGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAggcaggca                                  <  2:1983036/148‑1 (MQ=255)
                                                                                                                   gtcagtGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTTTAAGGGCAACACCTAAGCCTTCAGGCAGGCAc                                 >  1:2388407/1‑148 (MQ=255)
                                                                                                                    tcagtGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACt                                >  2:3225148/1‑148 (MQ=255)
                                                                                                                      agtGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGa                              <  2:2005945/148‑1 (MQ=255)
                                                                                                                       gtGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGaa                             >  2:232775/1‑148 (MQ=255)
                                                                                                                       gtGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGaa                             >  2:231592/1‑148 (MQ=255)
                                                                                                                               gTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCaa                     >  2:880506/1‑148 (MQ=255)
                                                                                                                                  gAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAgg                  >  1:2424968/1‑148 (MQ=255)
                                                                                                                                     tAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACAAAAGGGTa               >  2:2837295/1‑148 (MQ=255)
                                                                                                                                     tAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTa               <  1:2108960/148‑1 (MQ=255)
                                                                                                                                      aaaaGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAg              <  1:1875454/148‑1 (MQ=255)
                                                                                                                                      aaaaGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAg              >  1:3007875/1‑148 (MQ=255)
                                                                                                                                       aaaGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAGt             <  2:1916441/148‑1 (MQ=255)
                                                                                                                                         aGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAGTTc           >  1:857613/1‑148 (MQ=255)
                                                                                                                                          ggCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAGTTCa          >  2:1472058/1‑148 (MQ=255)
                                                                                                                                                 ccTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTGTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAGTTCAGATTaca   >  1:311285/1‑148 (MQ=255)
                                                                                                                                                  cTATCACGGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAGTTCAGATTacac  >  2:1029009/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGTAATGCCCTGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAGTTCAGATTACAC  >  CP009273/633931‑634224

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: