Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,385,578 A→T 6.7% intergenic (‑191/+239) tldD ← / ← yhdP putative peptidase/DUF3971‑AsmA2 domains protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,385,5780AT 6.7% 942.0 / 3.8 335intergenic (‑191/+239)tldD/yhdPputative peptidase/DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base A (193/96);  new base T (6/18);  total (199/136)
Fisher's exact test for biased strand distribution p-value = 8.50e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 340 total aligned reads displayed.

CGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACG  >  CP009273/3385430‑3385727
                                                                                                                                                    |                                                                                                                                                     
cGGCTGCCGGTTGCGGCGTGAACGCCATAGCCAGCCTACGGTTATGTTCGGTTTTGTCGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAg                                                                                                                                                      <  1:5423759/150‑1 (MQ=14)
cGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGGTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGCGGCGTTAACGCCTTATCCGACCTACGGTTATGTTCGGTGTGTAg                                                                                                                                                      >  1:2792486/1‑150 (MQ=14)
 ggCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAgg                                                                                                                                                     >  2:4052341/1‑150 (MQ=255)
   cTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGcc                                                                                                                                                   >  1:3243389/1‑150 (MQ=255)
    tGCCGGTGGCGGCGTGAACGCCCTAGCCAGCCTAGGGTTATGTTCCGGTTTGTTGGCCTGATAAGACGCACCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCt                                                                                                                                                  <  1:4527812/150‑1 (MQ=17)
    tGCCGGATGCGGCGTGACCGCCCAACCCAGCCTACGGTTATGTTCCGTTTTTTCGGCCTTATACGACGCACCGCGTCCCAACCGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCt                                                                                                                                                  <  1:2667212/150‑1 (MQ=14)
     gcggtttGCGGCGTTAACGCCCTACCCCGCCTAGGGTTAGGTTCGGGTTTGTCGGCCTTATACGCCGCACCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTg                                                                                                                                                 <  1:1342241/144‑1 (MQ=14)
      ccGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGa                                                                                                                                                >  2:3279481/1‑150 (MQ=255)
        ggATGCGGCGTGAACGCCCTATCCCGCCTACGGTTATGTTCAGGGTAGTAGGCCTGATAACACGAAAACCCACGAATCACCAACCGTCACACGAATCCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTGTGTAGGCCTGaga                                                                                                                                              >  1:5275151/1‑148 (MQ=9)
        ggATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATa                                                                                                                                              >  1:3119309/1‑150 (MQ=255)
         gATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATaa                                                                                                                                             <  2:5030156/150‑1 (MQ=255)
         gATGCGGCGTGAACGCCCTAGCCAGCCTACGGTTATGTACCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTCTCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATaa                                                                                                                                             <  2:2268311/150‑1 (MQ=255)
          agggggcgTGAACGCCCTACCCCGCCTAGGGTTTTGTTCGGGTTGGTTGGCCTTATAAGACGCACCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAg                                                                                                                                            <  2:2322278/146‑1 (MQ=14)
            gcggggTTAACGCCCTATCCAGCCAAGGTTTTTGTTCGGGTTTGTCGGCCTTATAAGACGCACCGCGTCGACTCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGAc                                                                                                                                          <  1:2164440/150‑1 (MQ=14)
            gcggcgTGAACGCCATAACCAGCCTAGGTTTATGTTCCGTTTTTTTGGCCTTATAATACCCACCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGAc                                                                                                                                          <  1:1177276/150‑1 (MQ=14)
               gcgTGACCGCCCTATCCAGCCTACGGTTTTGTTCCGGTTTGTAGGCCTTATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGca                                                                                                                                       <  1:3620643/150‑1 (MQ=255)
               gcgTGAACGCCCTACCCAGCCTACGTTTTTGTTCCGTTTTGTAGGCCTTAAAAGACGCACCGCGTCGCATCCGCCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGca                                                                                                                                       <  1:855664/150‑1 (MQ=14)
                 gTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcaca                                                                                                                                     >  2:3506434/1‑150 (MQ=255)
                  tGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAg                                                                                                                                    <  2:4137242/150‑1 (MQ=255)
                   gAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGc                                                                                                                                   >  2:4646035/1‑150 (MQ=255)
                   gAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCATATCCGACCTACGGTTATGTTCAGTGTGTAGGCCAGATAAGAAGCACAGc                                                                                                                                   >  1:2311265/1‑150 (MQ=14)
                    aaCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCg                                                                                                                                  <  2:4005325/150‑1 (MQ=255)
                    aaCGCCCAACCCAGCCTACGGTTATGTTCCGTTTTGTAGGCCTTATAAGACGCACAGCGTCGCATCAGCCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCg                                                                                                                                  <  1:4084246/150‑1 (MQ=14)
                     aCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGt                                                                                                                                 >  1:424322/1‑150 (MQ=255)
                      cGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                                                                                                                                <  2:3916676/150‑1 (MQ=255)
                       gCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTCTCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCg                                                                                                                               >  1:4441816/1‑150 (MQ=255)
                       gCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCATATCCGACCTACGGTTATGTTCAGTATGTATGCCAGATAAGACGCACAgagacg                                                                                                                               >  1:4179763/1‑145 (MQ=14)
                          cTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTGTGTAGGCCTGATAAGACGCACAGCGTCGCAt                                                                                                                            >  1:5430809/1‑150 (MQ=255)
                           tATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATc                                                                                                                           >  2:3387500/1‑150 (MQ=255)
                            ctCCAGCCTAGGGTTATGTTCGGGTTTGTAGGCCTGATAAGCCGCACAGCGTCGCAACAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCa                                                                                                                          <  2:607147/149‑1 (MQ=255)
                            aTCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCa                                                                                                                          >  2:2288349/1‑150 (MQ=255)
                             tCCCGCCTACGTTTTTGTTCCGGTTTTTAGCCCTTAAAAGACCCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg                                                                                                                         <  1:5144060/150‑1 (MQ=11)
                              cccgCCTAGGTTTTTGTTCGGGTTTGTGGGCGTGATACGACGCCCCGCGTCGCATCCGGCCACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAgg                                                                                                                        <  1:5186991/147‑1 (MQ=11)
                              ccAGCCTTGGTTTTTGTACCGTTTTTTAGGCCTTATAAGACCCACAGCGTCTAATCAGCCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAgg                                                                                                                        <  1:2415928/150‑1 (MQ=11)
                                aGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCa                                                                                                                      <  1:3631569/150‑1 (MQ=38)
                                      cGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCt                                                                                                                >  1:3742973/1‑150 (MQ=37)
                                          tttGTTCGGTTTTTTTGGCCTTATATGACGCACTGGGTCGCATATGGAAAGGGCTGTCGTAGGCGGCTTAAACGCCTTATCCGACCTACGGTTTTGTTCCGTTTGTTGGCCTTATAAGACGCACAGCGTCGCATCTGGCAACGGCTGCCg                                                                                                            <  2:1367548/148‑1 (MQ=2)
                                          tATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCg                                                                                                            >  2:823903/1‑150 (MQ=37)
                                           aTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCgg                                                                                                           >  2:3151167/1‑150 (MQ=37)
                                             gTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGAt                                                                                                         >  2:1248576/1‑150 (MQ=37)
                                                  ggTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcggc                                                                                                    >  1:709729/1‑150 (MQ=32)
                                                  ggTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTGTGTAGGCCTGATAAGACGCACAGCGACGCATCAGGCAACGGCTGCCGGATgcggc                                                                                                    >  2:4607771/1‑150 (MQ=17)
                                                      tGTAGGCCTTATAAGACCCACCGCGTCGAATCTGGCAAGGGCTGTGGGAGGCGGGGTAACCGCGTTACCCGACATCGGGTTTTGTCCCGTTTGTAGGCCTTATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTaa                                                                                                <  2:4783941/150‑1 (MQ=11)
                                                        tcggCCTGAAACGCCGCCCAGCGTCGCATCTGGCAACGGCTGTCGGTTGCGGCGTACCCGCCTTACCCGACCTAGGGTTTTGTTCGGTTTGTTGGCCTGATAAGACGCACAGCGTCGCAACAGGCAACGGCTGCCGGATGCGGCGTAAAc                                                                                              <  2:4994641/148‑1 (MQ=11)
                                                         aGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTTAACg                                                                                             >  2:1122090/1‑150 (MQ=17)
                                                             cTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCtt                                                                                         >  2:958033/1‑150 (MQ=35)
                                                                     aCGCCCCGCGTCGAATCAGGCAACGGCGGTCGTTTGCGGCGTAAACGCCTTATCCGCCCTACGGTTTTGTTCCGTTTGTCGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGAcc                                                                                 <  2:2961747/150‑1 (MQ=11)
                                                                       gCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCAGTTTGTAGGCCTGATAAGACGCACAGCGACGCATCAGGCGACGGCTGACGGAAGCGGCGTAAAAGCCTTATCAGAcaaa                                                                               >  2:1293698/1‑147 (MQ=11)
                                                                        cactGCGTCGATTCTGCCAACGGCTGTCGGATGCGGCTTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTTGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTAc                                                                              <  2:4551634/150‑1 (MQ=11)
                                                                         actgCGTCGCATAAGGAAACGGCTGTCGGTGGGGGCGTAAACCCGTTAGCCGACCTAGGGTTATGTTCCTTTTGTTGGCCTTATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACg                                                                             <  2:656630/147‑1 (MQ=11)
                                                                         acaGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTGATCCGACCTACg                                                                             >  2:4899732/1‑150 (MQ=37)
                                                                           aGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGGTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGt                                                                           >  1:1917710/1‑150 (MQ=37)
                                                                            gCGTCGCTACTGGCAAGGGCTGTGGGATGCGGCGTAAACGCCTTATCCGACATACGTTTTTTTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGtt                                                                          <  2:368370/150‑1 (MQ=11)
                                                                              gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCGTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAt                                                                        >  1:4734613/1‑150 (MQ=37)
                                                                                   aTCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTcc                                                                   >  1:5234215/1‑150 (MQ=38)
                                                                                   aTCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTcc                                                                   >  1:1060180/1‑150 (MQ=38)
                                                                                    tCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCg                                                                  >  1:1631950/1‑150 (MQ=38)
                                                                                     cAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGt                                                                 >  1:1976575/1‑150 (MQ=39)
                                                                                     cAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGt                                                                 >  1:1694470/1‑150 (MQ=39)
                                                                                      aGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGtt                                                                >  2:4549779/1‑150 (MQ=40)
                                                                                          aaCGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCAGACCTACGGTTATGTTCCGTATGTa                                                            >  2:440546/1‑150 (MQ=255)
                                                                                           ccGGCTGTCGGTTGCGGCTAAAACGCCTTAGCCGACCTCGGTTTATGTCCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAg                                                           <  2:3446083/149‑1 (MQ=14)
                                                                                             gggtGTCGGTGGCGGCTTAACCGCTTTACCCGACCTAGGTTTTTGTTGGTTTTGTCGGCTTCATACGACGCACCGCGTCGCCCCCGCCCACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGc                                                         <  1:2013042/147‑1 (MQ=9)
                                                                                             ggcgggcGGAGGCGGCGTAAACGCCTTATCCGACCTAGGTTTTTGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGc                                                         <  2:332467/144‑1 (MQ=14)
                                                                                               cTGTCGTATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCGGTTTGTTGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCt                                                       <  2:5383770/150‑1 (MQ=255)
                                                                                                 gTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGa                                                     >  1:2909946/1‑150 (MQ=255)
                                                                                                   gggATGCGGCGTAAACGACTTAGCCGACCTACGTTTATGTGCCTTTTTTAGGCCTGATAAGCCGCACAGCGTCGAAACAGGAAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATa                                                   <  2:1596078/149‑1 (MQ=14)
                                                                                                   cGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATa                                                   >  1:2209417/1‑150 (MQ=255)
                                                                                                   cGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATa                                                   >  1:4803571/1‑150 (MQ=255)
                                                                                                    ggATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCGTCTCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGAtca                                                  >  2:2140575/1‑148 (MQ=255)
                                                                                                    ggAGGCGGCGAAAACGCGTAATCCGACCTAGGTTTATGTTCCGTTTTTCGGCCTGAAACGACCCCCAGCGGCGCATCCTGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATaa                                                  <  2:3858000/150‑1 (MQ=14)
                                                                                                             gTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAg                                         >  2:5321839/1‑150 (MQ=255)
                                                                                                             gTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAg                                         >  1:2528857/1‑150 (MQ=255)
                                                                                                             gTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTGTGTAGGCCTGATAAGACGCACAg                                         >  2:4826702/1‑150 (MQ=255)
                                                                                                              tAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGc                                        >  1:4384705/1‑150 (MQ=255)
                                                                                                                  cgcgTTACCCGACATCGGGTTTTTTTCCGTTTTTGGGCCTGTTAACACGCACCGGTACCCATCTGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCg                                    <  1:4003121/150‑1 (MQ=11)
                                                                                                                  cGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCg                                    >  1:3556050/1‑150 (MQ=255)
                                                                                                                    ctttATCCGACATACGTTTTGGTTCCTTTTGTAGGCCTTATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCa                                  <  2:4823878/148‑1 (MQ=14)
                                                                                                                    ccTTACCCGACCTACGGTTATGTTCCGTTTGTCGGCCTGATAAGACGCACCGCTGCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCa                                  <  2:2061579/150‑1 (MQ=14)
                                                                                                                     cTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCAt                                 >  1:598019/1‑150 (MQ=255)
                                                                                                                       tATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCa                               >  1:3852318/1‑150 (MQ=40)
                                                                                                                        aTCCGCCCTAGGTTTATGTTCGGTTTGTAGGCCTGATATGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg                              <  2:4311273/150‑1 (MQ=11)
                                                                                                                        aTCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg                              >  2:1748974/1‑150 (MQ=40)
                                                                                                                          ccGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGc                            >  2:2466331/1‑150 (MQ=39)
                                                                                                                           cGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCa                           >  2:2907931/1‑150 (MQ=38)
                                                                                                                           cGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCa                           >  2:1809313/1‑150 (MQ=38)
                                                                                                                             aCCTACGTTTATGTACCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGAAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAAc                         <  2:1766622/150‑1 (MQ=255)
                                                                                                                              ccTAGGTTTATGTTCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACg                        <  2:779084/150‑1 (MQ=255)
                                                                                                                               cTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACgg                       >  1:2474233/1‑150 (MQ=35)
                                                                                                                                   ggTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCTTCAGGCAACGGCTga                   >  1:1097335/1‑149 (MQ=255)
                                                                                                                                   ggTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGc                   >  1:1832766/1‑150 (MQ=35)
                                                                                                                                      tATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCgg                >  2:3821138/1‑150 (MQ=35)
                                                                                                                                      tATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCgg                >  1:4150444/1‑150 (MQ=35)
                                                                                                                                         gTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCAAAGCGTCGCATCAGGCAACGGCTGCCGGCTg             >  2:1447953/1‑150 (MQ=25)
                                                                                                                                           tCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcg           >  2:3034376/1‑150 (MQ=34)
                                                                                                                                            ccGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcgg          >  2:2703133/1‑150 (MQ=33)
                                                                                                                                            ccGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcgg          >  2:1706933/1‑150 (MQ=33)
                                                                                                                                              gTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcggcg        >  1:3935592/1‑150 (MQ=32)
                                                                                                                                                   tAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAAc   >  1:4869963/1‑150 (MQ=34)
                                                                                                                                                    aGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACg  >  1:1259964/1‑150 (MQ=34)
                                                                                                                                                    |                                                                                                                                                     
CGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACG  >  CP009273/3385430‑3385727

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