Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,060 C→T 18.8% 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,0600CT18.8% 164.0 / 39.9 80intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (39/26);  new base T (5/10);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 8.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-02

TGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGAC  >  NC_000913/4295803‑4296298
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
tGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGTCGCTTCACGCGTCTTATCAGGCCTACGCTAGACAGCGCAATAGCCTGATTTAGCGTGCTTTTTTAGGCCGGATAAGGCGTTTACGCCGCATCCGACATTAACGCCTGATGCTACGcct                                                                                                                                                                                                                                 >  2:467729/1‑272 (MQ=255)
            ccACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCATCTTAt                                                                                                                                                                                                                     >  2:458740/1‑274 (MQ=255)
                 cGGGAATGTGGGAATTGCGGATGGTACGCTTGCGCGGGCTTATCAGGGCTACGGCAAGAGGCGAAATAGGCTGAGTTAGCTTGATTTTGTAGGTCGGTTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                <  4:86468/274‑1 (MQ=11)
                 cGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                <  2:340571/274‑1 (MQ=39)
                   cataggtgggAATGCCTGGTGGGAGTCTTGGCGCGTCTGTTCAGGCCTCACCCAGAAGGGGCAATAGTCTGCTTTGGCGTGATTTTGTAGGGCGGATCAGGCGTTTATGCCGCATCGGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCt                                                                                                                                                                                                              <  4:364999/269‑1 (MQ=11)
                      aTGTGTGAATGCCTGATGCGACGCTTGCGGGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACg                                                                                                                                                                                                           <  2:168173/274‑1 (MQ=38)
                          gtgCATGCTTGATGGGCCGCTGGCCGGGTCGTATCAGGCCTACGCCGGCAAGCGCAATTGCCTGATTTAGCGTGATTTTGTAGGTCGGGTAAGGCGTTTATGCCGCATCCGACATAAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                                                                                                                      <  2:353565/275‑1 (MQ=11)
                                 ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTACGCCAGACAgcgc                                                                                                                                                                                                >  2:288662/1‑274 (MQ=17)
                                   tGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                                                                                                                        >  1:18833/1‑248 (MQ=33)
                                          aCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATCAGTCGTTTACGCCGCATCCCACCTCATCGCCTTATGCAACGCTTGCCGCGTCTTCGCAGGCCTACTCCAGACCGTGCacaatccgta                                                                                                                                                                                      >  4:490423/1‑266 (MQ=11)
                                                          aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGt                                                                                                                                                                      >  2:124216/1‑275 (MQ=33)
                                                                         gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCCTCCGACATCAACGCCTGATGCGACCCTTTCCGCGTCTTATCAGGCCTACCCCCGACACCTCAACACCACGATGTCGCGTGGTTTTGTAgtcgcgtaagggac                                                                                                                                                       >  4:45481/1‑262 (MQ=11)
                                                                         gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                                                                    >  2:423357/1‑214 (MQ=34)
                                                                             gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                <  4:423356/214‑1 (MQ=34)
                                                                              cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCCTCTTATCAGGCCTACCCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTa                                                                                                                                                   >  2:55572/1‑274 (MQ=18)
                                                                               gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAt                                                                                                                                                  >  2:276684/1‑274 (MQ=33)
                                                                                   tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                                                                                                                   >  2:156067/1‑205 (MQ=18)
                                                                                       cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                                                               <  4:156066/205‑1 (MQ=18)
                                                                                         gATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGGTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCg                                                                                                                                        <  2:33531/274‑1 (MQ=37)
                                                                                                                            tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTCATTTGGTATGTCGGATACGGCGTTTATGCCGCATCCGACATCAACGCCTGTTGCGACCTTATCCGCGTCTTAt                                                                                                   >  4:165935/1‑275 (MQ=11)
                                                                                                                                               aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGCCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGACTTAGCGTGATTTTGTAGGTCGGATAAGTCGTGTATGCCGCATCCGACCTCAACGCCTGATGCGACGCTTAACGCGTCTTTGCAGCCCTACGCCAGACAgc                                                                                >  2:516237/1‑275 (MQ=255)
                                                                                                                                                cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCGGCATCCGACATCAACGCCTGATGCTACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAgcg                                                                               >  2:448728/1‑275 (MQ=25)
                                                                                                                                                    tgataagACGCGTTAAGCGTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCg                                                                                                                                                                                                                                       <  7:622091/114‑1 (MQ=11)
                                                                                                                                                                        cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                                                               <  4:274631/124‑1 (MQ=25)
                                                                                                                                                                        cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                                                               >  2:274632/1‑124 (MQ=38)
                                                                                                                                                                             cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                                                                                                                  >  6:555611/1‑116 (MQ=25)
                                                                                                                                                                             cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                                                                                                                  <  5:555611/116‑1 (MQ=37)
                                                                                                                                                                               ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                                          <  5:589282/106‑1 (MQ=255)
                                                                                                                                                                               ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                                          >  6:589282/1‑106 (MQ=255)
                                                                                                                                                                                gCCTACGCCAGACAGAGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATAAGGCCTACGCCAGCCAGAGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTATGCCGCATCGGCCATCAATGCATGGTGCGACGATTAACGCGTCTTCTCAGGCATACGCCCGACAGGACAAGAGCCTGCTTAAGTGTGATTTTGGAGGTc                                               >  1:154795/1‑275 (MQ=11)
                                                                                                                                                                                 ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                                                                                                                                                                                                           >  7:205085/1‑87 (MQ=33)
                                                                                                                                                                                 ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                                                                                                                                                                                                           <  8:205085/87‑1 (MQ=33)
                                                                                                                                                                                   tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCt                                                                                                                                                                                                                                  >  5:95247/1‑94 (MQ=21)
                                                                                                                                                                                   tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCt                                                                                                                                                                                                                                  <  6:95247/94‑1 (MQ=21)
                                                                                                                                                                                     cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACg                                                                                                                               >  4:313510/1‑191 (MQ=34)
                                                                                                                                                                                         aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTg                                                                                                                           <  2:313511/191‑1 (MQ=34)
                                                                                                                                                                                           aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                                                                                                                            >  6:913030/1‑107 (MQ=25)
                                                                                                                                                                                           aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                                                                                                                            <  5:913030/108‑2 (MQ=25)
                                                                                                                                                                                               cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                                                               <  4:397415/101‑1 (MQ=38)
                                                                                                                                                                                               cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                                                               >  2:397416/1‑101 (MQ=40)
                                                                                                                                                                                               cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                <  4:546267/100‑1 (MQ=40)
                                                                                                                                                                                               cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                  >  4:133719/1‑209 (MQ=34)
                                                                                                                                                                                                   aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                              <  2:133720/209‑1 (MQ=34)
                                                                                                                                                                                                       ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCaa                                                                                                                                                                                              >  4:569786/1‑110 (MQ=40)
                                                                                                                                                                                                          gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                                          <  2:221056/79‑1 (MQ=255)
                                                                                                                                                                                                          gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                                          >  4:221055/1‑79 (MQ=255)
                                                                                                                                                                                                          gATTTAGCGGGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCg                      <  2:333353/274‑1 (MQ=25)
                                                                                                                                                                                                           aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgcg                                                                                                                                                                                                                           <  4:39447/77‑1 (MQ=255)
                                                                                                                                                                                                           aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgcg                                                                                                                                                                                                                           >  2:39448/1‑77 (MQ=255)
                                                                                                                                                                                                                gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCACGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCCGCCGCATCCGACATCa                >  2:91841/1‑274 (MQ=25)
                                                                                                                                                                                                                 cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGTGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCACCCCTTAACGCGTCTTCTCAGGCCCACCCTTAACAGTGCTATAGCCTGATTTATCCTGATTTTGCTGGCCGGATCAGGCTTTAACGCCCCTTCCCAccacca               >  4:385928/1‑269 (MQ=11)
                                                                                                                                                                                                                     aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAc                                                                                                            >  4:539451/1‑177 (MQ=35)
                                                                                                                                                                                                                       tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGa                                                                                                                                                                                                                                           >  8:1018768/1‑49 (MQ=255)
                                                                                                                                                                                                                       tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGa                                                                                                                                                                                                                                           <  7:1018768/49‑1 (MQ=255)
                                                                                                                                                                                                                         tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGt                                                                                                        <  2:539452/177‑1 (MQ=35)
                                                                                                                                                                                                                              gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGCTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCAGCATCCGACATCAATGCCTGATGCGAc  >  2:310115/1‑274 (MQ=21)
                                                                                                                                                                                                                               tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAg                                                                                                                                                                                                   <  7:560076/81‑1 (MQ=255)
                                                                                                                                                                                                                               tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAg                                                                                                                                                                                                   >  8:560076/1‑81 (MQ=255)
                                                                                                                                                                                                                                   atgcggCGTAAACGCCTTATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATaa                                                                                                                                                           <  5:420250/113‑1 (MQ=11)
                                                                                                                                                                                                                                   atgcggCGTAAACGCCTTATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATaa                                                                                                                                                           >  6:420250/5‑117 (MQ=11)
                                                                                                                                                                                                                                   aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATaa                                                                                                                                                           >  5:132655/1‑117 (MQ=40)
                                                                                                                                                                                                                                   aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATaa                                                                                                                                                           <  6:132655/117‑1 (MQ=39)
                                                                                                                                                                                                                                         cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTg                                                                                                                                                                       <  6:760646/99‑1 (MQ=40)
                                                                                                                                                                                                                                         cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTg                                                                                                                                                                       >  5:760646/1‑99 (MQ=40)
                                                                                                                                                                                                                                         cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGt                                                                                                                                                                      >  6:332497/1‑100 (MQ=255)
                                                                                                                                                                                                                                         cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGt                                                                                                                                                                      <  5:332497/100‑1 (MQ=255)
                                                                                                                                                                                                                                           tttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                                                                                                                   >  7:952497/1‑53 (MQ=255)
                                                                                                                                                                                                                                           tttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                                                                                                                   <  8:952497/53‑1 (MQ=255)
                                                                                                                                                                                                                                           tttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGt                                                                                                                                                                      >  7:1139459/1‑98 (MQ=255)
                                                                                                                                                                                                                                           tttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGt                                                                                                                                                                      <  8:1139459/98‑1 (MQ=255)
                                                                                                                                                                                                                                             tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                                                                                                                          >  7:824459/1‑76 (MQ=255)
                                                                                                                                                                                                                                             tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                                                                                                                          <  8:824459/76‑1 (MQ=255)
                                                                                                                                                                                                                                                   gcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGttt                                  >  4:403170/1‑221 (MQ=32)
                                                                                                                                                                                                                                                      tCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATc                                                                                                                                          <  6:974378/115‑1 (MQ=37)
                                                                                                                                                                                                                                                      tCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATc                                                                                                                                          >  5:974378/1‑115 (MQ=38)
                                                                                                                                                                                                                                                       ccGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGc                              <  2:403171/221‑1 (MQ=17)
                                                                                                                                                                                                                                                          aCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgc                                                                                                                                                                                                  >  6:444400/1‑55 (MQ=255)
                                                                                                                                                                                                                                                          aCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgc                                                                                                                                                                                                  <  5:444400/55‑1 (MQ=255)
                                                                                                                                                                                                                                                             tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTc                                                                                                                                                                 <  6:798644/85‑1 (MQ=37)
                                                                                                                                                                                                                                                             tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTc                                                                                                                                                                 >  5:798644/1‑85 (MQ=37)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
TGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGAC  >  NC_000913/4295803‑4296298

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