Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	23,983	G→C	L531F (TTG→TTC)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	23,983	0	G	C	93.2%	212.5 / 6.3	73	L531L (TTG→CTC)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base C (27/41);  total (29/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-02

AGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGTGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGAT‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTG‑GATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAGGGTAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCA  >  NC_000913/23748‑24215                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          gggcGATGGTTGCGAACCGGCCTGACTGGGGTAACTCCAGAAATCGCACCAGGGGCGGACCGATGTCACTGTACGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATTCAGTCGTGGTGGGATCTCGATACGAAAGAGAGCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTggc                                                                                                                                                                                                                           <  2:249766/248‑3 (MQ=255)      aTGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCgg                                                                                                                                                                                                                      <  1:470153/251‑1 (MQ=255)      aTGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCgg                                                                                                                                                                                                                      <  2:435776/251‑1 (MQ=255)              tAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGGATGTGAAAGTGCCGGACACTCTC‑GCCGAATGGTTTGCCTCCGGATCTcccc                                                                                                                                                                                                              >  2:154678/1‑251 (MQ=255)               aaCCGGCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTGCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACGGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAGC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACGCTC‑GACGTATGGTTTGACTCCGGATCTACCCa                                                                                                                                                                                                             <  2:411794/251‑1 (MQ=255)                aCCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCAc                                                                                                                                                                                                            >  1:238582/1‑251 (MQ=255)                          cTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGCTGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGCGATCCTCGGCGACGAAGCTGACCAGTATGTGACAGTGCCGGATACTCTC‑GCCGTATGGTTTGACTCCGGATCTACCCCCTCTTCtgttg                                                                                                                                                                                                  >  2:369942/1‑251 (MQ=255)                          cTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGCCCAGTATGTGAAAGTGCCGGATCCTCTC‑GACGTATGGTTTGACTACGGATCTACCCACTCTTCtgttg                                                                                                                                                                                                  >  2:420372/1‑251 (MQ=255)                            ggTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCGTCtgttg                                                                                                                                                                                                  <  1:226231/249‑1 (MQ=255)                                tATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACg                                                                                                                                                                                            >  1:306020/1‑251 (MQ=255)                                   ctcCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGc                                                                                                                                                                                         >  1:46217/1‑251 (MQ=255)                                     cccGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGt                                                                                                                                                                                       >  1:174440/1‑251 (MQ=255)                                        gTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCg                                                                                                                                                                                    <  2:83428/251‑1 (MQ=255)                                           aGCGCAACTGGCGCGTACCAGTGTCCCTGTTCGTGCACAAAGGCACGGCAGAGCTGGATTCGAGTACCCTCGAACTGAGGGAAGAAGGTGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCGGTTGTTGACGTGCGTCCGGaa                                                                                                                                                                                 <  2:118366/251‑1 (MQ=255)                                            gcgcACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTTCCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAAt                                                                                                                                                                                <  1:355593/251‑1 (MQ=255)                                              gcACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAAttt                                                                                                                                                                              <  1:231656/251‑1 (MQ=255)                                               ccccTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACAAGGAAGAGCTGCAGCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGAGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTg                                                                                                                                                                             <  2:225673/249‑1 (MQ=255)                                               cACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGACGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTACGGAATTTg                                                                                                                                                                             >  2:415489/1‑251 (MQ=255)                                                 ccTGGGGCGTACCGATGGCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGcc                                                                                                                                                                           <  1:360500/251‑1 (MQ=255)                                                        cgtacctaggtcactgtgccggCACAAAGACACGGAAGAGCTGCATCCTCGGACCCGCGAACTGATGGAAGAAGTGGCAAAGCGCGTTGAAGT‑TATCGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACTAATCTGACCAGTAGGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTGGATTCCGGATCTTCCCACTTTTCTTTTGTTGACGTGCGTCCGGAATTTGCCGGTCACg                                                                                                                                                                    <  2:463813/230‑1 (MQ=255)                                                         gtACCGATGTCACTGCTCGGGCTCAAGGACTCTGTAGTGTGGCATCCGCGTACACCCGAACTGATGGGAGAAGTGGCAAAGCGCGTTGAAGTTAAC‑TGCATCCAGGCGGGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGAGCCTGATACTCTC‑TTCGTAGGGTGCGACTCCGGAGCTACCCACTCGTCTGTTGGTGACGTGCGTCCGGAATTTGCCggctacgc                                                                                                                                                                   <  2:458354/251‑7 (MQ=255)                                                             cGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCg                                                                                                                                                               >  1:119022/1‑251 (MQ=255)                                                                 gTCACTGTTAGAGCACAAAGAAACGGAAGAGATGCCTCCGCGTACCCTCGACCTGCTGGAAGAAGTGGCAAACCGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GTCGTCTGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGCCGTGAGTCCGGAATTTGCCGGTCACGCAGCGGAcc                                                                                                                                                           >  2:500538/1‑250 (MQ=255)                                                                  tCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACAt                                                                                                                                                          >  1:213132/1‑251 (MQ=255)                                                                   cACTGTTCGTGCACAAATGCACCGAAGATCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACAt                                                                                                                                                          <  1:30699/250‑1 (MQ=255)                                                                   cACTGTTCGTGCACAAAGACTCGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGAGCTCGATGCGAAAGAGATCCTCGGCGACGAAGTTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACGCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATg                                                                                                                                                         <  2:177918/251‑1 (MQ=255)                                                                       ggtCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCGCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCGTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTa                                                                                                                                                       <  2:137850/247‑1 (MQ=255)                                                                         tCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTg                                                                                                                                                   >  1:289218/1‑251 (MQ=255)                                                                                caAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGCCCGGAATTTGCCGGTAACGCAGCGGACATGTATCTGGAAGGtt                                                                                                                                            <  1:441694/251‑1 (MQ=255)                                                                                 aaaGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCCTCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGCAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAACGGCCATGTATCTGGAAGGTTc                                                                                                                                           >  2:131690/1‑251 (MQ=255)                                                                                            agagCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGTAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GTCGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGCATTTGCCGGTCACGCAGCGGCCATGTATCTGGTAGGTTCTGAcccgcacc                                                                                                                                >  2:150548/1‑245 (MQ=255)                                                                                             gagCTGCATCCGCGTACCCTCGAACTGATGTAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACGCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCTGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACcg                                                                                                                               <  1:23998/251‑1 (MQ=255)                                                                                                  gCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACcgtggtt                                                                                                                          <  1:262738/251‑6 (MQ=255)                                                                                                    aTCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACcgtggttgg                                                                                                                        <  1:385456/251‑8 (MQ=255)                                                                                                           tACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCCGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATgg                                                                                                                 >  2:11691/1‑250 (MQ=255)                                                                                                            aCCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGt                                                                                                                 <  1:371679/250‑1 (MQ=255)                                                                                                            aCCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTc                                                                                                                >  1:7771/1‑251 (MQ=255)                                                                                                             cccTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTc                                                                                                                <  1:502111/250‑1 (MQ=255)                                                                                                             cccTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCt                                                                                                               >  1:366720/1‑251 (MQ=255)                                                                                                               cgcgAACTGATGGAAGAAGGGGCCAAACGCGTTGCAGGTTAC‑GGCATCCAGGCGTGGTGGGATCTATATGCGAAAGAGATCCTCGGCGAAGAAGCTGACAAGTATGTGAAATTGCCGGATACTCTC‑GACGTATGGTTTGCCTCCGGATCTACCCACTCTGCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCTGCACCGTGGTTGGTTCATGTCtgc                                                                                                             <  2:373539/248‑3 (MQ=255)                                                                                                                 cgAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTccc                                                                                                           >  1:8109/2‑251 (MQ=255)                                                                                                                  gAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTAccc                                                                                                                                                                                                              <  1:470990/151‑1 (MQ=255)                                                                                                                  gAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTAcac                                                                                                                                                                                                              >  2:470990/1‑149 (MQ=255)                                                                                                                             agaagTGGCAAAACGCGTTGCAGTTAAC‑GGCATCCAGGCGTAGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGGGCCGGATACTCTC‑GACGTATGGATTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACAGTGGTTGGTTGATGTATTCCCTAATGATCTccc                                                                                               >  2:333189/1‑250 (MQ=255)                                                                                                                              gaagTGGCAAAACGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCAc                                                                                              <  1:73381/251‑1 (MQ=255)                                                                                                                                       aaaCGCGTTGAAGTTAAC‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGaa                                                                                     <  1:170297/251‑1 (MQ=255)                                                                                                                                           cgcggtgaaattaacggCATCCAGGCGTGGTGGGATCTCGATGCGAAATATATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCGTCTTTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTGCATGTCTTCCCTAATGATCTCCACCGCGATGaaagg                                                                                  <  2:420144/236‑4 (MQ=255)                                                                                                                                              ttGAAGTCGAT‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTG‑GATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGt                                                                                                                                             <  2:144464/188‑1 (MQ=255)                                                                                                                                              ttGAAGTCGAT‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTG‑GATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGt                                                                                                                                             >  1:144464/1‑188 (MQ=255)                                                                                                                                              gttgaagttaacggCATCCAGGCGGGGTGGGATCTCGATGAGAAAGCGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGCCATGGATCTGGAAGGTTCTGACCCGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGCAGAAgggca                                                                                >  2:405826/13‑248 (MQ=255)                                                                                                                                                    gttaacggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGcc                                                                         <  1:375991/245‑1 (MQ=255)                                                                                                                                                       aacggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGt                                                                       >  1:326704/4‑250 (MQ=255)                                                                                                                                                          ggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCg                                                                   <  1:426267/251‑1 (MQ=255)                                                                                                                                                           gCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGt                                                                  <  1:286426/251‑1 (MQ=255)                                                                                                                                                                    cGTGGTGGGATCGCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTGCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTg                                                         <  1:293340/251‑1 (MQ=255)                                                                                                                                                                          gggagcccGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCAc                                                   <  2:158392/244‑1 (MQ=255)                                                                                                                                                                            gATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGGGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGGCTTCCCTAATGATCTCAACCGCGATGAAAGGAAAAGCGCCGTATCGTCAGGTACTTACCAACgg                                                 >  2:368565/1‑251 (MQ=255)                                                                                                                                                                                  gATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTG‑GATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTc                                                                                                                                           <  1:201754/155‑1 (MQ=255)                                                                                                                                                                                  gATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGCCATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCtttcc                                           >  2:51790/1‑249 (MQ=255)                                                                                                                                                                                  gATGCGAAAGAGATCCTCGGCGAAGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTG‑GATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTc                                                                                                                                           >  2:201754/1‑155 (MQ=255)                                                                                                                                                                                        aaaGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGa                                     <  1:370176/251‑1 (MQ=255)                                                                                                                                                                                         aaGAGATCCTCGGCGACGAAGCTGCCGAGTATGTGAAAGTGCCGGATTCTCTC‑GACGTATGGTTTGACTCCGGACCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATTTATCTGGAAGGTTCTGACCAGCGCCGTGGTTGGTTCATGGCTTTCTTGATGATCGCCTCCGTGATGAAAGGCAAGGCGCCGTAGCGTCAGGATCTGATCCAGGGCTCTACCGTGGAt                                    <  2:211957/251‑1 (MQ=255)                                                                                                                                                                                          agagATCCTCGGCGACGAAGCTGACCAGTAAGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCTTCAGGTACTGATCCACGGCTTTACCGTGGATg                                   <  2:289218/251‑1 (MQ=255)                                                                                                                                                                                            agaTCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGt                                 <  1:481228/251‑1 (MQ=255)                                                                                                                                                                                                       gacgaAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGAGGATGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGAACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCa                       >  2:384546/1‑250 (MQ=255)                                                                                                                                                                                                        acgaAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCa                       >  1:259883/1‑249 (MQ=255)                                                                                                                                                                                                           aaGCTGACCAGTATGTGAAAGTGCCGGATAC‑TCGCGACGGATGGTTTGACTCCGGATCTACCGACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGGCTTCCCTAATGATCTCCACCGCTATGAAAGGCAAAGCGCCGTATCGGCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATg                  <  2:56897/251‑1 (MQ=255)                                                                                                                                                                                                             gCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTc                <  1:427841/251‑1 (MQ=255)                                                                                                                                                                                                                gaccAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAAGTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTTACCCACGGTTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTaa             <  1:191976/248‑1 (MQ=255)                                                                                                                                                                                                                   cAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACGCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTTTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCTGCAAGATGTCTAAATc          <  1:409289/251‑1 (MQ=255)                                                                                                                                                                                                                    aGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATcc         <  2:358757/251‑1 (MQ=255)                                                                                                                                                                                                                         gTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACAGGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGAGCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGTCCCACGGCTTTACCGTGGATGGTTAGGGCCGCAAGATGTCTAAATCCATCgg    <  1:444009/251‑1 (MQ=255)                                                                                                                                                                                                                           gAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGCCATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGcc  >  2:253612/1‑250 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          AGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGTGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGAT‑GGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTG‑GATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAGGGTAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCA  >  NC_000913/23748‑24215

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