Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	990,365	T→G	22.6%	S138S (TCA→TCC)	pncB ←	nicotinate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	990,365	0	T	G	22.6%	150.9 / 2.2	77	S138S (TCA→TCC)	pncB	nicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (18/36);  new base G (13/8);  total (32/45)
Fisher's exact test for biased strand distribution p-value = 3.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATC  >  NC_007779/990129‑990599                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           gtTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTg                                                                                                                                                                                                                                           <  1:1059553/238‑1 (MQ=255) gtTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTg                                                                                                                                                                                                                                           <  1:1387386/238‑1 (MQ=255)   tCCCCTCGGCGTGAAGGAAAGCGGGCGCGCCGGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTGGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGTCCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGag                                                                                                                                                                                                                                         <  2:1227458/238‑1 (MQ=255)   tCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGag                                                                                                                                                                                                                                         <  1:948646/238‑1 (MQ=255)       aTCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGGTAACGCCGCGAAAGCGACTAATTTTT‑TTTTCAG‑CGTGGCGGGGGCGTGC‑GC‑CACGT‑CGGCCTGCGGGGagggata                                                                                                                                                                                                                                    >  1:514735/1‑239 (MQ=255)        tCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGagcgatag                                                                                                                                                                                                                                   <  1:967367/239‑1 (MQ=255)          ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCTAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                 <  1:1867622/239‑1 (MQ=255)          ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                 <  2:733671/239‑1 (MQ=255)          ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                 <  1:156272/239‑1 (MQ=255)          ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCGGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                 <  1:126229/239‑1 (MQ=255)          ggCGTGAAGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                 <  1:1120748/239‑1 (MQ=255)           gCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGGCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTTTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTg                                                                                                                                                                                                                                           <  1:1916812/228‑1 (MQ=255)           gCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TGTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                 <  2:1232274/238‑1 (MQ=255)            cGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAAAGCCGCGAAGGAGACTCATTTGC‑TTTTCAGACGGGT‑GGGGGGGTGC‑GC‑CACGG‑CGGGCTGCGGGGggcggtagagat                                                                                                                                                                                                                               >  1:723279/1‑239 (MQ=255)             gTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACACACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATg                                                                                                                                                                                                                              <  1:966238/239‑1 (MQ=255)             gTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATg                                                                                                                                                                                                                              <  1:1048705/239‑1 (MQ=255)             gTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATg                                                                                                                                                                                                                              <  1:1392710/239‑1 (MQ=255)             gTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAAGCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATg                                                                                                                                                                                                                              <  1:127355/239‑1 (MQ=255)               gAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTa                                                                                                                                                                                                                            <  1:168327/239‑1 (MQ=255)                            cgcgCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑CACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGAt                                                                                                                                                                                                               >  1:1992816/1‑239 (MQ=255)                                 cAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTa                                                                                                                                                                                                                            <  1:1743788/221‑1 (MQ=255)                                  aGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGGTAACGCCCAAAATATGCATAATTTTC‑CTTTCAG‑CGGGTCGAGGGCCGTC‑GC‑AACCT‑CGGCCTTGGGGTAGCGTGAGCGATGTACCATTTCACTGACaaact                                                                                                                                                                                                          >  2:195094/1‑235 (MQ=255)                                                       acaAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGAGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAAAGCACAGAAGTCGACTAAATTGC‑TTTCCAG‑CGTGTCGAGGGCGTGC‑GC‑AACGT‑CGGCCTGCGGGGGGCGATAGCGATGTACCACTTCTCTGATAAACGCCCGCAAAGGAACCTCCCagag                                                                                                                                                                                    >  1:160297/1‑239 (MQ=255)                                                         aaaCCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGGTAAAGCCGAGAAGTGGACTAATTTGG‑TTTTCAG‑CGTGGCGAGGGGTTGC‑GCGAACGT‑‑GGCCCGCGGTGAGGGAGAGAGCAGTACCACTTCTCTGAAAAACCCCAGCCAAGGAACCTTCCCGAGGa                                                                                                                                                                                  >  2:1387799/1‑239 (MQ=255)                                                                 aTTCCTGTGGCGGCCGCTAAACGGTGGTTTCTTGGCCTCCGGGGAAAAACCGGCGCGGGGGGCAAAAACCCGTCGGATGGAGGCGCAAAAAATAA‑AGCCCGGCGGTTACCGCCGAGAAGGCGATTAATTGGC‑TTTCCGG‑GGGGTGGGGGGCTTGC‑GC‑AGCGT‑CGGCCTGGGGTGGGCGATAGAGGTGTCCCTTTCCACGGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTc                                                                                                                                                                           <  2:1222456/238‑1 (MQ=255)                                                                     cTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGGCGACTAATTTTC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGGGAGCGATAGCGATGTACCATTTCACTGATAAACGCCCGCAAAAGaa                                                                                                                                                                                              >  1:1317638/1‑215 (MQ=255)                                                                       gTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCCAGAAGTCGACTAATTTGC‑TTTTCAG‑CGGGTCGAGGGCTTTC‑GC‑AACGT‑CGGCCTGCGGGGAGCGATAGCGATGTACCCTTTTACTGATAACCGCCAGCCAAAGAACTTCCCCGAGGATGACTTTCCGccc                                                                                                                                                                     >  1:1353043/1‑237 (MQ=255)                                                                                gCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGGGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCtt                                                                                                                                                           >  1:88441/1‑239 (MQ=255)                                                                                    aaCGATGGTTTCTTGTGCTCCGGGA‑AAAAACGGGGGCGGGGGCCAAAAACCAGCGGAAGGGAGGGGCGCATGTCAGACACGGGCGTTTACCGCGGAGAGGGCGGCTAATTGGT‑TTCCCGG‑TGTG‑GGGGGGGTTGCGGC‑AGCGT‑CGGCCGGCGGGGAGGGAGAGCGTGCTCCGATTTCCCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAAc                                                                                                                                                       <  2:277131/239‑1 (MQ=255)                                                                                         tGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGGCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGG‑CGGCCTGCGGGGAGCGATAGCGATGTACCCTTTCACTGATAACCGCCAGCAAAGGaa                                                                                                                                                                                              >  1:1764203/1‑195 (MQ=255)                                                                                           gTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑GGACCGGGGGTTAACGCCGAGAAGTCGCCTAATTTGC‑TTTCCAG‑CGTGTGGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCTTTTCACTGATAACCGCCAGCAAAgg                                                                                                                                                                                                <  1:385115/191‑1 (MQ=255)                                                                                            tttCTTGTACTTCGCGAGAAAAACGGCGACGGGTGTCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAAAGCCGAGAAAGCGACCAATTTGA‑TTTTCAG‑CGTGTCGGGGGGTTGC‑GC‑AACCG‑CGGCCTGCGGGGGGGGAGAGCGATGTACCATTTTACTGATAACCGCCCGCAAAAGGACATTCCCGAGGAGGGAATCTCCCCACGGGCCGCCTTAACAAatcgcc                                                                                                                                               >  2:543117/1‑235 (MQ=255)                                                                                              tCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGGTAACGCCCAGAAGTCGACTCATTTGT‑TTTCCAG‑CGTGTCGAGGGGATGC‑CC‑AACAT‑CGCCCCCCGGCGGGCGATCGCGATGGATCGTTTCTCTTATAACCGCCCGCCAACGAAGTTCCCAGAGGGTGAATTCACACCCCTCACCGCCTTAACCAATATcac                                                                                                                                              >  2:1777142/1‑236 (MQ=255)                                                                                                     cTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCTTCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTTACGCCCAGCAATCGACTAAATTTC‑TTTCCAC‑CGTGTCGAGGGGCTGC‑GC‑CACGT‑CGGCCTGCGGGGAGCGCTAGCGATGTACCATTTCTCTCATAACCGCCAGCAAAGGGACATCCCCGAGGGTGACTTTACGCCCCAGGGCCCGTAAACGAAAATcccacctgcc                                                                                                                                       >  2:554552/1‑230 (MQ=255)                                                                                                      ttCGCGAGAAAACCGGCGACGGGTGCAAAAATCCTTCAGAGGGAAGCGCGAAATACCA‑AGCCCGTCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTCCCAG‑CGTGTCGAGGGTTTGC‑GA‑AACGT‑CGGCCTGCGGGGAGAGATAGCGATGTACCTTTTCACTAATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCAt                                                                                                                                     <  2:1353043/239‑1 (MQ=255)                                                                                                        cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCCAGAAGTCGAATAATTTTC‑TTTTCCG‑CGCGGCGAGGGGCTGC‑GC‑CACGT‑CGGCCTGCGGGGGGGGGTAGCGGTGGTACATTTCACTGGTAACCGCCCGCAAAGGGACTTCCCCGAGGATGACTTCACCCCACGGGCCGCCTAAACAAAAATCCACCCTGCCATTa                                                                                                                                   >  2:156904/1‑239 (MQ=255)                                                                                                        cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCCAGAAATCGAATAATTTTC‑TTTCCCA‑CCTGTCGAGGGGGTGC‑CC‑AACGT‑CGGCCTGCCGGGAGCGATAGCGATGTAACATTTCACTGATAACCGCCCGCAAAGGAACTCCCCCGAGGGTGACTTCACGCCCCGGGCCCCCTACACGAAAATACAACCTGcgcttt                                                                                                                                   >  2:731884/1‑234 (MQ=255)                                                                                                             gaAAAACGGCGAGGGGTGCCAAAGCCCATCAAAGGGAAGCGCTGAACAACA‑AGCCGGGGGTTTCTCGCCGGGAAGCCGCATATTTTGT‑TTGCCGG‑AGGGAGGGGGGTTTGC‑‑C‑ACTTTCCGGCCGGGGGGGGGCAATGCGGGTGAACCTTTTCATTGAAATCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCg                                                                                                                                <  2:1539586/237‑1 (MQ=255)                                                                                                                         cGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGCCCGGCGGTTACCGCCGAGAAGTCGACTATTTTGC‑TTCCCAG‑CGTGTCGAGGGCTTGC‑GC‑ACCGT‑CGGCCGGCGGTGAGCGATAGCGATGTCCCATTTCACTGATACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAc                                                                                                                  <  1:156904/239‑1 (MQ=255)                                                                                                                               ccAAAATCCATGAGGTGGAACGGCGACATATCA‑GGCCCGGCGGTTCAGGCCGAGAAGTCGTCTAATTTGC‑TTTCGCG‑CGTGTCGAGGCCTGGC‑GC‑AACGG‑CGGCCTGGGGGGAGCGATAGAGATGTACGATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTc                                                                                                            <  2:1850461/239‑1 (MQ=255)                                                                                                                                     tCCCTCAGAGGGAAGCGCGAAATATCA‑AGACCGGCGGTTAACGCGGAGAAGTCGCCTAATTTGT‑TTTCCAG‑CGTGTGGAGGGCTTGC‑GC‑AACGT‑CGGCCTGGGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGtt                                                                                                      <  1:731884/239‑1 (MQ=255)                                                                                                                                                  ggcgcGAATTATCA‑AGCCCGGCGGTTAACGCGGAGAGGTCGACTAATTTGT‑TTTCCAG‑CGTGTCGAGGGTTGGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCTTTTCACGGATACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTc                                                                                                                                                                           <  2:155283/156‑1 (MQ=255)                                                                                                                                                       gataaacga‑aGACCGGCGTTTACGGCGGAAGAGCCCATTATTTTGT‑TTCCCGG‑CGTGTCGAGGGTTGGC‑GA‑AACGT‑CGGCCGGCGGTAGGCTATGGCGATGACCCATTTCACAGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                    <  2:723279/231‑1 (MQ=255)                                                                                                                                                       gaaatatCA‑AGCCCGGGGGTGAACGCCAAGAGGTCGACTATTTTGT‑TTCGCGG‑CGTGTGGAGGGCTTGC‑GC‑AACGT‑CGGCCGGCGGGGAGCGATAGCGATGTCCCATTCCACTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                    <  2:1587940/236‑1 (MQ=255)                                                                                                                                                       gACATACCA‑GGACGGGCGGTTACGGCGGAGAAGCCGATTAATTTGC‑TTTCCAG‑GGTGTCGAGGGTTTGC‑GC‑AACGT‑CGGCCTGCGGTGGGGGATGGCGATGTACCATTTCCCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                    <  1:543117/239‑1 (MQ=255)                                                                                                                                                         ataatCG‑AGACCGGCGGGTAACGGCGAGAGTGCGACTATTTGGT‑TTTCCAG‑CGTGTGGAGTTCTTGC‑GA‑GACGT‑CGGCCTGGGGTGAGAGATAGCGATGCACCATTTCTCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAAcc                                                                                  <  2:160297/236‑1 (MQ=255)                                                                                                                                                            atCA‑AGACCGGCGGTTAACGCCGAGAAGTCGGCTAATTTTC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑CACGT‑CGGCCTGCGGGGAGCGATAGCGATGTAACCTTTCACTGATAACCCCCAGCCAAAGAACTTCCCCGAGGGTGACTTTACGcccccgg                                                                                                                                                                 >  1:155283/1‑152 (MQ=255)                                                                                                                                                            atCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGTTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTa                                                                                                                                   <  1:1869731/187‑1 (MQ=255)                                                                                                                                                                      ggcggTAAACGCCGAGAAGTCGACTATTGTGC‑TTTCCGG‑GGTGTCGCGGGCTTGC‑GC‑AACGT‑GGCCCGGCGGGGAGCGATAGCGATGTACCATTTCACTGAAACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                      <  1:1656320/239‑1 (MQ=255)                                                                                                                                                                       ggggtaaccGCCGAGAGGGCGCTTATTTGGC‑TTCCCAG‑CGTGTCGAGGGCTTGC‑GA‑AACGT‑CGGCCTGGGGTGAGCGATAGCGATGTCCCATTTCACTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                      <  1:506449/230‑1 (MQ=255)                                                                                                                                                                       gcgtTTAGCGCCGAGATGT‑GACCTATTTGCGTTTCCAG‑GGTGTCGGGGGCTTGC‑GC‑AACGT‑CGGCATGGGGGGGGCGGTAGCGATGTACCATTTCCCTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                      <  2:1240822/238‑1 (MQ=255)                                                                                                                                                                       gcggTTAACGCGGAGAAGTCGACTATTTTGC‑TTCCCGG‑GGTGTCGAGGGTTTGC‑GC‑AACGT‑GGGCCTGGGGTGAGCGATAGCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                      <  1:1759609/238‑1 (MQ=255)                                                                                                                                                                            gaccGCCGAGATCTGGACTAATTGGC‑TTCCCAG‑CGTGTCGGGGGCTTGC‑CC‑AACGT‑GGGCCTGCGGAGAGCAATAGCGATGCACCTTTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCt                                                                <  1:1732736/236‑1 (MQ=255)                                                                                                                                                                            gaaCGCCGAAAAGTGGACTATTTTTC‑TTTCCAG‑CGTGTGGAGGGTTTGC‑GC‑AACGT‑CGGCCGGCGGTGAGCGATAGCGGTGTACCATTTCCCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCt                                                                <  1:304979/238‑1 (MQ=255)                                                                                                                                                                                  cGAGATGTCGTGTATTTTGC‑TCCCCAG‑CGTGTGGCGGGCTTGC‑GA‑AACGC‑CGGCCGGGGGTGAGCATTAGCGATGTACCTTTTCATTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGGACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaag                                                           <  2:514735/238‑1 (MQ=255)                                                                                                                                                                                    agaAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaagaaa                                                        >  1:37923/1‑239 (MQ=255)                                                                                                                                                                                         tCGATTAATTTGT‑TTTCCGC‑GGTGTCGAGGGTTGGC‑GA‑AACCT‑CGGCCTGCGGTGAGGAATTGCGATGTCCCATTTCCCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggca                                                    <  1:1387799/238‑1 (MQ=255)                                                                                                                                                                                         tCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggcag                                                   >  1:302642/1‑239 (MQ=255)                                                                                                                                                                                          cGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCGTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAgccagg                                                  >  2:2055140/1‑239 (MQ=255)                                                                                                                                                                                                   tGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAg                                          >  1:62272/1‑238 (MQ=255)                                                                                                                                                                                                     tccttcccggggtgtGGAGGGTTTGC‑GA‑ACGGC‑CGGCCGGGGGTGAGGGATAGGGAGTAACCATTCAACAGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGc                                         <  2:1351703/228‑1 (MQ=255)                                                                                                                                                                                                     c‑TTTCCAG‑CGTTTGGAGGGCTTGC‑GC‑AACGG‑CTGCCTGGGGTGGGCGATAGCGATGTCCCTTTTCCGTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCa                                       <  1:1620775/239‑1 (MQ=255)                                                                                                                                                                                                     c‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGcc                                        >  1:219619/1‑238 (MQ=255)                                                                                                                                                                                                       gtccgaG‑CGTGTCGAGGGCTTGC‑GC‑ACGGA‑GGGCCTGGGGGGGGGGATAGCAATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCAc                                      <  1:1424197/234‑1 (MQ=255)                                                                                                                                                                                                         tCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAAc                                                                                   >  1:1296299/1‑192 (MQ=255)                                                                                                                                                                                                                 tgtCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATtca                             >  1:1850281/1‑239 (MQ=255)                                                                                                                                                                                                                   tCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATtcatc                           >  1:1901101/1‑239 (MQ=255)                                                                                                                                                                                                                      aGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTc                                                                                                                                                                                                                     <  2:17711/50‑1 (MQ=255)                                                                                                                                                                                                                      aGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTc                                                                                                                                                                                                                     >  1:17711/1‑50 (MQ=255)                                                                                                                                                                                                                           ttGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCa                   >  1:159560/1‑239 (MQ=255)                                                                                                                                                                                                                             gc‑gc‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTaa                                                                                                                                                         >  1:1831811/1‑103 (MQ=255)                                                                                                                                                                                                                             gc‑gc‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGACCTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                  >  2:597454/1‑238 (MQ=255)                                                                                                                                                                                                                                   caactccGGCCTGCGGTGGGCGTTACCGATGCCCCTTTTCACTAATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGgcgc              <  1:1324673/233‑1 (MQ=255)                                                                                                                                                                                                                                     cGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAgg           >  2:399894/1‑239 (MQ=255)                                                                                                                                                                                                                                              tGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa    >  2:979872/1‑238 (MQ=255)                                                                                                                                                                                                                                              tGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCAAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTAATCATCCTGCAGGCGCAGGAGCTGCa    >  1:980178/1‑238 (MQ=255)                                                                                                                                                                                                                                                cGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATc  >  2:2015006/1‑238 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           GTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCA‑AGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAG‑CGTGTCGAGGGCTTGC‑GC‑AACGT‑CGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATC  >  NC_007779/990129‑990599

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