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: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,0600CT16.0% 128.8 / 24.2 71intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
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

GAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGC  >  NC_000913/4295797‑4296319
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                   
gAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGCCGTTTACGCCGCATCCGACATCAACGCCTGGTGCGACGCTTg                                                                                                                                                                                                                                                    >  1:758509/1‑281 (MQ=255)
   gATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcc                                                                                                                                                                                                                                                  >  1:368884/1‑280 (MQ=255)
        ccagggATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGTTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGTCAGCGCAATAGCCTGATGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCTGCATCCGACACCAACGTCTGATGCGACGCTCGCcgcgcc                                                                                                                                                                                                                                             >  1:649523/4‑278 (MQ=255)
              atatCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCTGATAAGGCGTTTACGCCGCATCCGACCTCAACGCCTGATTCGACGCTTcgcgcgcggtatca                                                                                                                                                                                                                                       >  1:206049/1‑272 (MQ=255)
                taccACCCGGCAATGTGTGAATGCCTGATCGGACGCTTGCCCCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCTGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACTCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                                    <  2:366172/279‑1 (MQ=255)
                atcaccccGGCAATGTATGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                                    <  1:305597/276‑1 (MQ=255)
                 tCCACCGTGCAATGTGTGAATGCCTGATGCGACCCTTGCCGCGTCTTATCAGGCCTACGCCCGACAGCGCAATAGCCTGATTTGGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                                    <  1:63872/280‑1 (MQ=255)
                 tCCACCCGTCAATGTGTGAATGCCTGCTGCGACGCTTGCCGCGTCTTAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATAAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATTAGGc                                                                                                                                                                                                                                    <  1:809415/280‑1 (MQ=255)
                 acgaccagagaaTGTGTTAAGGCCTGATGGGACGCTTGCCGCGTCTTACCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                                    <  1:730258/270‑1 (MQ=255)
                  ccACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCCTCCCGCGTCTTATCAGGcc                                                                                                                                                                                                                                   >  1:363943/1‑280 (MQ=255)
                       cGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGc                                                                                                                                                                                                                              >  1:180892/1‑280 (MQ=39)
                                gtgAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGACTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCAGCATCCGACAGCAACGCCTGATGCGACGCTTGCCGCGGCTTATCAGGTCTACGCCAGACAgcgc                                                                                                                                                                                                                    >  1:488362/1‑281 (MQ=255)
                                       gcggATGCGAGGCTGGCCGCGTTTTATCCGGCCTACGCCAGCAGGCGCAATCGCCTATTTTGGCGGTATTTTGGAGGTCGGGTAAGGCGTTTATGCCGCATCCGAGATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCGAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                                                                                                                                              <  2:644173/277‑1 (MQ=11)
                                       ccttggggtaCGCTTGCGGCGGTTTATCAGGCCTAACCCGGAAGGCGAAATGGCCTGATTTGGCGTGATTTTGTAGGTCGGAGAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                                                                                                                                              <  2:971371/271‑1 (MQ=11)
                                         tGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTg                                                                                                                                                                                                           >  1:287873/1‑281 (MQ=33)
                                         tGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGT‑ATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCCGATTTAGCGTGATTTTGTAGGTTGGATAAGGCGTTTACGCCGCATACGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCATACCCCAACCAGCGCAATACCCTgc                                                                                                                                                                                                          >  1:156060/1‑280 (MQ=11)
                                              cgccGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTT‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGTTTGc                                                                                                                                                                                                                                                   <  1:806648/233‑1 (MQ=12)
                                                         cgttCTTGTCAGTCCTACGCCGGAAAGCACAATATCCTGGTTTAGGGTGATTTTGTAGATCGGATAAGGCGTTGATGCCGCGTCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGGGCAATAGGCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTg                                                                                                                                                                                           <  1:627418/278‑1 (MQ=11)
                                                            tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGGTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGc                                                                                                                                                                                                                                                   <  1:534548/222‑1 (MQ=17)
                                                            tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGGTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGc                                                                                                                                                                                                                                                   >  2:534548/1‑222 (MQ=17)
                                                                          cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCACGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAACCGACATCAAGGCCTGATGCGACGCTTGCCGCGTCTTACAAGGCCTACGCCAGACCGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTGGGCTAAgggg                                                                                                                                                                           >  1:763247/1‑278 (MQ=11)
                                                                               gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACATCAACGCCTGATGCAACGCATGCCGCGTCTTATCAGGCCTACGCTAGACAGCGCAATACCCTGCTTTAGCGTGATTTTGGAGGTCGGATCAGGCGTTTaa                                                                                                                                                                      >  2:725650/1‑279 (MQ=11)
                                                                                aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                                                              >  2:485791/1‑207 (MQ=18)
                                                                                aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                                                              <  1:485791/207‑1 (MQ=18)
                                                                                  aGCGCAATAGTCTGATTTAGCGTTATTTTTTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCg                                                                                                                                                                  <  1:842920/281‑1 (MQ=25)
                                                                                    cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                                                                                   >  2:993093/1‑214 (MQ=18)
                                                                                    cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                                                                                   <  1:993093/214‑1 (MQ=18)
                                                                                               gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGCGATTTTGTAGGTCGGATACGGCGTTTATGCCGCATCCGACAtaca                                                                                                                                                     >  1:577326/1‑278 (MQ=255)
                                                                                                        tGAGTTTGTTGGTCGGTTAAGGCGTGTATGCCGCATCCGAAATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATg                                                                                                                                            <  1:638336/281‑1 (MQ=255)
                                                                                                          aTTTTGTAGGTCGGTTAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGACGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCACACAGCGCAATAGCCTGATTTAGCGTCATTTTGTAGGTCGGATAAGGCGTTTATGCCTCATCCGACATCATCGCCTGATGc                                                                                                                                           >  1:645989/1‑280 (MQ=255)
                                                                                                                   cgcGGATATGGCGTTTAGCCCCCATCCCACATCAACGCCTG‑TTGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTaa                                                                                                                                  <  1:636627/278‑1 (MQ=255)
                                                                                                                       gATAAGGCGTTTATGCCCGAGCCGACATCAACGCCGG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGt                                                                                                                             <  1:358258/281‑1 (MQ=255)
                                                                                                                        aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAACCAGCGCCATAGCCTGATTTAGCGTGATTTTTTAGGTCGGTACAGGGGTTTATTCCGCATCCGACATCAACGCGTTGTTCGACGCttgatgcaa                                                                                                                             >  2:157443/1‑273 (MQ=11)
                                                                                                                            tgtataagcgcccgaTCCGAAATCTACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCGGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGAATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAt                                                                                                                        <  1:542522/267‑1 (MQ=11)
                                                                                                                                 ttgagCCGCATCCGGCATCGACGACGT‑ATGCGA‑GCTGTAACGGGTATTATCAGGGCTAAGCCAGAAAGCGAAATAGCCTGATTTCGGGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATACGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                   <  2:847049/277‑1 (MQ=11)
                                                                                                                                 tgtgcccGCGTCCTACATCAACGCCTG‑TTGCGACGCT‑TACCTCGTTTTATCAGGCCTAGGCCGGAAAGCGCATTAGCCTGATTTAGCGTGATTTTGAAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCTTCTTATCAGGCCTACGCCAGACAGCGCAGTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                   <  1:755062/276‑1 (MQ=11)
                                                                                                                                  tATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCCCAGACAGCGTAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACCCCTGATGCGACGCCTACCGCGTCTTATCAGGcc                                                                                                                  >  2:879454/1‑281 (MQ=255)
                                                                                                                                   aTGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                                                                                                                                              >  2:884401/1‑188 (MQ=35)
                                                                                                                                   aTGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                                                                                                                                              <  1:884401/188‑1 (MQ=35)
                                                                                                                                       cgcagctcaCATCACTGCCGG‑ATGCTACGTT‑TGAGGCTGCTTACCAGGCCTACGCGCAACAGCGCAGTGGCCGTATTTAGCGTGATTTTGTAGGTCGAGTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACg                                                                                                              <  2:180892/272‑1 (MQ=11)
                                                                                                                                         cATCCGATATCAACGCCTGTGTGCGAC‑CT‑TAGCGCGTCTTTTCAGGCCTACTCCATACAGCCCAATAGCCTTATTTAGCAGGATTTTTTAGGTCTGATAAGGCGTGTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCa                                                                                                           <  2:368884/281‑1 (MQ=11)
                                                                                                                                                     aCGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGTGCAATACCCTGATTGAGCGTGATTTTGTAGGTCGGCTAAGGCGTTTATGCCGGCTGCGACATCAAAGCCTGATGCGACGCTTAAAGCGGCTTATCAGGCCTACGctactcccttatcca                                                                                               >  2:542522/1‑267 (MQ=11)
                                                                                                                                                        ccTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGTCATCAACGCCTGATTCGACGCTTAACGACTCTTATCCAGCCTACGCCAGACAggtaaacagcc                                                                                            >  2:1009044/1‑271 (MQ=255)
                                                                                                                                                         cTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCAGCGTCTTATCAGTCCTACGACAACCAGCGCAATATCCTGATTTAGCGTGATTTTGTACGTCGGATACGGCGTTTTTGCCGCATCCGACATAAGCCCCTGATTCGTCCCTTAAAGCGTCTTAACAGCCCTACTCAAGACCGCGCATTaccgt                                                                                           >  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)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                   
GAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTG‑ATGCGACGCT‑TAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGC  >  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: ‑