Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,434,8820TG39.5% 58.2 / 13.6 39Y252S (TAC→TCC) truAtRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (4/19);  new base G (15/0);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 5.01e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTGCGCCACGGGGTTC  >  NC_000913/2434641‑2435122
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
aaCAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑GCAGGGGAATCCAccc                                                                                                                                                                                                                                           >  1:477553/1‑250 (MQ=255)
             aaaCCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGAAAATCATA‑CCCG‑GCAGGGGAATACAcccccacccaaacaa                                                                                                                                                                                                                               >  2:77597/1‑237 (MQ=255)
                          gTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGAGCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCtttt                                                                                                                                                                                                                 <  1:45316/251‑1 (MQ=255)
                                         ctctGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAAATCAAA‑CCGG‑TCAGGGGAATCCAACGCGACCAGATACAGCCCCTCCCCTTTTTCCGTTGCTgccgcc                                                                                                                                                                                                  >  1:109396/1‑251 (MQ=255)
                                                          cAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTgc                                                                                                                                                                                 <  1:77597/251‑1 (MQ=255)
                                                              gtgCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑CCGG‑GCAGGGGAATCCAACCCGACCAGAAACAGCCCCTCCCCCTTTTCCGTTGGTTCCCCCCGCGGTCTGTCCTTTtccgcc                                                                                                                                                                             >  1:268375/1‑251 (MQ=255)
                                                              gtgCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTgccgcc                                                                                                                                                                             <  1:400262/251‑1 (MQ=255)
                                                                        aaaTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCCG‑GCAGGGGAAACCACCGCGACCAGATACAACCCCTCCGCTTTTGCCGGTGCTGCCCCCCGCGGTCTGTCCTTTTCCGCCAGC‑CActctt                                                                                                                                                                  >  1:379335/1‑250 (MQ=255)
                                                                          aTCGATGAGGAAGTAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCAAA‑CCCG‑TCCGGGGAATCCCCCGCGACCAGAAACAACCCCTCCCCTTTTGCCGTTGCTTCCGCCACCGTTTTGTCCTTTTCCGCCCGC‑AAATCTGCt                                                                                                                                                                >  1:434029/1‑251 (MQ=255)
                                                                                 aGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACCCTGCTATCCAGc                                                                                                                                                         <  1:418567/251‑1 (MQ=255)
                                                                                      ataaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTcc                                                                                                                                                    <  2:336977/251‑1 (MQ=255)
                                                                                              ggTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCAAA‑CCCG‑TCAGGGTAAACCACCGCGACCAGATACAACCCTTCCGCTTTTGCCCTTGCTGCCGCCCGCGATCTGGCCTTTTCCCCGAAC‑AACCCTGCTATCCAGCCCTCCCGCTGGtt                                                                                                                                            >  1:236019/1‑251 (MQ=255)
                                                                                                     atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGTCTGGTTGTGGG‑c                                                                                                                                     <  2:502679/250‑1 (MQ=255)
                                                                                                      tttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATA‑CCGG‑GCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCCTTTTCCGGTGCTGCCGCCAGCGGTCTGTCCTTTGCCGCCCGC‑AACTCCGCTATCCCGCTCTCCGGCTGGGTGTGGG‑CAc                                                                                                                                   >  2:42074/1‑251 (MQ=255)
                                                                                                             aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTc                                                                                                                            <  2:365087/251‑1 (MQ=255)
                                                                                                                 cTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCCG‑TCAGGGGAATCCACCGCGACCAGATACAACCCCTCCCCCTTTGCCGTTGCTGCCCCCCGCGTTCTTTCCTTTTCCGCCCGC‑AACTCTCCTATCCAGCTCTCCGGCTGGTTGTGGG‑GACCGACTTCCATc                                                                                                                        >  1:215275/1‑251 (MQ=255)
                                                                                                                         cGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGcc                                                                                                                <  1:587637/251‑1 (MQ=255)
                                                                                                                             aTCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACATA‑CCCG‑GCCGGGGAATCCACCGCGACCAGATACAGCCCTTCCCCTTTTTCCCTTGCTGCCGCCCGCGGTCTGTCCTTTTCCGCCCGC‑AACTCTGCTTTCCAGCTCTCCCGCTGGTTGTGGG‑CACCGACTTccctccggcctcccac                                                                                                             >  1:245933/1‑236 (MQ=255)
                                                                                                                                  tCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGCCGACAatat                                                                                                        <  2:122591/250‑1 (MQ=255)
                                                                                                                                   cAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGGAATCCACCGCGACCAGATACAACCCTTCCGCTTTTGCCGTTGCTGCCCCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGGTGTGGG‑CACCGACTTCCATCAAGCTGCCGACAATATTc                                                                                                      >  1:348032/1‑251 (MQ=255)
                                                                                                                                    aGGCTGTAACTACGTTGGCTACAGCCTGTAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGCTTTGGAAGATCATA‑CCGGGCCAGGG‑AATACCCCCCGAACAGATAAAAAACTTCCGCTTTTTCCGGTGCTGCCCCCCACGTTATGTCCTTTGCCGCCAGC‑CAATCTTCTCTTCAGCCCTCCCGCCGGGTGTGGGTCACC‑ACTTCCATCAGGGTGCCCACaaccttcc                                                                                                     >  1:489280/1‑245 (MQ=255)
                                                                                                                                      gCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑CCGG‑GCAGGGGAAACCACCCCGCCCAGATACAACCCTTCCCCTTTTTCCGTTGCTGCCCCCAGCGTTCTTTCCTTTTCCGCCCGCCAAC‑CTTCTATCCAGCTCTCCGGCCGGGTGTGGG‑CACCGACTTCCCTCAGGCTGCCCACAATATTCCTg                                                                                                   >  1:564284/1‑251 (MQ=255)
                                                                                                                                                             cTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATA‑ACGG‑TCAGGGGAATCCACCCCGACCCGATACAGCCCTTTCCCTTTTTCCGTTGCTTCCGCCCGCGCTCTTTTCCTTGTCGCCCGC‑CACCCTGCTATCCCGCTCTCCCGCTGGGTGTGGG‑GACCGACTTCCCTCCGGGCGCCCGCAAAATTCCTGACCATatgagggaacaaaacaa                                                                             >  2:214710/1‑237 (MQ=255)
                                                                                                                                                               gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCAtt                                                                            <  1:42074/249‑1 (MQ=255)
                                                                                                                                                                         gAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATTCTA‑CCGG‑TAAGGGGAATCCACCCCGCACAGATACAAGCCTTCCGCCTTTTCCGTTGCTTCCCCCAGCGTTCTGTCCTTTTCCGCCAGC‑AAACCCGGTCTCCAGCTCTCCCGCTTGGTGGGGG‑GACCCACTTCCATCAGGGCGCCGACAAAATTCCTGGCCCTATTATGTACAAAAGAATTCGCTTTAATATc                                                                >  2:227385/1‑251 (MQ=255)
                                                                                                                                                                           gCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTa                                                              <  2:379335/251‑1 (MQ=255)
                                                                                                                                                                            cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGGCCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTac                                                             <  2:236019/251‑1 (MQ=255)
                                                                                                                                                                            cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTac                                                             <  2:215275/251‑1 (MQ=255)
                                                                                                                                                                                attCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCCCCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTaccac                                                          <  2:64753/250‑1 (MQ=255)
                                                                                                                                                                                         aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAAATTATAACCGG‑‑AAGGGGAATTCACCGCGAACAGAAACAGCCCCTCCGCTTTTGTCGTTGTTGTCGCCCGCGTTCTTTCCCTTGCCGCCCCC‑AACTCTTCTTTTCATATCTCCGGCTGGGTGGGGG‑GAACGAAGTCCATCAGGGCGTCGACAATATTCCTGACCATAAGATAGACAAAAACAATTGGTTTAATATCTactcacaccaaa                                                   >  2:416090/1‑238 (MQ=255)
                                                                                                                                                                                          gTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑GCAGGGGAATCCCCCCCCACCCGATACAACCCTTTCCCTTTTGCCGTTGCTGCCCCCAGCGGTCTGTCCTTTGCCGCCCGCCAAC‑CTGCTATCCAGCCCTCCGGCCGGGTGTGGG‑CACCGACTTCCCTCAGGCTGCCGACAATATTCCTGACCATATGATTGACAAAAACAATCGCCTTTATATTTACCACCACATaaaga                                                >  2:360497/1‑247 (MQ=255)
                                                                                                                                                                                                 aGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGAcg                                        <  2:245933/251‑1 (MQ=255)
                                                                                                                                                                                                                      ggTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAatat                                                                 <  1:562911/205‑1 (MQ=255)
                                                                                                                                                                                                                      ggTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAatat                                                                 >  2:562911/1‑205 (MQ=255)
                                                                                                                                                                                                                              aaGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTGCGcc           <  1:161555/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aTA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTGCGCCACgggg     >  2:537729/1‑251 (MQ=255)
                                                                                                                                                                                                                                      a‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTGCGCCACGGGGtt   >  2:581194/1‑251 (MQ=255)
                                                                                                                                                                                                                                      a‑CCGG‑TCAGGGTAACCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTGCGCCACgggg     <  2:477553/249‑1 (MQ=255)
                                                                                                                                                                                                                                        ccGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTGCGCCACGGGGTTc  <  1:214710/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
AACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGG‑TCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGC‑AACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGG‑CACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTGCGCCACGGGGTTC  >  NC_000913/2434641‑2435122

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑