Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	23,956	T→C	D522D (GAT→GAC)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	23,956	0	T	C	91.8%	156.9 / 3.0	64	D522D (GAT→GAC)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (22/34);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

GCGTGCAGTGGATCCCGGACTGGGGCCAGGCGCGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGTGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGATCA‑GTACGTGAA‑AGTGCCGGACACAT‑TGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAGGGTAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATG  >  NC_000913/23709‑24199                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       gCGTGCAGTGTGTCCCGGACTGGGTCCAGGTGCGTAGTGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTTGCTCCCGTCAGCGCTCCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGTAAGAAGTGGCAAAACGCGTGGAAGTCAACGGCATCCAGGTGTGGTGGGATCTCGATTCGAAAGATATCCTCGGC‑GACGAAGCTGACCA‑g                                                                                                                                                                                                                                                    <  1:330121/251‑1 (MQ=255)      cAGTGGATCCCGGACTGGGGCCAGGCGCGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTatgt                                                                                                                                                                                                                                               <  1:59856/251‑4 (MQ=255)               ccGGACTGGGGCCAGGCGCGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGcc                                                                                                                                                                                                                                     <  1:221693/251‑1 (MQ=255)                cGGACTGGGGCCAGGGGCGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGGACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCg                                                                                                                                                                                                                                    <  1:278266/251‑1 (MQ=255)                cGGACTGGGGCCAGGCGCGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTTCCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTAGATGCGAAAGAGATCCTCGGC‑GACGAAGTTGACCA‑GTATGTGAA‑AGTGCCg                                                                                                                                                                                                                                    <  1:137642/251‑1 (MQ=255)                  gACTGGGGCCAGGCGCGTATCGAGTTGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGa                                                                                                                                                                                                                                  <  1:13802/251‑1 (MQ=255)                  gACTGGGGCCAGGCGCGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGa                                                                                                                                                                                                                                  <  1:209828/251‑1 (MQ=255)                              gcgcGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGatactctcgacgt                                                                                                                                                                                                                      <  1:288056/251‑13 (MQ=255)                                  gTATCGTGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTTGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGatactctcgacgtatgg                                                                                                                                                                                                                  <  1:228158/251‑17 (MQ=255)                                   tctCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACCAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGCTGGAAGAAGTGGCAAAACGCGTTGAAGTTACCGGCATCCAGGCGTGGGGGGATCGCGTTGCGATAGAGATCCTCGGCAGA‑GAAGCTGACCA‑GTATGTGAA‑AGGGCCGGatactctcgaggtagggt                                                                                                                                                                                                                 >  2:282445/3‑234 (MQ=255)                                        aGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGAc                                                                                                                                                                                                           >  1:343441/1‑251 (MQ=255)                                           cggtGGTTTCTAACCGTCCTGACTGGTGTCACTCCCTGCAGCGCACCTGGGGCGTACCGATGTCACTGGTCGGTCACAAAGACGCGGAAGAGCTGCATCCGCGTCCCCTCGAACTGATGGAATAAGTTGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACTAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTcc                                                                                                                                                                                                        <  2:142012/248‑1 (MQ=255)                                            gctggtttatatccGTCCTGACTTGTGTATCTCCCGTCAGCGCACCTGGGGTGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGACAGAGATCCTCGGC‑GACGAAGCTGATCA‑GTACGTGAA‑AGTGCCGGACACAT‑TGGCAGTATGGTTTGAATCCg                                                                                                                                                                                                       >  2:251117/13‑251 (MQ=255)                                                      aaCCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTAcc                                                                                                                                                                                               <  1:43918/249‑1 (MQ=255)                                                               gACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTTCCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCtg                                                                                                                                                                                     <  1:97079/250‑1 (MQ=255)                                                                 cTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCC‑GTATGTGAA‑AGTGCCGGATACTC‑TCGCAGTATGGGTTGACTCCGGATCTACCCACTCTTCtgttg                                                                                                                                                                                  >  2:66719/1‑251 (MQ=255)                                                                              cGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTcc                                                                                                                                                                     <  1:130454/251‑1 (MQ=255)                                                                                    cgcACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAAtt                                                                                                                                                               <  1:9669/251‑1 (MQ=255)                                                                                        ccTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGcc                                                                                                                                                           >  2:88762/1‑251 (MQ=255)                                                                                          tGGGGCGTACCGCTGTCACTGTTCGTGCCCAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGCACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCACGGC‑GACGAAGCTGA‑CACGTATGTGTA‑AGAGCCGGATA‑ATCTCGACGTATGGGTTGACTACGGATCTACCCACTCTTCTGTTGTTGCCTTTCGTCCGGAATTTGCCgg                                                                                                                                                         >  2:200705/1‑251 (MQ=255)                                                                                           gtggcgtACCGATGTCACTGTTCGGGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGGGTGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGGATGGTTTGACTCCGGATCTACCCACTCTGCTGTTGTTGACGTGCGTCCGGAATTTGCCg                                                                                                                                                          <  1:151906/244‑1 (MQ=255)                                                                                             ggcgttccattttctcTGTACGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCATCGAACTGATGGAAGAATTGGCAAAACGAGTTGAAGTTAACGGCTTCCAGGTGTGGTGGGTTCTCGTTGCGAAATAGATCCGCGGA‑GACGAAGCTGACAA‑TTATGTGAA‑AGTGCCGGATGCTC‑TCGGCGTATGGTTTGACTTCGGATGTACTCACTCTTCTGTTTTTGTCGTGCGTCCGGAATTTGTCGGtta                                                                                                                                                      <  2:287993/236‑3 (MQ=255)                                                                                                gtACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGAGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAACAGT‑CCGGATACAC‑ACCAAGTATGGTTTGACTCCGGATCTACCACCTCTTATGTTGTTGACGGGCGTACGGCATTTGCCGGTCACGc                                                                                                                                                   >  2:338675/1‑251 (MQ=255)                                                                                                gtACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAAATGATGGCAGAAGTGGCAAAACGCGTTGCAGTTAACGGCATACAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑CGTGCCGGAGACTC‑TCGGCGAATGGTTTGACTCCGGATCTACCCACTCTTCTTTTGTTGACGTGCATCCGGAATTTGCCGGTctctc                                                                                                                                                   >  2:322366/1‑247 (MQ=255)                                                                                                 tACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGc                                                                                                                                                   >  1:369161/1‑250 (MQ=255)                                                                                                    cGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCg                                                                                                                                               >  1:60791/1‑251 (MQ=255)                                                                                                          cACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATg                                                                                                                                         >  1:189274/1‑251 (MQ=255)                                                                                                              gTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATc                                                                                                                                     <  1:25955/251‑1 (MQ=255)                                                                                                                 cGTGCACAAAGACACGGAAGAGCTGCATCCGCGTAACCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCC‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTgg                                                                                                                                  >  1:344671/1‑251 (MQ=255)                                                                                                                   ttcacaAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGACGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTTGTGGGATATCGATGCGAAAGAGATCCTCGGC‑GACGATGCTGACCA‑TAATGTGAA‑AGTGACGGATACTC‑TCGCCGTATGGTTTGACTCCGGATATACCCACTCTGCTGCTGTTGACGTGCGTCCGGAGTTTGCAGGTAACGCAGCGGAccagtatccggt                                                                                                                                 >  2:265887/3‑239 (MQ=255)                                                                                                                       caAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGGTGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGTCCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTGGACTCCGGATCTACCCACTCTTCTGTTGTTGACGGGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTACCTGGAAGGtt                                                                                                                            <  1:208257/251‑1 (MQ=255)                                                                                                                             cacGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGAcc                                                                                                                      <  1:94342/251‑1 (MQ=255)                                                                                                                                     agctgcatccgtgggacctcggagtGATGGAAGAATTGGCAAAATGCGCTGAAGGTAACGGCATCCAGGCTGGGTGTGAGCTCGCTGCGAAAGAGTTCCTCGTC‑GACGAAGCAGACCA‑GTATGTGAA‑AGTGCCGGAGACTC‑TCGACGTATTGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGAAGTGCGTCCGGAATTTGCCGTTCACGCAGCGGACATGTATCTGGAAGGTCCTGACCAGCACcgt                                                                                                              <  2:305355/227‑2 (MQ=255)                                                                                                                                     agCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTAACGCAGCGGACATGTATCTGGAAGGTTCTGAcccgcaccgt                                                                                                              >  1:221614/1‑243 (MQ=255)                                                                                                                                      gCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACcgtg                                                                                                             <  1:312788/251‑3 (MQ=255)                                                                                                                                        tGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACcgtggt                                                                                                           >  1:170413/1‑247 (MQ=255)                                                                                                                                               gcgTACCATTGAAATGATGGAAGAAGTGGAAAAAAGAGTTGAAGTAGATGGTATCAAGGAGTGGTGGGATATAGATGAGAAAGAGATAATAGGA‑GAAGAAGATGATCA‑GTAAGTGAA‑AGTGCAGGAAAAAT‑TGGATGTATGGTTTGAATACGGATATAAAAAATATTATGTTGTTGACGTGCGTACGGAATTTGCAGGTAACGCAGCGGACATGTATCTGGAAGGTTCTGACCAACACCGAGGCTGGTTCa                                                                                                    <  1:251117/251‑1 (MQ=255)                                                                                                                                                cgaaacctcgAACAGATGGAAGAAGTGGGAAAACGCTGTGAAGTTAACGGCATCCAGGCGTGGCGGGATCTCGATGCGAAATAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGGATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGGCCGGAATTTGCTGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACcgtggtttgttcga                                                                                                   <  2:13978/242‑13 (MQ=255)                                                                                                                                                         gAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCt                                                                                          <  1:99796/251‑1 (MQ=255)                                                                                                                                                             tGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GAAGAAGCTGATCA‑GTACGTGAA‑AGTGCCGGACACAT‑TGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATg                                                                                      >  1:229218/1‑251 (MQ=255)                                                                                                                                                              gATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGATCA‑GTACGTGAA‑AGTGCCGGACACAT‑TGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGTCATGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATGa                                                                                     <  1:16331/251‑1 (MQ=255)                                                                                                                                                                 ggAAGAAGTGGCAAAACGCGTTGAATTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGATGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGAtct                                                                                  <  1:242797/251‑1 (MQ=255)                                                                                                                                                                     gaagTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTTATCA‑GTACGTGAA‑AGTGCCGGACACAT‑TGGATGTATGGTTTGACTACGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGGAGGTTCTGACCAACAACGCGGCTGGTTCATGTCTTCCCTAATGATCTCCAc                                                                              <  1:156967/251‑1 (MQ=255)                                                                                                                                                                                aCGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGaa                                                                     >  1:187875/1‑249 (MQ=255)                                                                                                                                                                                 cgcgTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGaaagg                                                                  <  1:111146/251‑4 (MQ=255)                                                                                                                                                                                  gcgTTGAAGTTAACGGAATCCAGGCGTGGTGGGCTCTCGATGCGAAAGAGATCCTCGGC‑CACGAAGCTGACCA‑GTATGTAAA‑AGTGCCGGATACTC‑TAGAAGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGCAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGTTCTCCACCGCGCTGaaaggc                                                                 >  1:377387/1‑247 (MQ=255)                                                                                                                                                                                   cgTTGAAGTTAACGGCATCCAGGTGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACTAAGCTGACTA‑GTATGTGAA‑AGTGCCTGATACTC‑TCGACGTATGGTTTGGCTTCGGATCTACCCTCTCTTCTGGTGTTGACGTGCGTCCGGAATTTGCTGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGAGCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTGCACCGCGATGaaaggca                                                                <  1:18002/251‑6 (MQ=255)                                                                                                                                                                                       gaagttaacggCATCCAGGCGTGGTTGGATCTCGATGCGAAAGAGATCCTCGGC‑TACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAg                                                             <  2:118720/241‑1 (MQ=255)                                                                                                                                                                                       gaagttaacggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAg                                                             <  1:310494/241‑1 (MQ=255)                                                                                                                                                                                       gaagttaacgacATCCAGGCGCGGTGGGATCTCGACGCGAATGAGATCCTCGGC‑GACGAAGCTGACCA‑GAATGTGAA‑AGTGCCGGATACTC‑TCGGCGTATGTTTTGACTCCGGATCTACCCAATCTTCTGGGGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGGCATGTATCTGGCAGGTTCTGACCAGCACCGTGAGTGGAGCATGTCTTCCCTCATGATCTCCACCTCGAttaaaggcaaagg                                                            >  2:97874/12‑239 (MQ=255)                                                                                                                                                                                           taaacggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCACGGC‑GAAAAAGCTGCACA‑GTATGTGAA‑AGTGCCGAAGACTC‑TCGACGTATGGATTGACTCCGGATCTACCAACTCTTCAGGTGTTTCCGTGCGTCCGGACATTGGCGGTCACGCCACCGACATGTATATGGAAGGTTCTGACCCACACAGAGGTAGGTTCAAGTCTTACCTAATGAACTCCCCCGCGAGGaaagacagagcgccg                                                        >  2:311613/6‑238 (MQ=255)                                                                                                                                                                                                ggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGAGCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGAATGGTTTGACTCCGGATCTACCAACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGAAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGAACTCCAGCGCTATGAAAGGCAAAGCGCCGTATCg                                                   >  2:266672/1‑251 (MQ=255)                                                                                                                                                                                                                gggATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GAATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTCTCTGGAAGGATCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACAGCGATGAAAGGCAAAGCGCCGTCTCGTCAGGCACTGAACCAc                                   >  2:279759/1‑251 (MQ=255)                                                                                                                                                                                                                   atcgcGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGGCCT‑GTTTGTGAA‑AGTGCCGGATACTC‑TCGGCGTATGGTTTGACTGCGGATCTACCCACTCTTCTGTTGTTGACGTGCGGCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGAGCTCCACCGCGATGAAAGGCAGAGCGCCGTATCGTCAGGTACTGTCCCACGGc                                <  2:195622/247‑1 (MQ=255)                                                                                                                                                                                                                       cGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTa                            <  1:81009/251‑1 (MQ=255)                                                                                                                                                                                                                         atacGAAAGAGATCCTCGGC‑GACGATGCTGACCA‑GTATTGGAA‑AGGTCCGGATACTC‑TTGACGTATGGTTTGACGCCGGATCTACCCACTCTGATGTTGTTGACGTGCGTCCTGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGATCAGCATCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCTCCGCGATGAAAGGCAAAGCGCCGTTGCGTCAGGTACTGACTTACGGTTTTAcc                          <  2:369161/248‑1 (MQ=255)                                                                                                                                                                                                                         aTGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTAcc                          <  1:88762/251‑1 (MQ=255)                                                                                                                                                                                                                         aTGCGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGGACTGCCCCCCGGCTTTAcc                          >  2:143539/1‑251 (MQ=255)                                                                                                                                                                                                                            cGAAAGAGATCCTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCCGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAATCGCCGTATCGTCAGGTACTTACCCACGGCTTTACCGTg                       >  2:145152/1‑251 (MQ=255)                                                                                                                                                                                                                                      ccgcGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAgg             <  1:167965/248‑1 (MQ=255)                                                                                                                                                                                                                                      ccTCGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGCGCCGGATACTC‑TCGACGTATGGCTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCAAGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCTCCGTAGCGTCATGGAATGCCCCACGGCTTTACCGTTGATGGTccgg             >  2:372466/1‑248 (MQ=255)                                                                                                                                                                                                                                         cGGC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGcc          <  1:361058/251‑1 (MQ=255)                                                                                                                                                                                                                                          ggC‑GACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGCCGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACAGTGGATGGTCAGGGCCg         >  2:101510/1‑251 (MQ=255)                                                                                                                                                                                                                                                  aaGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATg  >  1:73456/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       GCGTGCAGTGGATCCCGGACTGGGGCCAGGCGCGTATCGAGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGTGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGC‑GACGAAGCTGATCA‑GTACGTGAA‑AGTGCCGGACACAT‑TGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAGGGTAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATG  >  NC_000913/23709‑24199

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