Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,670	T→A	33.5%	pseudogene (934/1617 nt)	glvC ←	pseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,862,670	0	T	A	33.5%	11.7 / 2.5	33	pseudogene (934/1617 nt)	glvC	pseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base T (9/8);  new base A (4/4);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAAT‑TCCCACCAGCATGGC‑GG‑TT‑AAGG‑TTGCGGGA‑ATCAGC‑AAG‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGA  >  NC_000913/3862435‑3862900                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                              aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAAAAGCACAGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAAT‑TCCCACCAGCATGGC‑GG‑TTAAAGG‑TTGGGGGA‑Ata                                                                                                                                                                                                                                              >  1:673974/1‑237 (MQ=255)           gCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGAAAATTCCAACGGTTCGGTAAT‑TTCCACCAGCATGGC‑GG‑GT‑TAGG‑TTTCGGGG‑AACAGG‑AAgacgcc                                                                                                                                                                                                                                  >  1:309457/1‑233 (MQ=255)                 tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGGGAAATGAACCGGAAGGTAAATTCCAGCGGTTCGGTAAT‑TCCCCCCAGCATGGC‑GG‑GT‑AAGGATTGGGGGA‑AT‑AGC‑AAA‑CCCC‑CCA‑ccctt                                                                                                                                                                                                                        >  2:240796/1‑235 (MQ=255)                      gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGAAAT‑TCCCCCCAACATGGCGGG‑CT‑AAGG‑TTGCGGGA‑AAACCA‑AAC‑CCCG‑CCC‑C‑TTTTA‑ccc                                                                                                                                                                                                                  >  1:452369/1‑239 (MQ=255)                           gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCCAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGAAAATTCCAGCGGTTCGGTAAT‑TCCCACCAGCATGGC‑GGTTT‑AAGG‑TTCGGGGA‑TACGCC‑AAA‑CCCC‑CCC‑A‑TTTTA‑CCCGAtt                                                                                                                                                                                                              >  1:1064030/1‑238 (MQ=255)                           gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGGGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAAT‑TCCCACCAGCATGGC‑GG‑GT‑AAGG‑GTTCGGGG‑ATCAGC‑AAA‑CCCC‑CCC‑CTTTTTA‑CCCGCttt                                                                                                                                                                                                             >  1:158819/1‑239 (MQ=255)                                    cATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGATCGGTAAT‑TCCCACCAGCATGGC‑GG‑GT‑AAGG‑GTGCGGGA‑ATCAGC‑AAG‑CCCG‑CCC‑C‑TTTTA‑CCCGATTTTCCGGTGCgg                                                                                                                                                                                                   >  1:874162/1‑239 (MQ=255)                                    cATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGGGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAAT‑TCCCACCAGCATGGC‑GG‑TT‑TAGG‑TTGCGGGA‑TACCAC‑AAGCCCCC‑CCC‑C‑TTTTT‑CCGCATTTTctgcgcgg                                                                                                                                                                                                    >  1:363822/1‑232 (MQ=255)                                                    tGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAAT‑TCCCACCAGCATGGC‑GG‑TT‑AAGG‑TTGCGGGA‑ATCAGC‑AAA‑CCCG‑CCC‑C‑ACTTT‑CCCGCCTTTCCCGTGCGGCAGGGaaagaaaagcc                                                                                                                                                                                   >  1:344111/1‑230 (MQ=255)                                                                 aGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAAT‑TCCCACCAGCATGGC‑GG‑TT‑TAAG‑GTGCGGGG‑ATCAAC‑AAA‑ACCG‑CCC‑C‑CCTTT‑TCCGGTTTTTCGGTTGGGGGGGGGAGGACATCGCTAaaaaaaaaccc                                                                                                                                                                      >  1:431863/1‑234 (MQ=255)                                                                   tCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAAT‑TCCCACCAGCATGGC‑GG‑TT‑AAGG‑GTTCGGGA‑ATCAGCAAAA‑CCCC‑CCC‑C‑TTTTA‑‑CCGATTTTTCGGGTCGGCAGTGAAGGACAATGCTAAAGAAAAGcccc                                                                                                                                                                     >  1:768177/1‑237 (MQ=255)                                                                    cAAACAGCCCCCCCCCCATGTTCCCCACAAACCAAAAGGAATAAATTAGGGCAGCAATGGGGGCCGCCAGCGCGGGGTGTACCGAAAACAGCAACGGTAAAAGGAACAGGAAGGTAAATCCCGGCGGTTCGGTAAT‑TCCCACCAGCATGGC‑GT‑TT‑AAGG‑TTGCGGGA‑ATCAGC‑AAC‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCa                                                                                                                                                                     <  1:71096/237‑1 (MQ=255)                                                                               ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAAT‑TCCCACCAGCATGGCGGG‑TT‑AAGG‑‑TGCGGGG‑ATCAGC‑AAA‑CCCC‑CCACC‑TTTTA‑CCC‑ATTTTTCGGTGCGGCCGGGAAAGAAATCGCCCAAAAAAAACACACAGCA‑CCcaaaaa                                                                                                                                                       >  2:136068/1‑236 (MQ=255)                                                                                      gtttccccACCACACCAAAGAGATACATGACGGTCGACATGGAGGCCGCCAGCACGGCTTGAACCGAAACCAGCAAGGGTGAAATGAACAGGAAGGTAAATTCCGGCGGTTCGTTAAT‑CCCCCCCAGCATGGC‑GT‑TT‑AAGG‑TGGCGGAA‑AGCAGC‑AAG‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAg                                                                                                                                                <  1:265633/234‑1 (MQ=255)                                                                                       gTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAAT‑TCCCACCAGCATGGC‑GG‑TT‑AAAG‑GTGGGGGA‑ATACGC‑AAG‑GCCC‑CCC‑C‑TTTTATCCCGA‑TTTCCGGTGCGGCGAGGCAGTGCATCGTCAAAGAAAATGCCACGGCC‑CCCAAAAGTTTTgtgt                                                                                                                                               >  1:543156/1‑239 (MQ=255)                                                                                              ccccaccaaaaGATATATATTGCGGTGGAGAGTGGGCCCGCCAGCCCGGGGTTTACAGAAAACAAAAAGGGGAAAGGGAAGGGGAAGGAAATGCCCGGGGGTTGGGTAAT‑CCCCACCAGCAGGGG‑GG‑TA‑AAGG‑TTGGGGAA‑ACCAGA‑AAG‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTg                                                                                                                                        <  2:344111/231‑1 (MQ=255)                                                                                                             tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAAT‑TCCCACCAGCATGGC‑GG‑TT‑AAAG‑ATGCGGGG‑ATCAGC‑AAA‑CCCC‑CCC‑C‑TTTTT‑CCCCCTTTTTCGGGTGGGGAGTGTAGTAAATTGCTCAAAAAATACCCCCGGCACCCCAAAATTTTTG‑GGTACCCTGCAGGGCAAAACAGc                                                                                                                          >  1:274777/1‑238 (MQ=255)                                                                                                                  tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAAT‑TCCCACCAGCATGGCGGG‑TT‑AAGG‑TTGGGGGA‑A‑CAGC‑CAA‑CCCGCCCA‑C‑TTTTT‑CCCGCTTTTCCGGGGGGCGGGGGAAGACAACCGCCAAAAAAATACCCCCCGC‑‑CCAAAAACTTTTTAGTTATCCCGCAGGGGaaaaaacccccccg                                                                                                                    >  2:548329/1‑229 (MQ=255)                                                                                                                                             gggcttccccccAAACAGCAAGGGTAAAAGGAAGGGGAAGGACAATCCCGGCGGTCGGGTACT‑CCCCCCCAGGAGGGT‑GT‑TT‑AGGG‑TTGCGGAA‑AGCGGC‑AAC‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGgcgc                                                                                         <  2:549155/228‑1 (MQ=255)                                                                                                                                                     ccGCAAACAGCAAGGGTGAAATGAACGGGAAGTTAAATTCCAGCGTTTCGGTACT‑CCCCCCCAGCATGGC‑GT‑TT‑AAGG‑TGGCGGGA‑ATAAGC‑AAG‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTg                                                                                 <  1:1065356/239‑1 (MQ=255)                                                                                                                                                                     gtaaagaacgaGGATGGT‑AATTCCCGTGGTTTGGTTATCTCCCACCACTTTGGG‑GG‑TT‑GAGT‑GCGGGGCA‑AACAGC‑AGC‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgc                                                                  <  2:497382/228‑1 (MQ=255)                                                                                                                                                                     aaaaaTGAACAGGAAGGTAAATTCCGGGGGTTCGGTAAT‑CCCCACCGGCAGGGC‑GG‑TA‑AAGG‑TGGCGGGA‑ACCAGC‑AAC‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                 <  1:104817/237‑1 (MQ=255)                                                                                                                                                                        agtGAAGAGGAAGGAAAATTCCAGCGGTTCGGTAAC‑TCCCCCCAGCGTGGG‑GG‑TT‑AAGG‑TTGCGGAA‑ATCAGC‑AAG‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                              <  1:136068/237‑1 (MQ=255)                                                                                                                                                                        aaTGAAAAGGAGGGTCAAGTCTGGTGGTGCGGAAAT‑CCCCCCCGGCGTGG‑‑GG‑TT‑AAGG‑GTGGGGGAGACAAGC‑AAC‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                              <  2:263413/239‑1 (MQ=255)                                                                                                                                                                        aaTAAACGGGAGGGAAAATCCCAGCGGTCGGGAAAT‑TCCCACCAGCAGGGC‑GG‑TT‑AAGG‑TGGCGGAA‑ATCAGC‑AAG‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                              <  1:334090/239‑1 (MQ=255)                                                                                                                                                                            aaggggaaggaaAATCGCGGGGGTTGGGTAAC‑CCCCACCAGGGGGGG‑GG‑TG‑AAGG‑TGGGGGAA‑AACAAC‑ACC‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                            <  2:1024350/233‑1 (MQ=255)                                                                                                                                                                             aaagagaGGTAAATTCGCGGGGTTCGATAAC‑TCCCTCCAGCAGGGC‑TG‑TT‑AAGG‑TTGCGGAA‑ATCAGC‑AAG‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                         <  1:182702/233‑1 (MQ=255)                                                                                                                                                                             aCAGGAGGGAAAATCCCGGCGGTTCGGAAAT‑TCCC‑CCAGGATGGCGGG‑TT‑AGGG‑TTGGGGAA‑ACCAGC‑AAG‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                         <  1:14295/239‑1 (MQ=255)                                                                                                                                                                             aCAGGAGGATAAATCCCAGCGGTCCGAAAAT‑CCCCCCCAGCTGGGC‑GG‑TT‑AAGG‑TGGGGGAA‑ATCAGC‑AAC‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                         <  1:70075/239‑1 (MQ=255)                                                                                                                                                                                         aTTCCGGGGTTTGGGAAAT‑CCCCACCCGTAGGGC‑GG‑TA‑AGGG‑TTGCGGGA‑CTCAGA‑ACC‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                             <  1:261763/239‑1 (MQ=255)                                                                                                                                                                                              gcggccgcGGTAAT‑CCCCCCTGGCGTGGC‑GG‑TT‑AAGG‑TTGGGGAA‑AACAGC‑AAC‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCGAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGAGTGCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgccg                                          <  1:821959/230‑1 (MQ=255)                                                                                                                                                                                                      cgttattccccACCGGCTTGGG‑TG‑TT‑ATGG‑GGGCGGGA‑AACAGC‑AAC‑CCCG‑CCA‑C‑GTTTA‑CCCTATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACa                                 <  1:1010277/231‑1 (MQ=255)                                                                                                                                                                                                                               ttt‑AAGG‑TTGGGGAA‑ACCACC‑AAG‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACt           <  1:388030/237‑1 (MQ=255)                                                                                                                                                                                                                                          gtgtggcaTCAGC‑AAC‑CCCG‑CCC‑A‑TTTTT‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAATATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGa  <  1:111014/232‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                              ATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAAT‑TCCCACCAGCATGGC‑GG‑TT‑AAGG‑TTGCGGGA‑ATCAGC‑AAG‑CCCG‑CCA‑C‑TTTTA‑CCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGA  >  NC_000913/3862435‑3862900

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