Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,848,300 2 bp→TT coding (376‑377/552 nt) ydjA ← putative oxidoreductase YdjA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,848,3000CT95.0% 133.8 / NA 40G126K (GGG→AAG) ydjAputative oxidoreductase YdjA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/16);  minor base G (0/2);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009131,848,3010CT100.0% 136.1 / NA 40G126K (GGG→AAG) ydjAputative oxidoreductase YdjA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/18);  total (22/18)

GATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAA  >  NC_000913/1848061‑1848535
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                          
gATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAAct                                                                                                                                                                                                                                  <  1:328900/251‑1 (MQ=255)
           gtgtCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCTACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAGGTTGATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAAGGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTg                                                                                                                                                                                                                       <  1:235324/251‑1 (MQ=255)
                tGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATcc                                                                                                                                                                                                                  >  2:278943/1‑250 (MQ=255)
                tGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCg                                                                                                                                                                                                                  >  2:438736/1‑251 (MQ=255)
                 gAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCg                                                                                                                                                                                                                  <  1:470324/250‑1 (MQ=255)
                  aaTTCCGGAGCCTTCGCTCAAAATCCAGACAATTTCGCGAAACTTATCAGCACTAAGTTACAAACGCCGTCGGGTACGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGAAGAGAAAACCGACAATTTTATCCTGCTCCCGGCAACCGAATGCTTCACGCACTACCGGACTTTCCGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAAATTCTGCCATTTGCCTCGcc                                                                                                                                                                                                                >  2:208997/1‑251 (MQ=255)
                       cTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGAc                                                                                                                                                                                                           >  1:191488/1‑251 (MQ=255)
                          cgcCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACcgc                                                                                                                                                                                                        <  1:19703/251‑1 (MQ=255)
                                     accaTCCAGACAGTTTCGCGACACTTATCAGAAATACGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGCCTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCTTCGCCATGACCGCGCATCCGGc                                                                                                                                                                                               >  2:202188/4‑249 (MQ=255)
                                       aaTCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGAc                                                                                                                                                                                           >  2:208777/1‑251 (MQ=255)
                                       aaTCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCAACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGAc                                                                                                                                                                                           >  1:51975/1‑251 (MQ=255)
                                                       cgACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAgc                                                                                                                                                                           >  1:308946/1‑251 (MQ=255)
                                                         aCAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGAAAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAAATTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAgcgc                                                                                                                                                                         >  2:248828/1‑251 (MQ=255)
                                                                  aGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACttt                                                                                                                                                                 >  1:159613/1‑250 (MQ=255)
                                                                      aTAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGAt                                                                                                                                                            <  2:437443/251‑1 (MQ=255)
                                                                                        cGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGccc                                                                                                                                          >  2:474860/1‑250 (MQ=255)
                                                                                          ggaccGGGACGTTAATCGACGTAGATGCTTTCAGCTGAGGCGTACAGAGGTAGAGAAAAACGACAATTTTATCCTGCTCACGGCAACCGAATGATTCACGCACTACCGGACTTACAGCTAATACGCCACTGCGCCAGATGCAGCCAAACTTCTGGGCAACTGCTGCCATTTGCGTAGCCATGACCGCGCATCCGGCAGACCTTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGccacc                                                                                                                                        >  2:507344/4‑251 (MQ=255)
                                                                                          ggTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGccacc                                                                                                                                        <  1:474860/251‑1 (MQ=255)
                                                                                            tCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGc                                                                                                                                                                                               >  1:53775/1‑194 (MQ=255)
                                                                                            tCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGc                                                                                                                                                                                               <  2:53775/194‑1 (MQ=255)
                                                                                                    gTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGtgatg                                                                                                                               >  2:164116/1‑250 (MQ=255)
                                                                                                        aTCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAg                                                                                                                          <  1:396777/251‑1 (MQ=255)
                                                                                                                            cTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCAt                                                                                                      >  1:473948/1‑251 (MQ=255)
                                                                                                                              gcggcgTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGGGCGCGGAACGGCGCATTa                                                                                                    <  1:78445/251‑1 (MQ=255)
                                                                                                                                        gAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCCCGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGt                                                                                          >  1:552782/1‑251 (MQ=255)
                                                                                                                                                    acagACAATTTTATCCTGCTCACGGCAAACGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGCTGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGAATTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTTCCGATAACTTTGt                                                                              >  1:487708/4‑251 (MQ=255)
                                                                                                                                                    aCCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGCTTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGt                                                                              >  1:487721/1‑251 (MQ=255)
                                                                                                                                                               tATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTAcc                                                                    <  1:16838/250‑1 (MQ=255)
                                                                                                                                                                     gCTCACGGCCACCGATTGCTTCACGCACTACCGGACGTTCAGTTAATGCGCCACGGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGgcagca                                                             <  2:308946/251‑1 (MQ=255)
                                                                                                                                                                     gCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGgcagc                                                              >  2:259725/1‑250 (MQ=255)
                                                                                                                                                                        cACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGAAGCaagc                                                          >  2:289209/1‑249 (MQ=255)
                                                                                                                                                                           ggCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGcc                                                       <  2:301416/251‑1 (MQ=255)
                                                                                                                                                                                      gCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACGTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAg                                            <  2:467573/251‑1 (MQ=255)
                                                                                                                                                                                         tCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTAc                                         >  2:189323/1‑251 (MQ=255)
                                                                                                                                                                                            cGCACTACCGGACTTTCAGTTAATGCGCCTCTGCGCCAGATGCCGCCAAACTTCTTGGCAACTGCTGCCAGTTGCAGCGCCAGGACCTCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGTCTTTATGATTCTCTTCGCATTTCGCCACCACCGTTATGATGAGCGGTGCGCGGGACGGCGCATTACGGGCTTTGGCGATAGCTTTTTCATCACTCCCGGCAGCAATCGCCCCCTGTTCCAGTACGgc                                      <  1:47030/251‑1 (MQ=255)
                                                                                                                                                                                                        cgttCAGGTAATGCGCCACTGCGCCAGATGCCGCCAAACGTCGGGGCAACGGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATGTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCACTACGGGCTGTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTc                          <  2:51975/249‑1 (MQ=255)
                                                                                                                                                                                                            cAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCgcgc                      >  1:388852/1‑251 (MQ=255)
                                                                                                                                                                                                                     cgcCACCGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCGGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGcccttcccct             <  2:457274/251‑1 (MQ=255)
                                                                                                                                                                                                                                cAGATGCCGCCAAACGTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGAGGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGGCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCGTCCCCTTCAATCACaaa  <  2:535185/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cgccAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTGCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGAGAGCTTTGtcatc                                                                          <  2:457011/172‑1 (MQ=255)
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                          
GATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAA  >  NC_000913/1848061‑1848535

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: