Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,190 A→G 19.5% intergenic (+396/+246) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,1900AG19.5% 65.4 / 12.9 41intergenic (+396/+246)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base A (16/15);  new base G (3/5);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.06e-01

ATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCC  >  NC_000913/4295903‑4296471
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                            
aTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCCCAGACAGCGCAATAGCCTGATTAAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCCCATCCGGAATCAACCCCCGATGGCACGCTTa                                                                                                                                                                                                                                                                                              >  2:29793/1‑288 (MQ=255)
    tGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATAAACGCCTGATGGGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTGGCGTGATTTTGTAGGTCGGAGAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgc‑g                                                                                                                                                                                                                                                                                        <  2:405373/289‑1 (MQ=19)
      tgagTCGGATAAGTCGGTGATGCCCCATACGAAATCAACGCATGATGCGACGCTTACCGCGTCTTATCAGGCCTACGGCAGACGGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑Gt                                                                                                                                                                                                                                                                                       <  2:254700/285‑1 (MQ=11)
           cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTAt                                                                                                                                                                                                                                                                                  >  2:104805/1‑288 (MQ=37)
           cGGATAAGGCGGTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTATGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTA‑TGCCGCATCCGACATCCACGCTTGATGGGACGATTACCGC‑GTCTTAt                                                                                                                                                                                                                                                                                  >  1:706332/1‑288 (MQ=11)
           cGGAAAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTAt                                                                                                                                                                                                                                                                                  <  1:123811/288‑1 (MQ=37)
               tAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACTCTTATCGGGTCTTATCGGGCCTAAGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGc                                                                                                                                                                                                                                                                             <  2:95668/289‑1 (MQ=255)
                 cggggTTTATGCCGAGACCGACATCAACTCTTGATGTGACGATTACGGCGACTGATCAGGGCTACGCGAGACAGGGCAATGGTCTGATTTAGCGGGATTTTCTAGGTCGGATAAGGGGTTTACGCCGCATCCGAGATCAACGGCTGATGCGACGCTTGCGGCGTCTTGTCATGCCTACGCCAAACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGc                                                                                                                                                                                                                                                                             <  2:352705/284‑1 (MQ=11)
                   gtgTTTTTGCCCAACCCGATATCAATGGCTGATCCAACGCATAACGCGTCTTATCAGGACGACGCCAGACAGCGCAATAGCCTAATTTAGCGTTATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGc                                                                                                                                                                                                                                                                             <  2:661862/283‑1 (MQ=11)
                         aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGa                                                                                                                                                                                                                                                                   <  1:104805/289‑1 (MQ=37)
                             cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCCACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGc                                                                                                                                                                                                                                                                             >  2:285754/1‑275 (MQ=25)
                             cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGc                                                                                                                                                                                                                                                                             <  1:285754/275‑1 (MQ=35)
                               cATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATATCCTGATTTAGCGTGATTTTGGAGGTCGGATAAGGCGTTTA‑GGCCGCATCTGACATCAACGCCTGACGCCACGCTTTACGC‑GTCTTATCCGGCCTACGCCAGACAgccc                                                                                                                                                                                                                                                             >  2:181018/1‑287 (MQ=255)
                               cATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACCTCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAgcgc                                                                                                                                                                                                                                                             >  1:99009/1‑289 (MQ=25)
                                 tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCaa                                                                                                                                                                                                                                                           <  1:370589/289‑1 (MQ=35)
                                         cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATCAGGCGTTTA‑TGCCGCATCCGACCTCAACGCCTGATTCGACGCTTAACGC‑GTCTTACCAGGCCTACCCCAGACAGCGCAATATCCTg                                                                                                                                                                                                                                                    >  2:720649/1‑288 (MQ=255)
                                                       aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCTCCGCTTAACGT‑GTCTTATCCGGCCTACGCTCGACAGCGCAATAGCCTGATTTAGCGTGAgttt                                                                                                                                                                                                                                     >  2:629602/1‑289 (MQ=255)
                                                       aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAtttt                                                                                                                                                                                                                                     >  1:302861/1‑289 (MQ=35)
                                                                               tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGACTTAGCCTGATTTTGTAGGTCCGATAAAGCTTTTacg                                                                                                                                                                                                              >  2:553792/1‑286 (MQ=12)
                                                                                 cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGTCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCCAACAGCGCAATAGCCTGTATTAGCTTGCTTTTGTATGTCGGATCAGGCGTTTACGcc                                                                                                                                                                                                            >  2:440833/1‑288 (MQ=11)
                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACATCGCACTCGCCTGATTTAGCTTTATTTCGTAGGTCGGATAAGGCGTTTACGCCGCACCCGAcc                                                                                                                                                                                                  >  2:528212/1‑287 (MQ=255)
                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAc                                                                                                                                                                                                                                                                                            >  1:53445/1‑199 (MQ=34)
                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACTCCAGACCTCGCAATAGCCTGATTTAGCCTGATTTTGTAGTTCGGATAAGGCGTTTCCGCCCCATCCGaac                                                                                                                                                                                                  >  2:350646/1‑286 (MQ=255)
                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCTCGACAGCGCAATAGCCTGATTGCGCGTTATTTTGTATGTCCGATCATGCCTTTACGCCGCATTCGACAt                                                                                                                                                                                                 >  2:540342/1‑289 (MQ=11)
                                                                                               aTAGCCTGATTTAGCGGGAGTTTGTGGGTCGGTTAAGGCGTTTTCGCCGGATCCTACGTCAACGCCTGATGCGACGCTTACCGGGTCTTATCCGGCATACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑CGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaat                                                                                                                                                                                             <  2:587549/289‑2 (MQ=11)
                                                                                                                          gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                                                                    >  1:223955/1‑287 (MQ=17)
                                                                                                                              gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACa                                                                                                                                                                                                                                                                 <  1:507251/190‑1 (MQ=35)
                                                                                                                                        ttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                                                                                                                                                                                       >  2:507251/1‑190 (MQ=35)
                                                                                                                                                                cTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGa                                                                                                                                                                                  <  1:390839/235‑1 (MQ=14)
                                                                                                                                                                cTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGAATCCGACATCAACGCCTGATGCGa                                                                                                                                                                                  >  2:390839/1‑235 (MQ=14)
                                                                                                                                                                                                  cgcgagagcgcgcAATAGACTGATTTAAGGGAATTGTGTGGGTCGGATACAGCGTTTA‑TGCCCCGACCGACATGAACTCCTGAGGCGACGCTTATCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑cgccg                                                                                        <  2:296055/280‑5 (MQ=11)
                                                                                                                                                                                                            cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCaa                                                                                                                                          <  1:44176/231‑1 (MQ=33)
                                                                                                                                                                                                            cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCaa                                                                                                                                          >  2:44176/1‑231 (MQ=33)
                                                                                                                                                                                                                        aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑cg                                                                                           <  2:623955/264‑2 (MQ=14)
                                                                                                                                                                                                                                                 tAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                           <  1:609434/289‑1 (MQ=14)
                                                                                                                                                                                                                                                           tacgCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTaa                                  <  1:553792/286‑1 (MQ=2)
                                                                                                                                                                                                                                                                 gCATCCGACATCAACGCCTGATGCGACGCTTGCCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                            >  1:581922/1‑289 (MQ=12)
                                                                                                                                                                                                                                                                   aTCCGACATCAACGCCTGATGCGACGCTTGCCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACTCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                          >  2:700199/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                               ggCCTGGTGCGACGCTT‑GCGCGGTCTGATTAGGCCGAAGCCAAGTAGCTAAATAGCCTGATTTAGGGTGATTTCGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCGGATGCGACGCTTAACGCGTCTTTTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa              <  2:430554/288‑1 (MQ=11)
                                                                                                                                                                                                                                                                               cGCCTGATGCGACGCTTGCCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa              <  2:302861/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                               cGCCTGATGCGACGCTGGCCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa              <  1:528212/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           cGCTTGCCGC‑GTCTTATCAGGCCTACGCCAGACGGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc  <  1:181018/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                            
ATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCC  >  NC_000913/4295903‑4296471

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑