Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,060 C→T 12.9% intergenic (+266/+376) 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,0600CT12.9% 145.5 / 16.5 78intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (35/31);  new base T (3/7);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

CAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAA  >  NC_000913/4295823‑4296285
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                   
cAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAAc                                                                                                                                                                                                                                     >  1:130327/1‑238 (MQ=37)
  aTGTGGGAATGCCTGATGCGACGCTTGCCGGGTCTTATCAGGCCTACGGCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCTCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGcc                                                                                                                                                                                                                                  <  1:52738/239‑1 (MQ=19)
               tGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCACTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGCCGTTTACGCCGCAT‑CCGCCATCAATGCCTGATCCGACGCt                                                                                                                                                                                                                      >  1:377180/1‑238 (MQ=9)
                gATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAATTCCTGCTGCGACGCtt                                                                                                                                                                                                                     >  1:113134/1‑238 (MQ=14)
                  ggCGACTGTTGTCGCGTTTTACCGGGCCTGAGCCAGACAGCGCAATTGCCTTATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAATGCCTGATGCGACGCTTGcc                                                                                                                                                                                                                  <  1:305514/238‑1 (MQ=11)
                     gcggtttGCCGCGTTTTATCGGGCCGACGCCAGGCAGCGCATTAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTAGGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCcgcg                                                                                                                                                                                                               <  2:611938/234‑1 (MQ=11)
                       cGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGAAATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                              <  1:376522/238‑1 (MQ=18)
                       acagtGCCGCGTCTTATCAGGCCTGCGCCAGACACGGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                              <  1:459422/234‑1 (MQ=255)
                                tttttttttgggCCTGCGCCAGACGTCGTAATAGCATGATTTAGGGTGTTTTTGTAGGTCGGATAAGGCTTTTATGCCGCATCCGACATCAACGCCTGATGCGACTCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                    <  2:348562/229‑1 (MQ=11)
                                 gTCTTATCAGGCCTACGGCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACAGCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                    <  1:817675/238‑1 (MQ=18)
                                 gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                    <  1:561622/238‑1 (MQ=18)
                                  tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGCCCCTTGCGCCGGCTTTTCAGGCCt                                                                                                                                                                                                  >  1:1069500/1‑239 (MQ=255)
                                      aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTACCTTGATTTTGTAGTTCGGATAAGGCGTTTA‑GGCGCATCCCGCCATCATGGCCTGTTGCGCCGCTGGCCCCGGCTTATccgcccaactt                                                                                                                                                                                              >  2:509343/1‑229 (MQ=11)
                                      aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTTTAGGTCGGATAAGGCGGTTACGCCGACT‑CCGCCACCAATGCCTGATGGCACCCTTGCCGCGGCTTACAAGGccctgt                                                                                                                                                                                               >  2:459422/1‑234 (MQ=11)
                                             cgaggCCAGACAGCGGAATAGCCTTGTTTGGCGTGATTTTGCAGGTCGGATAAGTCGTTTATGCCGCCTCCGACATCAACGGCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACa                                                                                                                                                                                        <  1:44047/234‑1 (MQ=11)
                                                      aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                                                                                                              <  1:914969/239‑1 (MQ=33)
                                                          cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCCAGGCTTACCGCGTCTTATCAGGCCTACGCCGGACACAGAAAAAGCCTgg                                                                                                                                                                          >  2:1069373/1‑238 (MQ=11)
                                                               tggCCTGATTTAGCGTGATTTTGTAGGTCGGATGAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAg                                                                                                                                                                     <  1:998268/237‑1 (MQ=18)
                                                                  ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACGTCAACGCCTGCTGCGCCGCTTACCGCGTCTTATCAGGCCTACCCCACACCGCGCAAAGGCTTGTTTTacctg                                                                                                                                                                  >  2:569245/1‑235 (MQ=11)
                                                                               gtttttttTGGGCGGGTAAGGGTTTTTTTGCGCGTCCGACATCAACCGCTGATGGGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACGGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGt                                                                                                                                                     <  2:923251/236‑1 (MQ=11)
                                                                                   ttgggggTCGGATAAGGCGTTTGCTCCCGTTCCGCCATCAACGCCGGATGGGACGTTTAACGCGGTTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGAt                                                                                                                                               <  2:130327/234‑1 (MQ=11)
                                                                                      tAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                    <  1:343050/185‑1 (MQ=35)
                                                                                      tAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                    >  2:343050/1‑185 (MQ=35)
                                                                                           cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTACCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGt                                                                                                                                                     >  2:226246/1‑226 (MQ=21)
                                                                                           cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTACCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGt                                                                                                                                                     <  1:226246/226‑1 (MQ=21)
                                                                                                       ttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCCGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACCCCAACCAGGCCACTAGCCTGTTTTAG‑GTGGTGTTTTTAGGTC‑GGGAAAGGGGGTTAGGcccgaacc                                                                                                                           >  2:561622/1‑233 (MQ=11)
                                                                                                        tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACACCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGCC‑GGATAAGCCGTTTATGCCGCAccc                                                                                                                           >  1:453842/1‑238 (MQ=255)
                                                                                                              gCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCTATAGCCTGGTTTAGCTGCAT‑TTTCGCGGTC‑GGAACTGCCGTTTATGCCCCCTCCCACATCa                                                                                                                    >  2:471778/1‑239 (MQ=11)
                                                                                                                aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACGTCAAACAGCGCAATAGCCTTCTTTAGCGTCAT‑TTCTTGGTTC‑GGAAAAGTCGCTTACACCGCATCCGACATaaaa                                                                                                                  >  2:14422/1‑238 (MQ=11)
                                                                                                                  ccGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                    <  1:788579/157‑1 (MQ=35)
                                                                                                                           atgCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCg                                                                                                        >  1:1071756/3‑238 (MQ=14)
                                                                                                                            cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                                                                                                          >  2:788579/1‑157 (MQ=35)
                                                                                                                                 gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                                                                                                          >  2:128194/1‑136 (MQ=37)
                                                                                                                                 gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGACTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                                                                                                          <  1:128194/136‑1 (MQ=25)
                                                                                                                                  gtgggACGGTTAACGCGTCTTATCAGGCCTACGCCGGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTaa                                                                                                <  1:385412/235‑1 (MQ=255)
                                                                                                                                  aTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGAGGCTTa                                                                                                 >  1:278993/1‑238 (MQ=37)
                                                                                                                                   tGCGACGATTAACGGGTCTTATCAGGCCTAAGCCAGACAGCGGAATTGCCTGAATTAGCGTAATTTTGTAGGCCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTaa                                                                                                <  1:75618/238‑1 (MQ=11)
                                                                                                                                             tcccctttttttCAGGCCTACCCCAGACGCGGCAATGGCCTGATTTAGCGTGTTTTGGTAGGTGGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                     <  2:927175/228‑1 (MQ=11)
                                                                                                                                                 cttCTTTTCAGGCCTACGCCAGGCGGCGAAATGGCCTTATTTGGCTTGATTTTGTAGGTCGGATAAGGCTTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                 <  2:458111/237‑1 (MQ=11)
                                                                                                                                                    cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGCCTTATCAGGCCt                                                                               >  1:125978/1‑238 (MQ=18)
                                                                                                                                                            gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                    <  1:165501/115‑1 (MQ=39)
                                                                                                                                                            gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGCGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATCTGACGCTTAACGCGTCTTATCAGGCTCACGCCAgcc                                                                      >  1:96697/1‑237 (MQ=255)
                                                                                                                                                             ccTACGCCAAGCAGGGAAATTGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCCTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTACCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGa                                                                       <  1:1004403/237‑1 (MQ=255)
                                                                                                                                                                   ccAGACTGCTGAATTGCCGGATTTAGCGGGATTTTGTAGGGCGGAGAAGGCGGTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCaa                                                               <  1:536956/239‑1 (MQ=11)
                                                                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCGACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                    >  2:916664/1‑100 (MQ=255)
                                                                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                    <  1:916664/100‑1 (MQ=255)
                                                                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                    >  1:246002/1‑100 (MQ=40)
                                                                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                    <  1:143737/100‑1 (MQ=40)
                                                                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                                                                                                          >  1:371175/1‑110 (MQ=40)
                                                                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCaa                                                                                                                                                                                  <  1:610066/118‑1 (MQ=39)
                                                                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATaa                                                                                                                                             >  2:729622/1‑153 (MQ=35)
                                                                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATg                                                                                                                                   <  1:658692/163‑1 (MQ=35)
                                                                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATg                                                                                                                                   >  2:658692/1‑163 (MQ=35)
                                                                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAAAGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTTGCGTGTT‑TTTGTAGACC‑GGCTAAGGCCTTTACGCCGCCTCCGCCATCCACGCCTTGTCCGACGCTTAACCCGCTTTTTCAGGCCCACGCCACCCAGAGCCACAACCTGa                                                       >  2:97669/1‑239 (MQ=11)
                                                                                                                                                                             ctattgCCTGGTTTGGGGTGATTTTGGAGGTCGGATAAGGCGTTTACGGCGCAT‑CAGGCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAtt                                                     <  2:1060228/234‑1 (MQ=11)
                                                                                                                                                                              aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGACGCGTCTTATCAGGc                                                                                                                                                                                                    >  1:546697/1‑97 (MQ=18)
                                                                                                                                                                                tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGt                                                                                           >  2:99307/1‑198 (MQ=34)
                                                                                                                                                                                tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGt                                                                                           <  1:99307/198‑1 (MQ=34)
                                                                                                                                                                                     tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATACGGCGTTTATGCCGCATCCGACCTCAACGCCTGATGCGACGCTTAACTCGTCTTATCAGCCCTACCCCAGACACCGAAAAAGCCTGCTTTACCGGGa                                             >  2:867043/1‑239 (MQ=11)
                                                                                                                                                                                             cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGAGGCTTAACGCGTCTTTTCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTCACCGAGCTTTTGAAg                                     >  2:280193/1‑239 (MQ=11)
                                                                                                                                                                                                gtatgtGTAGGGTGGAGAAGGCGTTTCCGCCGCAT‑CCGAGATCAACGCCTGATGCGACGGTTGCCGGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTCGGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGt                                   <  1:1075558/234‑1 (MQ=11)
                                                                                                                                                                                                           tCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCaa                                                               <  1:627602/199‑1 (MQ=34)
                                                                                                                                                                                                            cGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAATGCCTGATGCGCCGCTTAACGCGTCTTATCAGGCCTACGCCAGCCGGCGCATAAGCCTGATTTCGGGTGGTTTTTTAGGCCGTAAAAGGCgat                      >  2:709676/1‑237 (MQ=11)
                                                                                                                                                                                                            cGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGCCGCTTAACGCGTCTTATCAGGCCTAAGCCAGACAGCGCATTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAACGCGtt                      >  1:304860/1‑239 (MQ=16)
                                                                                                                                                                                                            cGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGa                                             <  1:254901/216‑1 (MQ=34)
                                                                                                                                                                                                            cGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGa                                             >  2:254901/1‑216 (MQ=34)
                                                                                                                                                                                                            cGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGGCTTATCAGCCCTACGCCCAACAGCGGAATACCCTGATTTAGCGCGATTTTGAAGGTTGGATAAGGAGtt                      >  2:204183/1‑239 (MQ=11)
                                                                                                                                                                                                                     cGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAtt                                                     >  2:627602/1‑199 (MQ=34)
                                                                                                                                                                                                                         tACGCCGCAT‑CCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                              <  2:490927/44‑1 (MQ=12)
                                                                                                                                                                                                                         tACGCCGCAT‑CCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                    <  2:946679/54‑1 (MQ=9)
                                                                                                                                                                                                                         tACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                              >  1:490927/1‑44 (MQ=255)
                                                                                                                                                                                                                         tACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                    >  1:946679/1‑54 (MQ=255)
                                                                                                                                                                                                                         tACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTCCGCCCTCCAGCCCAAAAGCCTTATTCTGCGTGATTTTGTATGTCGGATAAGGCGTTTAGGCACCATccc         >  2:796671/1‑238 (MQ=11)
                                                                                                                                                                                                                         tACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCACCGCAATAGCCTACGTTAGCGGGATTTTGTAGGTCCGATAAGGGGGTTCAGCCGCATCCg         >  2:1087986/1‑239 (MQ=11)
                                                                                                                                                                                                                         tACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACACCGCAATAGCCTTATTTAACGTGATTTTGTAGGTCGGATAAGGCGGTTACGCCgaactcg         >  2:369136/1‑233 (MQ=255)
                                                                                                                                                                                                                         tACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTg                                                                                                                                                          <  1:668385/95‑1 (MQ=18)
                                                                                                                                                                                                                               cgcacgcggcATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATc    <  1:792649/230‑1 (MQ=9)
                                                                                                                                                                                                                                cAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCTCGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATccccaataaa  >  2:931910/1‑231 (MQ=255)
                                                                                                                                                                                                                                            aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGa                                                                                                       <  2:23182/127‑1 (MQ=32)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                   
CAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAT‑TTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAA  >  NC_000913/4295823‑4296285

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