Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,060 C→T 16.4% 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.4% 491.0 / 80.6 208intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (93/80);  new base T (19/15);  total (113/95)
Fisher's exact test for biased strand distribution p-value = 8.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

Only 100 of 213 total aligned reads displayed.

AGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCC  >  NC_000913/4295799‑4296320
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                    
         gggATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgc                                                                                                                                                                                                                                                 >  4:413791/1‑274 (MQ=255)
                 gacccctcAATGTGTGAATGCCTGATGCGGCGCTTGCCGCGTCTTATCAGGGCTAAGCCAGACAGCGCGATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATAAGCGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                                                                                         <  2:427177/267‑1 (MQ=255)
                     cGGCAATGGGGGAAAGCGTGGTGGGAAGCTTGCAGCGTCTTATCAGGCCTAGGGCAGACGGCGAAATAGCCTGATTTAGCGTGAATTTGTAGGTCGGATAAGGCGTTTATGCCGCAGCCGACATCAAGGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                                     <  4:506768/274‑1 (MQ=11)
                              gtgAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCaaa                                                                                                                                                                                                                           >  2:53657/1‑273 (MQ=25)
                              gggAATGCCGTATGCGACGCGTGCCGCGTCTTATCAGGCCTACGCGAGACAGCGCAATAGCATGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAAGGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                                                                                                                                           <  2:374718/273‑1 (MQ=255)
                                       tGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAAt                                                                                                                                                                                                                  >  2:523791/1‑275 (MQ=33)
                                                gCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAtt                                                                                                                                                                                                         <  2:412685/275‑1 (MQ=33)
                                                     ccGGGTTTTATAAGGACGACTCCAGACAGCGGAGTGGCCGGTTTTAGCGTGAATTTGTGGGTCGGGTAAGTCGTTTATGCCCCATCCGACATCAACGCCTGATGCGACGCTTACAGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCg                                                                                                                                                                                                    <  4:650981/275‑1 (MQ=11)
                                                                        gtccAGACAGAGCGATAGCTGGATTGAGCGTGATTTAGTGGGTGGTGTATGGCGTTTATTCCGCATCAGAAATAAACGCTTGATGCGAGGCTTACGGC‑GTCTTATCAGGGTTACGCCAGAGGGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCACCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAt                                                                                                                                                                                  <  4:545151/272‑1 (MQ=11)
                                                                                                                 gTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAc                                                                                                                                                                                                                                                                      >  6:137972/1‑149 (MQ=32)
                                                                                                                  tCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAc                                                                                                                                                                                                                                                                      >  8:775796/1‑148 (MQ=32)
                                                                                                                   cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                                                                                                                                          >  1:96990/1‑175 (MQ=21)
                                                                                                                    ggATAAGGCGTTTATGCAGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCCCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGc                                                                                                                                                                                                                                                                    <  8:217553/148‑1 (MQ=17)
                                                                                                                        aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                                                                                                                                                                                                                               <  6:1034968/149‑1 (MQ=17)
                                                                                                                        aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                                                                                                                                                                                                                               <  6:620487/149‑1 (MQ=32)
                                                                                                                          ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc                                                                                                                                                                                                                                                             >  6:476490/1‑149 (MQ=32)
                                                                                                                           gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc                                                                                                                                                                                                                                                             >  7:85111/1‑148 (MQ=32)
                                                                                                                           gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc                                                                                                                                                                                                                                                             >  7:220730/1‑148 (MQ=32)
                                                                                                                                tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGTATAAAGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACCCTTAGCGCGTCTTAt                                                                                                                         >  4:123736/1‑275 (MQ=255)
                                                                                                                                tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTNATGCGACGCTTGCCGCGTCTTATCCGGCCTACGCCAGCCACCGCAATAGCTTGATTTAGCGTTACTTTGTAGGTNGGATCCGGCGTTTACGCCTCATCTCACATCNNACCCTTATCCGCCGCTTacccccctctgtca                                                                                                                       >  3:98257/1‑264 (MQ=11)
                                                                                                                                  tgcccgacccgcaaTCAACGCCTGAGGCGAAGTTTACCGC‑GCCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                                        <  2:570790/261‑1 (MQ=11)
                                                                                                                                    ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcc                                                                                                                                                                                                                                                   >  5:675677/1‑149 (MQ=35)
                                                                                                                                     gcactcctgaaTCAATGCCGAATGGGACGCTT‑ACGCTGTCTTATCAGGCCTACGCCAGAAGGAGAAATGGCCTGATTTAGCGTGATTTTGTGGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                    <  4:608300/265‑1 (MQ=11)
                                                                                                                                      gCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgc                                                                                                                                                                                                                                                 <  5:603216/149‑1 (MQ=25)
                                                                                                                                        aTCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                                                               >  6:426828/1‑149 (MQ=21)
                                                                                                                                             aCATCAACGACTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGAATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                                                                                                                                           >  8:688454/1‑148 (MQ=255)
                                                                                                                                               aTCAACGCCTGATGCGACGCTTACCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                                                                                                                                        <  6:251566/149‑1 (MQ=21)
                                                                                                                                                tCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                                                                                                                                       >  6:562182/1‑149 (MQ=38)
                                                                                                                                                 cAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                                                                                                                                       >  7:546468/1‑148 (MQ=38)
                                                                                                                                                  aaCGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                                                                                                                                      <  8:222362/148‑1 (MQ=39)
                                                                                                                                                   aCGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                                                                                                                                   >  7:736374/1‑118 (MQ=18)
                                                                                                                                                   aCGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                                                                                                                                   <  8:736374/118‑1 (MQ=18)
                                                                                                                                                       cTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                                                                                                                                                >  5:846481/1‑147 (MQ=17)
                                                                                                                                                        tGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                                                                                                                                                >  7:232896/1‑146 (MQ=18)
                                                                                                                                                        ctgataagacgcgttaagcgtcttatcaggcCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                                                                                                                                              <  6:407827‑M2/104‑1 (MQ=255)
                                                                                                                                                        ctgataagacgcgttaagcgtcttatcaggcCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                                                                                                                                              >  5:407827‑M2/32‑135 (MQ=255)
                                                                                                                                                          aTGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGcc                                                                                                                                                                                                                              <  7:507628/148‑1 (MQ=39)
                                                                                                                                                          aTGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGcc                                                                                                                                                                                                                              <  8:841512/148‑1 (MQ=38)
                                                                                                                                                              gACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAc                                                                                                                                                                                                                                                                      >  6:903695/1‑104 (MQ=32)
                                                                                                                                                              gACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACa                                                                                                                                                                                                                         <  5:562182/149‑1 (MQ=38)
                                                                                                                                                                    tAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCa                                                                                                                                                                                                                    <  7:456645/148‑1 (MQ=38)
                                                                                                                                                                     aaCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCCGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgc                                                                                                                                                                                                                       >  2:170167/1‑144 (MQ=37)
                                                                                                                                                                         c‑gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCaa                                                                                                                                                                                                                   <  4:170167/144‑1 (MQ=37)
                                                                                                                                                                           gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCt                                                                                                                                                                                                             >  6:536401/1‑149 (MQ=38)
                                                                                                                                                                                 tcaggcCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                                                                                                                                       >  5:875185‑M2/7‑149 (MQ=255)
                                                                                                                                                                                 tcaggcCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                                                                                                                                       <  6:675677‑M2/143‑1 (MQ=255)
                                                                                                                                                                                          cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                                                                                                                                                                                                                                 >  2:270280/1‑114 (MQ=39)
                                                                                                                                                                                          cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAtt                                                                                                                                                                                               <  7:910332/148‑1 (MQ=38)
                                                                                                                                                                                          cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGAt                                                                 >  2:637913/1‑274 (MQ=18)
                                                                                                                                                                                             cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTg                                                                                                                                                                                            <  8:379805/148‑1 (MQ=37)
                                                                                                                                                                                                aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAg                                                                                                                                                                                         >  7:105351/1‑148 (MQ=39)
                                                                                                                                                                                                  aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTc                                                                                                                                                                                      <  5:700480/149‑1 (MQ=38)
                                                                                                                                                                                                  aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTc                                                                                                                                                                                      <  6:271851/149‑1 (MQ=37)
                                                                                                                                                                                                    cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCgg                                                                                                                                                                                    >  5:100460/1‑149 (MQ=38)
                                                                                                                                                                                                    cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCg                                                                                                                                                                                     >  8:86907/1‑148 (MQ=38)
                                                                                                                                                                                                    cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCCACGCTTAACGCGTCGTATCAGGCCTACGCCAGACAGCGCAAAAGCCTGAATTAGCGTGATTTTGAGGGTCGGATAAGGCGTTTa                                                       >  4:359486/1‑274 (MQ=16)
                                                                                                                                                                                                            ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCg                                                                                                                                                                            >  6:490651/1‑149 (MQ=38)
                                                                                                                                                                                                             cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                                                                                                                                  >  7:569801/1‑94 (MQ=255)
                                                                                                                                                                                                             cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                                                                                                                                  <  8:569801/94‑1 (MQ=255)
                                                                                                                                                                                                              tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                                                                                                                                             <  6:837910/66‑1 (MQ=25)
                                                                                                                                                                                                                   tAGCGTGATTTGGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATcc                                             <  2:319727/269‑1 (MQ=25)
                                                                                                                                                                                                                       gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                                                                                                                                                                                              <  5:313198/104‑1 (MQ=255)
                                                                                                                                                                                                                       gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                                                                                                                                                                                              >  6:313198/1‑104 (MQ=40)
                                                                                                                                                                                                                        tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCa                                                                                                                                                                 >  7:317896/1‑148 (MQ=38)
                                                                                                                                                                                                                          aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCTGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                 >  2:201378/1‑274 (MQ=21)
                                                                                                                                                                                                                          aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGACTGATGAGACGCTTAACGCGTCTTATCAGGCATACGCAAGAAAGCGCAATAGCATGATTTCGCGTGATTTTGTAGGTCGGATAAGGCGTTTACTCTGCATCCTACATCAATGCCAGAta                            >  1:32001/1‑278 (MQ=11)
                                                                                                                                                                                                                           ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCg                                                                                                                                                             >  6:38584/1‑149 (MQ=39)
                                                                                                                                                                                                                             ttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGAc                                                                                                                                                           <  6:875185/149‑1 (MQ=37)
                                                                                                                                                                                                                                  ggTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCa                                                                                                                                                       <  7:751291/148‑1 (MQ=38)
                                                                                                                                                                                                                                     cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACg                                                                                                                                                    <  8:474866/148‑1 (MQ=37)
                                                                                                                                                                                                                                     cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGc                                                                                                                                                   <  5:490651/149‑1 (MQ=37)
                                                                                                                                                                                                                                        aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTg                                                                                                                                                >  5:290608/1‑149 (MQ=38)
                                                                                                                                                                                                                                            ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAg                                                                                                                                                                                         <  6:701314/104‑1 (MQ=40)
                                                                                                                                                                                                                                             gCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                                                               <  3:167514/49‑1 (MQ=255)
                                                                                                                                                                                                                                             gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                                                               >  1:167514/1‑49 (MQ=255)
                                                                                                                                                                                                                                             gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTCTCAGGCCTACGCCAGACAGCGCAATAGCCTGAGTTAGCGTGATTTTGTTGGTCGGATAAGGCGTTTACGACCCCTCCGACCTCAATGCCTGAGTCGActcgtgccgcg              >  4:258343/1‑264 (MQ=11)
                                                                                                                                                                                                                                              cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCg                                                                                                                                           >  8:268522/1‑148 (MQ=37)
                                                                                                                                                                                                                                                tttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGAc                                                                                                                                         >  8:596642/1‑148 (MQ=37)
                                                                                                                                                                                                                                                 ttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGcc                                                                                                                                                                    <  6:916679/120‑1 (MQ=39)
                                                                                                                                                                                                                                                 ttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGcc                                                                                                                                                                    >  5:916679/1‑120 (MQ=40)
                                                                                                                                                                                                                                                  tACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                                     <  2:290019/54‑1 (MQ=9)
                                                                                                                                                                                                                                                  tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                                                                     >  4:290019/1‑54 (MQ=255)
                                                                                                                                                                                                                                                  tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAt                                                                                                                         <  4:41173/162‑1 (MQ=35)
                                                                                                                                                                                                                                                  tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCCGCCCCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCttc         >  2:126210/1‑273 (MQ=255)
                                                                                                                                                                                                                                                  tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa         >  4:140172/1‑274 (MQ=17)
                                                                                                                                                                                                                                                  tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt        >  4:683235/1‑275 (MQ=17)
                                                                                                                                                                                                                                                    cgccgcATCCGACATCAACGCCTGATTCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGGAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAAGGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc  <  1:104619/279‑1 (MQ=9)
                                                                                                                                                                                                                                                     gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTa                                                                                                                                    <  7:82016/148‑1 (MQ=34)
                                                                                                                                                                                                                                                              gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAt                                                                                                                                                                       <  8:508248/104‑1 (MQ=38)
                                                                                                                                                                                                                                                                cATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAt                                                                                                                         <  8:317896/148‑1 (MQ=34)
                                                                                                                                                                                                                                                                 aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                                        <  8:105351/148‑1 (MQ=34)
                                                                                                                                                                                                                                                                 aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCa                                                                                                                       <  5:268090/149‑1 (MQ=34)
                                                                                                                                                                                                                                                                  tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCt                                                                                                                                                 <  5:153440/122‑1 (MQ=34)
                                                                                                                                                                                                                                                                  tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCt                                                                                                                                                 >  6:153440/1‑122 (MQ=34)
                                                                                                                                                                                                                                                                  tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAg                                                                                                                      >  6:984623/1‑149 (MQ=33)
                                                                                                                                                                                                                                                                    aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAgg                                                                                                                                                                              >  5:773428/1‑91 (MQ=33)
                                                                                                                                                                                                                                                                    aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAgg                                                                                                                                                                              <  6:773428/91‑1 (MQ=33)
                                                                                                                                                                                                                                                                    aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAt                                                                                                                                              >  5:412337/1‑123 (MQ=32)
                                                                                                                                                                                                                                                                    aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAt                                                                                                                                              <  6:412337/123‑1 (MQ=32)
                                                                                                                                                                                                                                                                     aCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                    >  7:122460/1‑148 (MQ=31)
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                    
AGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCC  >  NC_000913/4295799‑4296320

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑