Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,665,698	0	T	C	32.4%	65.9 / 11.9	37	D38G (GAT→GGT)	gadX	DNA‑binding transcriptional dual regulator GadX
Reads supporting (aligned to +/- strand):  ref base T (21/4);  new base C (0/12);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 9.82e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAACGCAAAGATTAATGCTCTTTTACGTAGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAACAGGATGATAATGTGTGATGTTC  >  NC_000913/3665448‑3665935                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               caacGCAAAGATTAATGCTCTTTTACGTAGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCAt                                                                                                                                                                                                                                                <  2:677582/251‑1 (MQ=255)  aacGCAAAGATTAATGCTCTTTTACGTAGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCTCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATGCAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATc                                                                                                                                                                                                                                               <  2:957589/251‑1 (MQ=255)             tAATGCTCTTTTACGTAGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAAcc                                                                                                                                                                                                                                     >  2:357228/1‑250 (MQ=255)                             aGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTaa                                                                                                                                                                                                                    >  2:584592/1‑250 (MQ=255)                                 tttCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAAt                                                                                                                                                                                                                 >  2:478929/1‑249 (MQ=255)                                          aGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTc                                                                                                                                                                                                       >  1:152787/1‑250 (MQ=255)                                                aaTTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGAAATATCTCTACTACAGTGATAAACGAGTAATGAAGAAACTTTTTTATGTGAATGAAGTTTTAAGCCCTCCTCCTGG‑AGCAGCGGGAAAAATAATAAAAGCGTGCCCC‑TTGACACAAATAAAAAATTTTC‑AACACCCTTATCTCCTCGAATTTGG‑CTTGCACCCGAAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTa                                                                                                                                                                                                <  1:709021/250‑1 (MQ=255)                                                       aTTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATgg                                                                                                                                                                                         >  1:265794/1‑251 (MQ=255)                                                         ttagcagttgtggcccatcttcgtaaaaaaCAGGAATATCTCTACTAAGGTAATGAACGATTAATGAAAAAATTTTTTTATGTGATTAAAGATTTAATGCCTCCTCTATG‑AGCATAGGAAAAAATAATAAAGGCTTCTCCC‑TTGCCAAGAATA‑AAAATTTTCCAACACCCTTATCTATACGAATTTGT‑CTTGCACCCGCAAAAACCGGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGt                                                                                                                                                                                        <  2:514058/225‑1 (MQ=255)                                                           tGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                                                                                                                                                                                       >  1:84137/1‑249 (MQ=255)                                                                  gTTGCGCACCTTCCTAAAAAACAGGAATTTCTTTATTAAAGGAATAAAGGAGAAATAAAGAAA‑TTTTTTAAGTGTAAGAAATTTTAAATCCCCCCCCTTGGAGCATGGGGAAAAATAA‑GAAAGCCGGTCCCCTTAACAAGAATAAAAAATTTCA‑AACCCATTTATCTACTAAATTTTGT‑TTTGCACCCGAAAAACCCGGGTA‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAatatat                                                                                                                                                                              <  2:84137/251‑1 (MQ=255)                                                                           attccaaaaaaaCAGAAATACCTCTAATGCAGTGATGAAAGAGAAAGGAAGAATCTTTTTTAGGTGCAAGAAGATTTAACCCCTCCCCCTGG‑AGCAGCGGGAAAAATAATGAAAGCGTGCCCC‑TTGACAAGAAAAAAAAATTTAC‑AACACACTTTTTTA‑TTAAATTTGGGTTTGCACCCAAAAAAACCGGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCtt                                                                                                                                                                     <  2:167875/244‑1 (MQ=255)                                                                             tcataaaaaaCAGAAATATCTCTACTACAGTGATGAACGAGTAATGAAAAAATTTTTTTATGTGCATGAAGATTTAATGCCCCCCCCTTG‑AGCATCGGGAAAAATAATAAAGGCGTGTCCC‑TTGACACAAATACAAAATTTTC‑AACACACTTATCTACTCGATTTTGG‑TTTCCACCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTg                                                                                                                                                                    <  1:117774/245‑1 (MQ=255)                                                                               ataaaaaaaAGGAATATCTCTACTACAGGGAAGAACGAAAAATGAAGAAACTTTTTTAGGTGAAAGAAGATTTAATGCCTCCTCCTTG‑AGCATGGGGAAAAAAAATAAAAGCTTGCCCT‑TTGACACAAATAAAAAATTTCC‑ACCCCCCTTATCTCCTCAAATTTGG‑CGTGCCCCCAAAAAAACCAGGAC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCAt                                                                                                                                                                 <  2:432521/248‑1 (MQ=255)                                                                                 gAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc                                                                                                                                                               >  1:980771/1‑251 (MQ=255)                                                                                     aaCAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCAt                                                                                                                                                                                                                                                >  2:556477/1‑167 (MQ=255)                                                                                     aaCAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCAt                                                                                                                                                                                                                                                <  1:556477/167‑1 (MQ=255)                                                                                                ctctACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTaa                                                                                                                                                                                               >  1:979681/1‑204 (MQ=255)                                                                                                ctctACTACAGTGATGAACGAGAAATGAAAAAACTTTTTTATGTGCATGAATTTTTAATGCCCCCCCCTTG‑ACAATCGGAAAAAAAAAAAAACGCGTGCCCC‑TAGACACAAATAAAAAATTTCC‑AACACCCTTATCTACACAATTTTGG‑TTTCCCCCCGAAAAACCCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTaa                                                                                                                                                                                               <  2:979681/204‑1 (MQ=255)                                                                                                              tgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTg                                                                                                                                  >  1:910255/1‑251 (MQ=255)                                                                                                                 agaagaaaTGAAAAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑ACCTTCGGGAGAAATAATGAAAGCTTGCCCC‑TTAACACAAATAAAAAATTTCC‑ACCACACTTTTCTACTCGATTTTGG‑CTTCCACCCAAAAAAACCGGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCAt                                                                                                                               <  2:980771/245‑1 (MQ=255)                                                                                                                           aagaaACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACtt                                                                                                                      >  1:22188/1‑250 (MQ=255)                                                                                                                                            gCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACatat                                                                                                    >  1:59517/1‑251 (MQ=255)                                                                                                                                                       tAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTttatt                                                                                         >  1:981163/1‑251 (MQ=255)                                                                                                                                                       tAATGCCTCCCCCTGG‑AGCATCGGGAAAAATAAGAAAAGCGTGTCCC‑TTGACACGAATAAAAAATTTCC‑AACACCTTTATCTCCTCGATTTTGG‑TTTGCACCCGAAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTttat                                                                                          <  1:105263/250‑1 (MQ=255)                                                                                                                                                        aaTGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTttatt                                                                                         >  1:812978/1‑250 (MQ=255)                                                                                                                                                               cctcctTG‑TCCATCGGGAAAAATAAAAAAGGCGTCCCCC‑TTGACACAAATAAAAAATTTCC‑ACCACTCTTATCTCCTCGATTTTGG‑CTTGCACCCGAAAAACCCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGc                                                                                 <  1:168166/251‑1 (MQ=255)                                                                                                                                                                ctcctTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGt                                                                                                                                       >  2:1006911/1‑196 (MQ=255)                                                                                                                                                                           aTCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCctc                                                                      >  1:24832/1‑251 (MQ=255)                                                                                                                                                                            gggggAAAAATAATAAAGGCGTGTCCC‑TTAAAAAAAAAAAAAAATTTCA‑AACACCCTTATCTACCGGAATTGGG‑CTGCCACCCAAAAAAACCAGTTA‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCctct                                                                     <  2:24832/249‑1 (MQ=255)                                                                                                                                                                             ggggAAAAAAAATGAAGGCGTGCCCC‑TTGCAAAAAATAAAAAATTTCA‑AACACATTTATCTCCTAAATTTTGG‑TTTGCATCC‑AAAAAACCCGGTCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCtt                                                                    <  2:910255/250‑1 (MQ=255)                                                                                                                                                                                      taatGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGACTTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATAtt                                                                                                                >  1:292850/1‑198 (MQ=255)                                                                                                                                                                                      aaatGAAAGCGTGTCCC‑TTGACCCGAATACAAAATTTTA‑AACACACTTATCTACTCGTATTTGG‑TTTGCATCCGAAAAACCCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAat                                                           <  1:357228/250‑1 (MQ=255)                                                                                                                                                                                                                aaTACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGa                                   >  2:945741/1‑250 (MQ=255)                                                                                                                                                                                                                                     cTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAACAGGATGATa               >  2:105405/1‑250 (MQ=255)                                                                                                                                                                                                                                       tATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc                                                                                                                                                               >  2:1113884/1‑104 (MQ=255)                                                                                                                                                                                                                                            acccGAATTTGG‑CTTGCATCCAAAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAACAGGATGATAAtgtgt         <  2:265794/246‑1 (MQ=255)                                                                                                                                                                                                                                            aCTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAACAGGATGATAAtgtgtg        >  1:976932/1‑250 (MQ=255)                                                                                                                                                                                                                                             cTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAACAGGATGATAATGTGTGAt      >  1:77911/1‑251 (MQ=255)                                                                                                                                                                                                                                                 aaTTTGG‑TTTGCACCCGAAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAACAGGATGATAATGTGTGATGTTc  <  2:77911/251‑1 (MQ=255)                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               CAACGCAAAGATTAATGCTCTTTTACGTAGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTG‑AGCATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTC‑ACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAACAGGATGATAATGTGTGATGTTC  >  NC_000913/3665448‑3665935

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