Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,060 C→T 16.0% intergenic (+266/+376) gltP → / ← yjcO glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,0600CT16.0% 128.8 / 24.2 71intergenic (+266/+376)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base C (36/23);  new base T (2/9);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

GAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGC  >  NC_000913/4295797‑4296319
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                   
gAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGCCGTTTACGCCGCATCCGACATCAACGCCTGGTGCGACGCTTg                                                                                                                                                                                                                                                    >  1:758509/1‑281 (MQ=255)
   gATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcc                                                                                                                                                                                                                                                  >  1:368884/1‑280 (MQ=255)
        ccagggATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGTTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGTCAGCGCAATAGCCTGATGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCTGCATCCGACACCAACGTCTGATGCGACGCTCGCcgcgcc                                                                                                                                                                                                                                             >  1:649523/4‑278 (MQ=255)
              atatCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCTGATAAGGCGTTTACGCCGCATCCGACCTCAACGCCTGATTCGACGCTTcgcgcgcggtatca                                                                                                                                                                                                                                       >  1:206049/1‑272 (MQ=255)
                taccACCCGGCAATGTGTGAATGCCTGATCGGACGCTTGCCCCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACTCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                                    <  2:366172/279‑1 (MQ=255)
                atcaccccGGCAATGTATGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                                    <  1:305597/276‑1 (MQ=255)
                 tCCACCGTGCAATGTGTGAATGCCTGATGCGACCCTTGCCGCGTCTTATCAGGCCTACGCCCGACAGCGCAATAGCCTGATTTGGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                                    <  1:63872/280‑1 (MQ=255)
                 tCCACCCGTCAATGTGTGAATGCCTGCTGCGACGCTTGCCGCGTCTTAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATAAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATTAGGc                                                                                                                                                                                                                                    <  1:809415/280‑1 (MQ=255)
                 acgaccagagaaTGTGTTAAGGCCTGATGGGACGCTTGCCGCGTCTTACCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                                    <  1:730258/270‑1 (MQ=255)
                  ccACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCCTCCCGCGTCTTATCAGGcc                                                                                                                                                                                                                                   >  1:363943/1‑280 (MQ=255)
                       cGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGc                                                                                                                                                                                                                              >  1:180892/1‑280 (MQ=39)
                                gtgAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGACTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCAGCATCCGACAGCAACGCCTGATGCGACGCTTGCCGCGGCTTATCAGGTCTACGCCAGACAgcgc                                                                                                                                                                                                                    >  1:488362/1‑281 (MQ=255)
                                       gcggATGCGAGGCTGGCCGCGTTTTATCCGGCCTACGCCAGCAGGCGCAATCGCCTATTTTGGCGGTATTTTGGAGGTCGGGTAAGGCGTTTATGCCGCATCCGAGATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCGAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                                                                                                                                              <  2:644173/277‑1 (MQ=11)
                                       ccttggggtaCGCTTGCGGCGGTTTATCAGGCCTAACCCGGAAGGCGAAATGGCCTGATTTGGCGTGATTTTGTAGGTCGGAGAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                                                                                                                                              <  2:971371/271‑1 (MQ=11)
                                         tGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTg                                                                                                                                                                                                           >  1:287873/1‑281 (MQ=33)
                                         tGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCCGATTTAGCGTGATTTTGTAGGTTGGATAAGGCGTTTACGCCGCATACGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCATACCCCAACCAGCGCAATACCCTgc                                                                                                                                                                                                          >  1:156060/1‑280 (MQ=11)
                                              cgccGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTTATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGTTTGc                                                                                                                                                                                                                                                   <  1:806648/233‑1 (MQ=12)
                                                         cgttCTTGTCAGTCCTACGCCGGAAAGCACAATATCCTGGTTTAGGGTGATTTTGTAGATCGGATAAGGCGTTGATGCCGCGTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGGGCAATAGGCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTg                                                                                                                                                                                           <  1:627418/278‑1 (MQ=11)
                                                            tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGGTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGc                                                                                                                                                                                                                                                   <  1:534548/222‑1 (MQ=17)
                                                            tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGGTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGc                                                                                                                                                                                                                                                   >  2:534548/1‑222 (MQ=17)
                                                                          cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCACGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAACCGACATCAAGGCCTGATGCGACGCTTGCCGCGTCTTACAAGGCCTACGCCAGACCGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTGGGCTAAgggg                                                                                                                                                                           >  1:763247/1‑278 (MQ=11)
                                                                               gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACATCAACGCCTGATGCAACGCATGCCGCGTCTTATCAGGCCTACGCTAGACAGCGCAATACCCTGCTTTAGCGTGATTTTGGAGGTCGGATCAGGCGTTTaa                                                                                                                                                                      >  2:725650/1‑279 (MQ=11)
                                                                                aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                                                              >  2:485791/1‑207 (MQ=18)
                                                                                aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                                                              <  1:485791/207‑1 (MQ=18)
                                                                                  aGCGCAATAGTCTGATTTAGCGTTATTTTTTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCg                                                                                                                                                                  <  1:842920/281‑1 (MQ=25)
                                                                                    cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                                                                                   >  2:993093/1‑214 (MQ=18)
                                                                                    cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                                                                                   <  1:993093/214‑1 (MQ=18)
                                                                                               gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGCGATTTTGTAGGTCGGATACGGCGTTTATGCCGCATCCGACAtaca                                                                                                                                                     >  1:577326/1‑278 (MQ=255)
                                                                                                        tGAGTTTGTTGGTCGGTTAAGGCGTGTATGCCGCATCCGAAATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATg                                                                                                                                            <  1:638336/281‑1 (MQ=255)
                                                                                                          aTTTTGTAGGTCGGTTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGACGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCACACAGCGCAATAGCCTGATTTAGCGTCATTTTGTAGGTCGGATAAGGCGTTTATGCCTCATCCGACATCATCGCCTGATGc                                                                                                                                           >  1:645989/1‑280 (MQ=255)
                                                                                                                   cgcGGATATGGCGTTTAGCCCCCATCCCACATCAACGCCTGTTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTaa                                                                                                                                  <  1:636627/278‑1 (MQ=255)
                                                                                                                       gATAAGGCGTTTATGCCCGAGCCGACATCAACGCCGGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGt                                                                                                                             <  1:358258/281‑1 (MQ=255)
                                                                                                                        aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAACCAGCGCCATAGCCTGATTTAGCGTGATTTTTTAGGTCGGTACAGGGGTTTATTCCGCATCCGACATCAACGCGTTGTTCGACGCttgatgcaa                                                                                                                             >  2:157443/1‑273 (MQ=11)
                                                                                                                            tgtataagcgcccgaTCCGAAATCTACGCCTGATGCGACGCTTAACGCGTCTTATCGGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGAATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAt                                                                                                                        <  1:542522/267‑1 (MQ=11)
                                                                                                                                 ttgagCCGCATCCGGCATCGACGACGTATGCGAGCTGTAACGGGTATTATCAGGGCTAAGCCAGAAAGCGAAATAGCCTGATTTCGGGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATACGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                   <  2:847049/277‑1 (MQ=11)
                                                                                                                                 tgtgcccGCGTCCTACATCAACGCCTGTTGCGACGCTTACCTCGTTTTATCAGGCCTAGGCCGGAAAGCGCATTAGCCTGATTTAGCGTGATTTTGAAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCTTCTTATCAGGCCTACGCCAGACAGCGCAGTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                   <  1:755062/276‑1 (MQ=11)
                                                                                                                                  tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCCCAGACAGCGTAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACCCCTGATGCGACGCCTACCGCGTCTTATCAGGcc                                                                                                                  >  2:879454/1‑281 (MQ=255)
                                                                                                                                   aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                                                                                                                                              >  2:884401/1‑188 (MQ=35)
                                                                                                                                   aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                                                                                                                                              <  1:884401/188‑1 (MQ=35)
                                                                                                                                       cgcagctcaCATCACTGCCGGATGCTACGTTTGAGGCTGCTTACCAGGCCTACGCGCAACAGCGCAGTGGCCGTATTTAGCGTGATTTTGTAGGTCGAGTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACg                                                                                                              <  2:180892/272‑1 (MQ=11)
                                                                                                                                         cATCCGATATCAACGCCTGTGTGCGACCTTAGCGCGTCTTTTCAGGCCTACTCCATACAGCCCAATAGCCTTATTTAGCAGGATTTTTTAGGTCTGATAAGGCGTGTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCa                                                                                                           <  2:368884/281‑1 (MQ=11)
                                                                                                                                                     aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGTGCAATACCCTGATTGAGCGTGATTTTGTAGGTCGGCTAAGGCGTTTATGCCGGCTGCGACATCAAAGCCTGATGCGACGCTTAAAGCGGCTTATCAGGCCTACGctactcccttatcca                                                                                               >  2:542522/1‑267 (MQ=11)
                                                                                                                                                        ccTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGTCATCAACGCCTGATTCGACGCTTAACGACTCTTATCCAGCCTACGCCAGACAggtaaacagcc                                                                                            >  2:1009044/1‑271 (MQ=255)
                                                                                                                                                         cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCAGCGTCTTATCAGTCCTACGACAACCAGCGCAATATCCTGATTTAGCGTGATTTTGTACGTCGGATACGGCGTTTTTGCCGCATCCGACATAAGCCCCTGATTCGTCCCTTAAAGCGTCTTAACAGCCCTACTCAAGACCGCGCATTaccgt                                                                                           >  2:111650/1‑277 (MQ=11)
                                                                                                                                                                                 ttAGCAGGCCTCACCCGGACAGCGTAATGGGCTAGTTTGGCGTTATTGTGTAGGTCGGCTAAGGCGTTTACGCCGCATCGTCCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATGGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGt                                                                     <  2:537471/281‑1 (MQ=11)
                                                                                                                                                                                             cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                                                              <  1:34752/100‑1 (MQ=255)
                                                                                                                                                                                             cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgc                                                                                                                                                                                                                                                >  2:187830/1‑98 (MQ=39)
                                                                                                                                                                                             cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgc                                                                                                                                                                                                                                                <  1:187830/98‑1 (MQ=39)
                                                                                                                                                                                             cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                                                              >  2:34752/1‑100 (MQ=40)
                                                                                                                                                                                             cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAAGGCCTGATGCGACGGTTACCGCGTCTTATCAGGCCTACGCCAACCAGCGGAATATCTCGACTAAGCTGGATCTTGTAGGTTGGATAAGGCg                                                          >  2:251108/1‑280 (MQ=11)
                                                                                                                                                                                                 aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                                                                                                                                          <  2:641825/116‑1 (MQ=39)
                                                                                                                                                                                                 aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                                                                                                                                          >  1:641825/1‑116 (MQ=39)
                                                                                                                                                                                                   aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTCATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCTGACATCAACGCCTGATGTGACGCTTAACGCTTCTTATCAGGCCTACGCAACGCATGGCAATAGCCTGATTAATCTTGCTTTTGTTGCTCGGACAAGGCGTTTAAGc                                                   >  2:831400/1‑281 (MQ=11)
                                                                                                                                                                                                       cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                                                                                                                                         <  1:241391/95‑1 (MQ=255)
                                                                                                                                                                                                       cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                                                                                                                                         >  2:241391/1‑95 (MQ=40)
                                                                                                                                                                                                       cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                                    >  1:875714/1‑100 (MQ=40)
                                                                                                                                                                                                       cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCTGCCGCCTTCGACATCAACGCCTGATGCGACGCTTCCCGCGTCCTATCAGGCCGACGTCCGACAGCGGCATAGCCTCATTTAGCGTGCTTTCGTAGCTAGGATACGGTGttgcccgcga                                                >  2:84423/1‑272 (MQ=11)
                                                                                                                                                                                                                 tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                                                              <  2:978537/80‑1 (MQ=255)
                                                                                                                                                                                                                 tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                                                              >  1:978537/1‑80 (MQ=255)
                                                                                                                                                                                                                 tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                         >  1:719791/1‑101 (MQ=40)
                                                                                                                                                                                                                  gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAtgttgt                                                                         <  1:792409/244‑1 (MQ=18)
                                                                                                                                                                                                                  gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGt                                                                         >  2:792409/1‑244 (MQ=33)
                                                                                                                                                                                                                             aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTCTACGCCGCATCCGACATCAATGCCTGATGc                          >  1:562850/1‑280 (MQ=21)
                                                                                                                                                                                                                                  gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGc                     <  1:902340/280‑1 (MQ=21)
                                                                                                                                                                                                                                        cGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                                    <  1:205882/67‑1 (MQ=255)
                                                                                                                                                                                                                                        cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                                    >  2:205882/1‑67 (MQ=255)
                                                                                                                                                                                                                                        cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCtt                   >  2:902340/1‑276 (MQ=21)
                                                                                                                                                                                                                                                gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCAGACATCAAAGCCTGATGCGACGCTTGCCGCCACTTATc      >  1:537471/1‑281 (MQ=255)
                                                                                                                                                                                                                                                   tttACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGc                                                                                                                                                                                                                                                   <  1:589887/41‑1 (MQ=255)
                                                                                                                                                                                                                                                   tttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGc                                                                                                                                                                                                                                                   >  2:589887/1‑41 (MQ=255)
                                                                                                                                                                                                                                                     tACGCCGCATCCGACATCAACTCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATAAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACCCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGACCGCTACGCCGTTTACGCCGCATCCGACCTCACAGCCTCATGCAACTTCTCCCGCGTCTTATCAgac  >  2:246928/1‑278 (MQ=11)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                   
GAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGC  >  NC_000913/4295797‑4296319

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