Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	23,986	T→C	D532D (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,986	0	T	C	92.0%	125.6 / 2.1	54	D532D (GAT→GAC)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (21/25);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.84e-02

AGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGTGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCA‑GTACGTGAA‑AGTGCCGGACACAT‑TGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAGGGTAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATG  >  NC_000913/23748‑24199                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                      aGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGAc                                                                                                                                                                                                           >  1:343441/1‑251 (MQ=255)    cggtGGTTTCTAACCGTCCTGACTGGTGTCACTCCCTGCAGCGCACCTGGGGCGTACCGATGTCACTGGTCGGTCACAAAGACGCGGAAGAGCTGCATCCGCGTCCCCTCGAACTGATGGAATAAGTTGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACTAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTcc                                                                                                                                                                                                        <  2:142012/248‑1 (MQ=255)     gctggtttatatccGTCCTGACTTGTGTATCTCCCGTCAGCGCACCTGGGGTGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGACAGAGATCCTCGGCGACGAAGCTGATCA‑GTACGTGAA‑AGTGCCGGACACAT‑TGGCAGTATGGTTTGAATCCg                                                                                                                                                                                                       >  2:251117/13‑251 (MQ=255)               aaCCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTAcc                                                                                                                                                                                               <  1:43918/249‑1 (MQ=255)                        gACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTTCCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCtg                                                                                                                                                                                     <  1:97079/250‑1 (MQ=255)                          cTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCC‑GTATGTGAA‑AGTGCCGGATACTC‑TCGCAGTATGGGTTGACTCCGGATCTACCCACTCTTCtgttg                                                                                                                                                                                  >  2:66719/1‑251 (MQ=255)                                       cGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTcc                                                                                                                                                                     <  1:130454/251‑1 (MQ=255)                                             cgcACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAAtt                                                                                                                                                               <  1:9669/251‑1 (MQ=255)                                                 ccTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGcc                                                                                                                                                           >  2:88762/1‑251 (MQ=255)                                                   tGGGGCGTACCGCTGTCACTGTTCGTGCCCAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGCACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCACGGCGACGAAGCTGA‑CACGTATGTGTA‑AGAGCCGGATA‑ATCTCGACGTATGGGTTGACTACGGATCTACCCACTCTTCTGTTGTTGCCTTTCGTCCGGAATTTGCCgg                                                                                                                                                         >  2:200705/1‑251 (MQ=255)                                                    gtggcgtACCGATGTCACTGTTCGGGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGGGTGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGGATGGTTTGACTCCGGATCTACCCACTCTGCTGTTGTTGACGTGCGTCCGGAATTTGCCg                                                                                                                                                          <  1:151906/244‑1 (MQ=255)                                                      ggcgttccattttctcTGTACGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCATCGAACTGATGGAAGAATTGGCAAAACGAGTTGAAGTTAACGGCTTCCAGGTGTGGTGGGTTCTCGTTGCGAAATAGATCCGCGGAGACGAAGCTGACAA‑TTATGTGAA‑AGTGCCGGATGCTC‑TCGGCGTATGGTTTGACTTCGGATGTACTCACTCTTCTGTTTTTGTCGTGCGTCCGGAATTTGTCGGtta                                                                                                                                                      <  2:287993/236‑3 (MQ=255)                                                         gtACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGAGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAACAGT‑CCGGATACAC‑ACCAAGTATGGTTTGACTCCGGATCTACCACCTCTTATGTTGTTGACGGGCGTACGGCATTTGCCGGTCACGc                                                                                                                                                   >  2:338675/1‑251 (MQ=255)                                                         gtACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAAATGATGGCAGAAGTGGCAAAACGCGTTGCAGTTAACGGCATACAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑CGTGCCGGAGACTC‑TCGGCGAATGGTTTGACTCCGGATCTACCCACTCTTCTTTTGTTGACGTGCATCCGGAATTTGCCGGTctctc                                                                                                                                                   >  2:322366/1‑247 (MQ=255)                                                          tACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGc                                                                                                                                                   >  1:369161/1‑250 (MQ=255)                                                             cGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCg                                                                                                                                               >  1:60791/1‑251 (MQ=255)                                                                   cACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATg                                                                                                                                         >  1:189274/1‑251 (MQ=255)                                                                       gTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATc                                                                                                                                     <  1:25955/251‑1 (MQ=255)                                                                          cGTGCACAAAGACACGGAAGAGCTGCATCCGCGTAACCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCC‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTgg                                                                                                                                  >  1:344671/1‑251 (MQ=255)                                                                            ttcacaAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGACGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTTGTGGGATATCGATGCGAAAGAGATCCTCGGCGACGATGCTGACCA‑TAATGTGAA‑AGTGACGGATACTC‑TCGCCGTATGGTTTGACTCCGGATATACCCACTCTGCTGCTGTTGACGTGCGTCCGGAGTTTGCAGGTAACGCAGCGGAccagtatccggt                                                                                                                                 >  2:265887/3‑239 (MQ=255)                                                                                caAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGGTGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGTCCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTGGACTCCGGATCTACCCACTCTTCTGTTGTTGACGGGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTACCTGGAAGGtt                                                                                                                            <  1:208257/251‑1 (MQ=255)                                                                                      cacGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGAcc                                                                                                                      <  1:94342/251‑1 (MQ=255)                                                                                              agctgcatccgtgggacctcggagtGATGGAAGAATTGGCAAAATGCGCTGAAGGTAACGGCATCCAGGCTGGGTGTGAGCTCGCTGCGAAAGAGTTCCTCGTCGACGAAGCAGACCA‑GTATGTGAA‑AGTGCCGGAGACTC‑TCGACGTATTGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGAAGTGCGTCCGGAATTTGCCGTTCACGCAGCGGACATGTATCTGGAAGGTCCTGACCAGCACcgt                                                                                                              <  2:305355/227‑2 (MQ=255)                                                                                              agCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTAACGCAGCGGACATGTATCTGGAAGGTTCTGAcccgcaccgt                                                                                                              >  1:221614/1‑243 (MQ=255)                                                                                               gCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACcgtg                                                                                                             <  1:312788/251‑3 (MQ=255)                                                                                                 tGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACcgtggt                                                                                                           >  1:170413/1‑247 (MQ=255)                                                                                                        gcgTACCATTGAAATGATGGAAGAAGTGGAAAAAAGAGTTGAAGTAGATGGTATCAAGGAGTGGTGGGATATAGATGAGAAAGAGATAATAGGAGAAGAAGATGATCA‑GTAAGTGAA‑AGTGCAGGAAAAAT‑TGGATGTATGGTTTGAATACGGATATAAAAAATATTATGTTGTTGACGTGCGTACGGAATTTGCAGGTAACGCAGCGGACATGTATCTGGAAGGTTCTGACCAACACCGAGGCTGGTTCa                                                                                                    <  1:251117/251‑1 (MQ=255)                                                                                                         cgaaacctcgAACAGATGGAAGAAGTGGGAAAACGCTGTGAAGTTAACGGCATCCAGGCGTGGCGGGATCTCGATGCGAAATAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGGATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGGCCGGAATTTGCTGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACcgtggtttgttcga                                                                                                   <  2:13978/242‑13 (MQ=255)                                                                                                                  gAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCt                                                                                          <  1:99796/251‑1 (MQ=255)                                                                                                                      tGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGAAGAAGCTGATCA‑GTACGTGAA‑AGTGCCGGACACAT‑TGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATg                                                                                      >  1:229218/1‑251 (MQ=255)                                                                                                                       gATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCA‑GTACGTGAA‑AGTGCCGGACACAT‑TGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGTCATGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATGa                                                                                     <  1:16331/251‑1 (MQ=255)                                                                                                                          ggAAGAAGTGGCAAAACGCGTTGAATTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGATGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGAtct                                                                                  <  1:242797/251‑1 (MQ=255)                                                                                                                              gaagTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTTATCA‑GTACGTGAA‑AGTGCCGGACACAT‑TGGATGTATGGTTTGACTACGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGGAGGTTCTGACCAACAACGCGGCTGGTTCATGTCTTCCCTAATGATCTCCAc                                                                              <  1:156967/251‑1 (MQ=255)                                                                                                                                         aCGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGaa                                                                     >  1:187875/1‑249 (MQ=255)                                                                                                                                          cgcgTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGaaagg                                                                  <  1:111146/251‑4 (MQ=255)                                                                                                                                           gcgTTGAAGTTAACGGAATCCAGGCGTGGTGGGCTCTCGATGCGAAAGAGATCCTCGGCCACGAAGCTGACCA‑GTATGTAAA‑AGTGCCGGATACTC‑TAGAAGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGCAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGTTCTCCACCGCGCTGaaaggc                                                                 >  1:377387/1‑247 (MQ=255)                                                                                                                                            cgTTGAAGTTAACGGCATCCAGGTGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACTAAGCTGACTA‑GTATGTGAA‑AGTGCCTGATACTC‑TCGACGTATGGTTTGGCTTCGGATCTACCCTCTCTTCTGGTGTTGACGTGCGTCCGGAATTTGCTGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGAGCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTGCACCGCGATGaaaggca                                                                <  1:18002/251‑6 (MQ=255)                                                                                                                                                gaagttaacggCATCCAGGCGTGGTTGGATCTCGATGCGAAAGAGATCCTCGGCTACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAg                                                             <  2:118720/241‑1 (MQ=255)                                                                                                                                                gaagttaacggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAg                                                             <  1:310494/241‑1 (MQ=255)                                                                                                                                                gaagttaacgacATCCAGGCGCGGTGGGATCTCGACGCGAATGAGATCCTCGGCGACGAAGCTGACCA‑GAATGTGAA‑AGTGCCGGATACTC‑TCGGCGTATGTTTTGACTCCGGATCTACCCAATCTTCTGGGGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGGCATGTATCTGGCAGGTTCTGACCAGCACCGTGAGTGGAGCATGTCTTCCCTCATGATCTCCACCTCGAttaaaggcaaagg                                                            >  2:97874/12‑239 (MQ=255)                                                                                                                                                    taaacggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCACGGCGAAAAAGCTGCACA‑GTATGTGAA‑AGTGCCGAAGACTC‑TCGACGTATGGATTGACTCCGGATCTACCAACTCTTCAGGTGTTTCCGTGCGTCCGGACATTGGCGGTCACGCCACCGACATGTATATGGAAGGTTCTGACCCACACAGAGGTAGGTTCAAGTCTTACCTAATGAACTCCCCCGCGAGGaaagacagagcgccg                                                        >  2:311613/6‑238 (MQ=255)                                                                                                                                                         ggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGAGCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGAATGGTTTGACTCCGGATCTACCAACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGAAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGAACTCCAGCGCTATGAAAGGCAAAGCGCCGTATCg                                                   >  2:266672/1‑251 (MQ=255)                                                                                                                                                                         gggATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GAATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTCTCTGGAAGGATCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACAGCGATGAAAGGCAAAGCGCCGTCTCGTCAGGCACTGAACCAc                                   >  2:279759/1‑251 (MQ=255)                                                                                                                                                                            atcgcGATGCGAAAGAGATCCTCGGCGACGAAGCTGGCCT‑GTTTGTGAA‑AGTGCCGGATACTC‑TCGGCGTATGGTTTGACTGCGGATCTACCCACTCTTCTGTTGTTGACGTGCGGCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGAGCTCCACCGCGATGAAAGGCAGAGCGCCGTATCGTCAGGTACTGTCCCACGGc                                <  2:195622/247‑1 (MQ=255)                                                                                                                                                                                cGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTa                            <  1:81009/251‑1 (MQ=255)                                                                                                                                                                                  atacGAAAGAGATCCTCGGCGACGATGCTGACCA‑GTATTGGAA‑AGGTCCGGATACTC‑TTGACGTATGGTTTGACGCCGGATCTACCCACTCTGATGTTGTTGACGTGCGTCCTGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGATCAGCATCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCTCCGCGATGAAAGGCAAAGCGCCGTTGCGTCAGGTACTGACTTACGGTTTTAcc                          <  2:369161/248‑1 (MQ=255)                                                                                                                                                                                  aTGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTAcc                          <  1:88762/251‑1 (MQ=255)                                                                                                                                                                                  aTGCGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGGACTGCCCCCCGGCTTTAcc                          >  2:143539/1‑251 (MQ=255)                                                                                                                                                                                     cGAAAGAGATCCTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCCGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAATCGCCGTATCGTCAGGTACTTACCCACGGCTTTACCGTg                       >  2:145152/1‑251 (MQ=255)                                                                                                                                                                                               ccgcGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAgg             <  1:167965/248‑1 (MQ=255)                                                                                                                                                                                               ccTCGGCGACGAAGCTGACCA‑GTATGTGAA‑AGCGCCGGATACTC‑TCGACGTATGGCTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCAAGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCTCCGTAGCGTCATGGAATGCCCCACGGCTTTACCGTTGATGGTccgg             >  2:372466/1‑248 (MQ=255)                                                                                                                                                                                                  cGGCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGcc          <  1:361058/251‑1 (MQ=255)                                                                                                                                                                                                   ggCGACGAAGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGCCGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACAGTGGATGGTCAGGGCCg         >  2:101510/1‑251 (MQ=255)                                                                                                                                                                                                          aaGCTGACCA‑GTATGTGAA‑AGTGCCGGATACTC‑TCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATg  >  1:73456/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                      AGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGTGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCA‑GTACGTGAA‑AGTGCCGGACACAT‑TGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAGGGTAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATG  >  NC_000913/23748‑24199

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