Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,385,702 A→C 6.5% intergenic (‑315/+115) 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,7020AC 6.5% 594.8 / 7.5 232intergenic (‑315/+115)tldD/yhdPputative peptidase/DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base A (143/66);  new base C (2/19);  total (146/86)
Fisher's exact test for biased strand distribution p-value = 1.54e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 241 total aligned reads displayed.

CGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAAC  >  CP009273/3385553‑3385850
                                                                                                                                                       |                                                                                                                                                    
 gACCTCGGGTTCTGTTCCGTTTTTAGGCCTTATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCaa                                                                                                                                                     <  2:3475918/150‑1 (MQ=11)
 gACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTAGCCGACCTACGGTTATGTTCC‑GTGTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCaa                                                                                                                                                     >  2:3241420/1‑150 (MQ=255)
  acctcgggTTATGTTCCTTTTGTAGGCATTTTAAGACCCACACCGTCGCATCAGCCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAAc                                                                                                                                                    <  2:1122901/144‑1 (MQ=11)
    cTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACgg                                                                                                                                                  >  2:352806/1‑150 (MQ=35)
    cTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACgg                                                                                                                                                  >  2:425671/1‑150 (MQ=35)
        ggTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTATGCC‑TGATAAGACGCACAGAGACGCATCAGGCAACGGCTGc                                                                                                                                              >  2:3046344/1‑150 (MQ=255)
        ggTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTga                                                                                                                                              >  2:2438507/1‑149 (MQ=25)
          ttATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCATATCCGACCTACGGTTATGTTCA‑GGATGTAGGCC‑TGATAAGACGCACAGAGTCGCATCAGGCAACGGCTGCCg                                                                                                                                            >  2:995629/1‑150 (MQ=11)
           tATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCgg                                                                                                                                           >  2:3946423/1‑150 (MQ=35)
            aTGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGa                                                                                                                                          >  1:3955123/1‑150 (MQ=35)
            aTGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGa                                                                                                                                          >  2:1584175/1‑150 (MQ=35)
               ttCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgc                                                                                                                                       >  1:2786242/1‑150 (MQ=35)
                  cGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGGTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcggc                                                                                                                                    >  2:2292572/1‑150 (MQ=17)
                   ttttttgggCCTTATACGCCGCACCGGGTCGCATCTGCCAACGGCTGCGGTAGGCGGCGTAAACGCCTTACCCGACCTCCGGTTATGTTCG‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGAATCTGGCAACGGCTGCCGGATgcggcg                                                                                                                                   <  2:2066265/143‑1 (MQ=11)
                     ttGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTCTCCGACCTACGGTTAGGTTCC‑GGTGGGAGGCC‑TGATACGACGCACAGCGTCGCGTCAGGCCACGGCTGCCGGCGGCGGCGTg                                                                                                                                 >  2:2475578/1‑150 (MQ=11)
                     ttGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCCACGGCTGCCGGATGCGGCGTg                                                                                                                                 >  2:301195/1‑150 (MQ=17)
                      tGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGa                                                                                                                                >  2:757431/1‑150 (MQ=33)
                      tGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGa                                                                                                                                >  1:1736666/1‑150 (MQ=33)
                      tGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGa                                                                                                                                >  2:2150572/1‑150 (MQ=33)
                        tAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAAc                                                                                                                              >  1:473321/1‑150 (MQ=34)
                        tAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAAc                                                                                                                              >  2:915910/1‑150 (MQ=34)
                               gATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTAt                                                                                                                       >  1:3101067/1‑150 (MQ=37)
                                aTAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATc                                                                                                                      >  2:784839/1‑150 (MQ=37)
                                 tAACACGCCCTGCGACGAATCTGGCAACGGCTGCGGTATGCGGTTTAAAGGCCTTATCCGACCTCGGTTTATTTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCTTCGCATCAGGAAACGGCTGCCGGATGCGGCGTGAACGCCCTATcc                                                                                                                     <  2:1356108/150‑1 (MQ=11)
                                  aaGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGACTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCg                                                                                                                    >  2:2273360/1‑150 (MQ=37)
                                     aCGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCACCGGCTGCCGGATGCGGCGTGAACGCCCTATCCGAcc                                                                                                                 >  2:22990/1‑150 (MQ=37)
                                      cccactGCGTCGCTCCTGGCACGGGCTGCGGGATGGGCCGTAAACGCCTTACCCGACCTACGTTTTGGTGCG‑GTTTGTAGGCC‑TGATAAGACGCACAGGGTCGCTACTGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCt                                                                                                                <  2:3341986/148‑1 (MQ=2)
                                      cGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCt                                                                                                                >  1:2988161/1‑150 (MQ=37)
                                         acaGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTCAACGCCTTATCCGACCTACGGTTATGTTCC‑GGTTGTAGGCC‑TGATAAGACGCACCGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCGATCCGACCTACg                                                                                                             >  2:4133970/1‑150 (MQ=255)
                                                 gCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGtt                                                                                                     >  1:2133138/1‑150 (MQ=37)
                                                   aTCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTCCGgttctgttct                                                                                                   >  2:4269771/1‑150 (MQ=255)
                                                       ggCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCt                                                                                               >  1:4096055/1‑150 (MQ=40)
                                                           aCGGCTGCCGGATGCGGCTTAACCGCCTTATCCGACCTACGTTTTTTTTCG‑GTTTGTAGGCC‑TGATAAGACGCACCGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGtt                                                                                           <  2:965416/150‑1 (MQ=14)
                                                            cGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTg                                                                                          >  2:2540479/1‑150 (MQ=255)
                                                             gggggCCGGAGGCGGCTTAAACGCCTGACCCGACCTACGGTTATGTTCC‑GATTTTTGGCC‑TGATAAGACGCACAGCTGCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTgg                                                                                         <  2:3604430/146‑1 (MQ=14)
                                                              ggtGCGGGAGGCGGCTTAAACGCTTTAGCCGCCCTCGGTTTTTGTTCC‑GTTTGTCGGCC‑TGATCAGACGCACCGCGACGCTTCCGGCCACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGc                                                                                        <  1:4347960/148‑1 (MQ=14)
                                                              gCTGCCGGATGCGGCTTAAACGCCTTATCCGACCTCCGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGATGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGc                                                                                        <  1:1154868/150‑1 (MQ=255)
                                                                tGCCGGATGCGGCGTAAACGCCTTATCCGACCCACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCt                                                                                      >  2:3098988/1‑150 (MQ=255)
                                                                tGCCGGAGGCGGCGTAAACGCCTTTACCGACCCACGGTTTTTGTCC‑TTTGTTAGCCC‑GCAAAAGACGCACCGCGCCGCCACAGTCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCt                                                                                      <  2:977256/150‑1 (MQ=17)
                                                                 gggggATGCGGCTAAACCGCCTTACCCGACCTAGGTTTTTTTTCG‑TTTTGTTGGCA‑TTATAAGACCCACAGCGTCGAACCAGCCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCtt                                                                                     <  2:50382/147‑1 (MQ=14)
                                                                 gcgggAGGCGGCTTAACCGATTTATCCGCCCTCGGTTTTTGTTCG‑TTTTGTAGGCG‑TGATAAGACCCACCGCGTCCAATCTGCCCACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCtt                                                                                     <  1:4224338/147‑1 (MQ=9)
                                                                 gCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCtt                                                                                     >  2:1024779/1‑150 (MQ=255)
                                                                  ccGGAGGCGGCGTAAACGCCTGACCCGACCGCCGGTTATTTTCC‑GTTTGTCGGCC‑TGTTAAGACGCCCAGCGTCGCATCAGGCTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCttt                                                                                    <  2:2844098/150‑1 (MQ=17)
                                                                   ggttgggggctTAAACCCTTTCCCCGACCTAGGTTTATTTTCT‑TTTTGTCGGCC‑TTATACGACCCACAGCGCCGCATCCGCCCACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTg                                                                                   <  1:1491417/143‑1 (MQ=12)
                                                                    ggATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTgg                                                                                  >  1:931942/1‑150 (MQ=255)
                                                                       tGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCaa                                                                               >  1:517727/1‑150 (MQ=255)
                                                                       tGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCaa                                                                               >  2:1946632/1‑150 (MQ=255)
                                                                        gcggctAAACCGCATTACCCGACCTA‑GGTTTTTTTTC‑GTTTTTTTGCCATTATATGACGCACAGCGTCGCATCAGCCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAAc                                                                              <  1:4184920/150‑1 (MQ=17)
                                                                         cggcgTAACCGCCTTAGCCGCCCTCCGGTTTTGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACCGCGTCGCATCAGCCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACg                                                                             <  2:2655530/150‑1 (MQ=17)
                                                                         cggcgTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACg                                                                             >  2:4157859/1‑150 (MQ=255)
                                                                            ggTACACGCCTACTCCGCCCTACGTTAATTTCCC‑GTTTGTAGGCC‑TAACAAGAGTCACAGCTTGTCTTATGTCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGAtt                                                                          <  2:4324925/149‑1 (MQ=12)
                                                                                aaCGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATcc                                                                      >  1:659876/1‑150 (MQ=255)
                                                                                aaCGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATcc                                                                      >  2:3192915/1‑150 (MQ=255)
                                                                                aaCGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATcc                                                                      >  1:1539116/1‑150 (MQ=255)
                                                                                aaCGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATcc                                                                      >  1:865734/1‑150 (MQ=255)
                                                                                aaCGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATcc                                                                      >  1:772193/1‑150 (MQ=255)
                                                                                 aCGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCt                                                                     >  1:2856833/1‑150 (MQ=255)
                                                                                   gCCTTATCCGACTTCGGGTTTTGTACG‑TTTTTTTGGCC‑TGATATGACGCACCGCGGCGCATCAGCCCACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTAt                                                                   <  1:2355708/150‑1 (MQ=18)
                                                                                     cTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATga                                                                 >  1:1998974/1‑150 (MQ=255)
                                                                                      ttATCCGACCTTGGGTTTTGTCCG‑TTTTGTAGGCC‑TTAAATGACGCACCGCGGCGAATCCGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATgag                                                                <  1:3201790/150‑1 (MQ=18)
                                                                                         tCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGAtt                                                             >  1:1144706/1‑150 (MQ=255)
                                                                                         tCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGAtt                                                             >  1:2588118/1‑150 (MQ=255)
                                                                                            gACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTggg                                                          >  1:1905753/1‑150 (MQ=255)
                                                                                                   ggTTATGTTCCGGTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCACCCAATCTCCTCAGAGATTCCGGCAAtt                                                    >  2:2599585/1‑150 (MQ=17)
                                                                                                    gTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACCCACCGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTAc                                                  <  2:2166175/150‑1 (MQ=255)
                                                                                                     ttttGTCCC‑TTTTTTAGCCC‑TCATACTACCCACAGCGTCCCAACAGCCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTAcg                                                 <  2:5975/147‑1 (MQ=25)
                                                                                                         gTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGcc                                             >  1:3746562/1‑150 (MQ=255)
                                                                                                          ttCC‑GTTTGTAGGCC‑TTATAATCCCCACAGCGTCCCACCCGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGccc                                            <  2:830397/150‑1 (MQ=25)
                                                                                                              ggTTTGTAGGCC‑TTATAAGACGCCCAGCTTCGCATCTGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCg                                         <  2:170738/149‑1 (MQ=255)
                                                                                                                 ttttAGGCA‑TTAAAACACCCACAGCGGCCAAACAGCCCACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCa                                      <  1:2438507/147‑1 (MQ=25)
                                                                                                                  tGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTCTGAGATTGGGGCAATTACGCGCCCTCGTCaa                                     >  1:3943007/1‑150 (MQ=255)
                                                                                                                  tGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCaa                                     >  1:3941561/1‑150 (MQ=255)
                                                                                                                     aGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATc                                  >  1:230264/1‑150 (MQ=255)
                                                                                                                     aGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATc                                  >  1:326984/1‑150 (MQ=255)
                                                                                                                      ggCC‑TTATACGACGCACTGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCa                                 <  2:452179/150‑1 (MQ=255)
                                                                                                                      ggCC‑TTATAAGACGCACAGCGTCGAATCTGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCa                                 <  2:1732573/150‑1 (MQ=255)
                                                                                                                        cc‑TGATTCGACGCACAGGGTCGAATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCAtt                               <  1:1046233/150‑1 (MQ=255)
                                                                                                                         c‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTg                              >  2:2023092/1‑150 (MQ=255)
                                                                                                                          gttaTACGACGCACCGCGTCGCATCAGGCCACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTg                              <  1:3507509/147‑1 (MQ=255)
                                                                                                                           ttaTAAGCCGCACCGCGTCGCATAAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTgc                             <  2:1092179/148‑1 (MQ=255)
                                                                                                                           tGATAAGACGAACCGCGCCGACGCCGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTgc                             <  1:2762431/150‑1 (MQ=37)
                                                                                                                            gATAAGACGCACCGCGTCGCATCAGGCCACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTgcg                            <  2:3446808/150‑1 (MQ=255)
                                                                                                                             aTAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCCAATCATTgcgc                           >  1:2916472/1‑150 (MQ=255)
                                                                                                                              tAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCt                          >  1:71577/1‑150 (MQ=255)
                                                                                                                               aaGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCtt                         >  2:3633425/1‑150 (MQ=255)
                                                                                                                                   cGCACCGCGTCGAATCAGGCAACGGCTGTCGGATGCGGCGTCAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTAAGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTc                     <  2:2513015/150‑1 (MQ=37)
                                                                                                                                   cGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTg                                                                                                  <  2:815314/73‑1 (MQ=255)
                                                                                                                                   cGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTc                     >  2:724920/1‑150 (MQ=255)
                                                                                                                                     caccGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTCTGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCtt                   <  2:4254483/150‑1 (MQ=255)
                                                                                                                                     cacaGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCtt                   >  2:2829554/1‑150 (MQ=255)
                                                                                                                                     cacaGCGTCGCAACAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCtt                   <  2:960296/150‑1 (MQ=255)
                                                                                                                                        aGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTAc                >  2:153539/1‑150 (MQ=255)
                                                                                                                                         ggggcGCATCCGGCCACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTAcg               <  1:2996351/146‑1 (MQ=255)
                                                                                                                                           gTCGCATCAGGCCACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTAcgcg             <  2:3955123/150‑1 (MQ=255)
                                                                                                                                           gTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTAcgcg             >  1:2822227/1‑150 (MQ=255)
                                                                                                                                            tCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCgg            >  2:3588738/1‑150 (MQ=255)
                                                                                                                                               cATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTg         <  1:2150572/150‑1 (MQ=255)
                                                                                                                                               cATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTg         <  2:2887975/150‑1 (MQ=255)
                                                                                                                                                aTCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTgg        >  1:2305610/1‑150 (MQ=255)
                                                                                                                                                   aGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGgcgc     >  2:3453051/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
CGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC‑GTTTGTAGGCC‑TGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAAC  >  CP009273/3385553‑3385850

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