Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,268 T→C 22.3% intergenic (+474/+168) 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,2680TC22.3% 104.9 / 24.9 57intergenic (+474/+168)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base T (16/28);  new base C (8/5);  total (24/33)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

CGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCAC  >  NC_000913/4295994‑4296551
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                               
cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACCTCAGACAGCGGAAAAGCATGATATACCATCATTCTTGT‑GGTCTGTAAAGGCTTTTATGccatctccgcgc                                                                                                                                                                                                                                                                                    >  2:36214/1‑278 (MQ=11)
    atggCCTGATTTGGCGTGATTTTGGAGGTCGTATAGGGCGTTTACGCCGCATGAGACATCAACGCCGGATGGGACGCTTGCCGCGTCTTATCAGGCATACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaat                                                                                                                                                                                                                                                                               <  2:1167191/286‑2 (MQ=11)
            aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTCGCGTGATT‑TTGTAGGTCCGATAAGGCTTTTCCGCCGCATCCGACATCAATGCCTGAt                                                                                                                                                                                                                                                                        >  2:1011338/1‑288 (MQ=255)
            aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                                                                                                                                                                                                                                        >  1:759616/1‑288 (MQ=21)
                 gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCATCCCACATACATGCCTGATGCTA‑‑c                                                                                                                                                                                                                                                                 >  2:847449/1‑288 (MQ=255)
                  cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCCTCTTATCAGGCCTCCGCCAGACCGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACGCCACTGCCTGATTCGA‑‑c                                                                                                                                                                                                                                                                 >  2:509491/1‑287 (MQ=11)
                  cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGA‑‑c                                                                                                                                                                                                                                                                 >  1:154238/1‑287 (MQ=21)
                      aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGA‑‑CGCttg                                                                                                                                                                                                                                                            >  1:305425/1‑287 (MQ=21)
                                tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGTCGCTTAACGCGTCTTATCAGGCCTACGCCAGCCACCGCAATAGCCTGCTTTAGGGTGATT‑CTGAAGGTCGGGTAAGGCGTTTACGCAGCATCCGACATTAATGCCTGATGCGA‑‑CGGttgccctgtcctc                                                                                                                                                                                                                                                  >  2:743066/1‑277 (MQ=11)
                                 cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGA‑‑CGCTTGCCGCGTCTTAt                                                                                                                                                                                                                                                 <  1:427727/288‑1 (MQ=17)
                                  ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑ATGTAGGTAGGATAAGGCGTTTACGCCGCATCCGACATAAATTCCGGATGCGA‑‑CGTTTTCCGCTGCTTATCa                                                                                                                                                                                                                                               >  2:808222/1‑289 (MQ=11)
                                     tAAGGCGTTTACGCCGAATCAGAAATCAACGCCGGATGCGACGCTTGCCGCGTCGTATCAGGCATAAGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGA‑‑CGCTTGCCGCGTCTTATCAgg                                                                                                                                                                                                                                             <  1:1226428/288‑1 (MQ=11)
                                      aaGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGA‑‑CGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                                                                                           <  1:936687/289‑1 (MQ=17)
                                                                                                       cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑cgcc                                                                                                                                                                         >  2:1090268/1‑285 (MQ=14)
                                                                                                                                  gCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGa                                                                                                                                                                                                                                                                    <  2:1095350/174‑1 (MQ=31)
                                                                                                                                                                  gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                                                               <  1:1090268/289‑1 (MQ=25)
                                                                                                                                                                    cGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                                                             >  1:936073/1‑289 (MQ=12)
                                                                                                                                                                                   cGCTTGATGCGTCCCTTTCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                                                              <  1:352631/289‑1 (MQ=255)
                                                                                                                                                                                   cGCCTGATGCGCCGCTTGCCGCTTCTTTTCAGGCCTACGCGAGACTGCGAAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                                                              <  1:850573/289‑1 (MQ=11)
                                                                                                                                                                                   cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTTATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                                                              <  1:847449/289‑1 (MQ=255)
                                                                                                                                                                                         aTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                                                         >  1:662000/1‑288 (MQ=255)
                                                                                                                                                                                         aTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                                     <  2:493686/246‑1 (MQ=2)
                                                                                                                                                                                         aTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                                     >  1:493686/1‑246 (MQ=2)
                                                                                                                                                                                              actcttgccgcgTCGTATCAGGCCAACGCCAGACAGAGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGGAGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCATCACAcc                                                                                   <  1:467/281‑1 (MQ=11)
                                                                                                                                                                                              aCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACCTAAATGCATGATGCGATGCTCGCCGCTTCTTCTCAGGCCTAGCTTAACCGTTGGTTACTCTTTCACACTCTTcaact                                                                                    >  2:491037/1‑285 (MQ=11)
                                                                                                                                                                                               cGCTTGCCGCGTCTTATGAGGCTTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                                                                  <  1:431307/289‑1 (MQ=255)
                                                                                                                                                                                               cGCTTGCCGCGTCTTATCGGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TGGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                                   <  1:376316/288‑1 (MQ=255)
                                                                                                                                                                                               cGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                                   <  1:36214/288‑1 (MQ=255)
                                                                                                                                                                                                    gccgcgTCTTATCAGGCCTACCCCAGTCGTCTCAATAGCCTTATTGAACTGGATT‑TTGTAGGTCGGATAAGCCGTTTATTCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt                                                                             <  1:875207/287‑1 (MQ=14)
                                                                                                                                                                                                         gTCTTATCAGGCCTAGTCCAGGCAGCGCAATTGCCTGTTTTTGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCGAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                                                         <  1:1195167/288‑1 (MQ=14)
                                                                                                                                                                                                         gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                                                         <  1:87626/288‑1 (MQ=255)
                                                                                                                                                                                                         gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                                                         <  1:122161/288‑1 (MQ=255)
                                                                                                                                                                                                         ctatttgcAGGCCTAGGCAAGCCAGGGAAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCCACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                                                         <  1:487886/281‑1 (MQ=14)
                                                                                                                                                                                                             tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                                                                    <  1:1016421/289‑1 (MQ=255)
                                                                                                                                                                                                                   gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                                                               <  1:1094861/288‑1 (MQ=255)
                                                                                                                                                                                                                   gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                                                               <  1:357071/288‑1 (MQ=255)
                                                                                                                                                                                                                        cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCCTCCGACATCAATGCCTGCTGCGACGCTTGCCGCGTCTTATCCGGCCTACCTTAACCGTTGGTTAATTTTTCAAACTCTTCCAACCCGGTAACAAAACCTTCCATACAGCt                                                          >  2:1181061/1‑288 (MQ=14)
                                                                                                                                                                                                                             gACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                                                              <  1:462028/247‑1 (MQ=255)
                                                                                                                                                                                                                             gACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                                                              >  2:462028/1‑247 (MQ=255)
                                                                                                                                                                                                                             gACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                                                              <  1:1259869/263‑1 (MQ=255)
                                                                                                                                                                                                                             gACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                                                              >  2:1259869/1‑263 (MQ=255)
                                                                                                                                                                                                                                  cgcAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAATTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc                                                >  1:30390/1‑288 (MQ=255)
                                                                                                                                                                                                                                  cgcAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAcc                                               >  2:842709/1‑289 (MQ=255)
                                                                                                                                                                                                                                   gcAATAGCCTGATTTACCGTGGTT‑TTGTAGGTCGGATAAGGCTTTTATGCCGCATCCGACATCACCGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc                                                <  1:1181061/287‑1 (MQ=255)
                                                                                                                                                                                                                                          ccTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                                   <  1:650993/245‑1 (MQ=255)
                                                                                                                                                                                                                                          ccTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                                   >  2:650993/1‑245 (MQ=255)
                                                                                                                                                                                                                                             gATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACCTCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCCTCACACCCGGTACCACAACCTTCCATACAGCTCAGGTTTAACCAGTGCAgcgc                                     >  1:849972/1‑288 (MQ=255)
                                                                                                                                                                                                                                                    cGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCt                             >  2:824025/1‑289 (MQ=255)
                                                                                                                                                                                                                                                              gTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTtcttct                    <  1:1028107/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTg                           <  1:112073/277‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTg                           >  2:112073/1‑277 (MQ=255)
                                                                                                                                                                                                                                                                              gTTGATGCCGCATCCGACATCAACCCCGGAGGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCTTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttctt     <  1:1151750/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                               tttATGCCGCATCCGACATCAACTGCTGGTGCGA‑‑CGCTTAACGCGTGTTATCAGGCCTAGGCCAGACAGGGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAAATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCa   <  2:179191/289‑1 (MQ=21)
                                                                                                                                                                                                                                                                                ttATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCAc  <  2:154238/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                ttACGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTGCCGCGTCTCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCAc  <  1:884237/289‑1 (MQ=18)
                                                                                                                                                                                                                                                                                 tATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCa   >  2:303762/1‑287 (MQ=255)
                                                                                                                                                                                                                                                                                 tATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCa   <  1:303762/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                               
CGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGA‑‑CGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCAC  >  NC_000913/4295994‑4296551

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