Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,622	A→T	E274D (GAA→GAT)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,622	0	A	T	94.3%	136.9 / ‑1.1	53	E274D (GAA→GAT)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base T (14/36);  total (15/38)
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.21e-01

CAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGAC  >  NC_000913/3380‑3855                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            cAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑GCAAGCCGGATACCGCCCa                                                                                                                                                                                                                                     <  1:452226/251‑1 (MQ=255)     aTCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAgcg                                                                                                                                                                                                                                  <  1:36261/250‑1 (MQ=255)         ccgccAGGATTGCATTGCGCACGGGCGTCTTCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGGTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTg                                                                                                                                                                                                                             <  2:229175/251‑1 (MQ=255)           gacaggagtgcatttctcaatggCGACATCTGGCTGGCTTCATTCCCGCCTGCTAGTCCCGTCAGCCTGATCTTGCCGCGAATCTGATGAAAGATGTTATCGCGGAACCCTACCGTGAACGGTTACTTCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCTGTGCGTTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGcc                                                                                                                                                                                                                           <  1:248046/230‑1 (MQ=255)                  ttgcattgcGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGACAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCACGACCTTGTTCGCTCTGGGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTggg                                                                                                                                                                                                                    >  1:273976/1‑250 (MQ=255)                      attgcGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTggg                                                                                                                                                                                                                <  2:93988/251‑1 (MQ=255)                         gcgcACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATTGGCGCGGGAGCGAGCGGTGTCTCCGGCTCCGGGCCGACCCTTTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCACGCGCGTTGCCGACCGGtttgggt                                                                                                                                                                                                              >  2:248164/1‑245 (MQ=255)                             ccGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAAc                                                                                                                                                                                                         <  1:41443/250‑1 (MQ=255)                              cGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATCCGCGCGGGGGCGAGCGGTATCTCCGGCCCCGGCCCGCCCTTGTTCGCTCTGTGGG‑ACAAGCCGGGTACCGCCCCGCGCGCTGCCCGCTGGGTGGGGAAGAAc                                                                                                                                                                                                         >  2:396522/1‑250 (MQ=255)                                        ctggctggCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGGATCTCCGGCTCCCGCCCCACCTTTGTCGCTCTGTGGGGACA‑GCCGGATAACGCCCAGCGCGTTGTCGACTGGGTGGGGTAGAACTACCTCCaaa                                                                                                                                                                                               >  2:89776/1‑250 (MQ=255)                                         tggctggcTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGCAATCGGCGCGGTAGCGAGCGGTATCTCCCGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGGCCAGCGCGTTGCCGTCTGGTTGGGgaagaaataactggcaaat                                                                                                                                                                                             >  1:59038/1‑238 (MQ=255)                                          ggctggCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCaaa                                                                                                                                                                                               <  1:193852/244‑1 (MQ=255)                                             tggCTTCATTCACGCCTGATATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGCTGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCCGGCGGTTGAGGAAATCGGCGCGGTAGCGAGCGGTATATCCGGCTCCGGGCCGACCCTGTTCGCTCTTTGTG‑ACAAGCCGGTACCCGCCCCGCGCGTTGGCGCCCGGTTGGGGAAGAACTAACTCCAAACTCacg                                                                                                                                                                                         >  2:338859/2‑249 (MQ=255)                                              ggCTTCATTCACGCCGGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAg                                                                                                                                                                                          <  1:333563/249‑1 (MQ=255)                                                   cattcaCGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGACAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATCGGCCCGGGAGGGAGGGGGATCTCCGGCTCCGGCCCGCACTTGTTCTCTCTGTGGG‑ACAAACCAGATAACGCCCAGGGCGTTGCCGCATGGTTGGGGCAGAACTACCAGCCAACTCAGGGAGG‑tt                                                                                                                                                                                  >  2:221367/1‑251 (MQ=255)                                                   cattcaCGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTAGCGCTGAGCCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGTCCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑tt                                                                                                                                                                                  <  1:47213/251‑1 (MQ=255)                                                       caCGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTg                                                                                                                                                                               <  1:248164/250‑1 (MQ=255)                                                           ccTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTGATCGCGGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCat                                                                                                                                                                          <  2:129051/251‑1 (MQ=255)                                                                 aTTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACAACCTGCAAAATCAGGAAGG‑TTTTGTTCATAttt                                                                                                                                                                      <  1:30623/249‑1 (MQ=255)                                                                 aTTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCCGCGCGGGAGCGAGCGGGATCTCCCGCTCCGGCCCGACCCTGTTCGCTCTGTGGG‑ACAAGCCGGATACCGCCCCGCGCGTTGCCGGCCGGGTGGGGAAGAAACACCCGCCAAATCAGGAAGGGTTTTTTCCATATTTg                                                                                                                                                                     >  2:395888/1‑251 (MQ=255)                                                                    cccGGCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCgg                                                                                                                                                                 <  2:469198/251‑1 (MQ=255)                                                                            ccTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATCGGCGCGGTAGCGAGCGGGATCTCCGGCTCCGGCCCGCCCTTGTTCGCCCTGTGGG‑ACAAGCCGGATACCGCCCAGCGCGGTGCCCGCCGGGTGGGGAAAAACTTCCCGCCAAAACAGGGAGG‑GTTTGTTCATATTTTCCggacgggatc                                                                                                                                                         >  2:108617/1‑243 (MQ=255)                                                                             cTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGGAGCGAGCGGTATCTCCGGCTCCGGCCCGCCCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGGTGCCGACTGGTTGGGGAAGAACTACCTGCAAAATCAGGAAGG‑GTTTGTTCATATTTGCCGGCTGGATAcg                                                                                                                                                        >  1:282959/1‑251 (MQ=255)                                                                                   ttGCCGCGAAGCTGATGAAAGATGTTAGCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACggcggg                                                                                                                                                   <  2:10184/250‑1 (MQ=255)                                                                                    tGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACggcggg                                                                                                                                                   <  1:34373/249‑1 (MQ=255)                                                                                    tGCCGCGAAGCGGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGGG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACggcgggc                                                                                                                                                  <  2:4161/250‑1 (MQ=255)                                                                                         cgAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGa                                                                                                                                            <  1:383077/251‑1 (MQ=255)                                                                                          gAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGGATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGGTGCCGAATGGGTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑GTTTGTTCATATTTGCCGGCTGGATACggcgggaggcccgg                                                                                                                                           >  1:264586/1‑243 (MQ=255)                                                                                            aGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTa                                                                                                                                         <  1:232557/251‑1 (MQ=255)                                                                                             gCTGATGAAAGATGTTATCGGTGAACCCTACCGTGAACGGTTACGGCCTGGCGTCCGGCAGGCGCGGCTGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTAc                                                                                                                                        <  2:368541/251‑1 (MQ=255)                                                                                             gCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGGCTGGGTGGGGAAGAAATAACTGCAAAATAAGGAAGG‑GTTTGTTCATATTTGCCGGCTGGATACggcgggggccccgagac                                                                                                                                        >  1:197126/1‑240 (MQ=255)                                                                                                gatgaAAGATGTTAGCGCTGAACCCTTCCGTGAACGGTTACGGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACt                                                                                                                                       <  2:235354/249‑1 (MQ=255)                                                                                                           ttATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑GTTTGTTCATATTTGCCGGCTGGATACggcgggggggcggggagctggaaacacaaa                                                                                                                           >  1:98969/1‑226 (MQ=255)                                                                                                             aTCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACggcgggcggtcggcgtattggaaaataaaaca                                                                                                                         >  1:187649/1‑226 (MQ=255)                                                                                                               cGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGaaa                                                                                                                       <  1:396522/250‑1 (MQ=255)                                                                                                                    atcccgaccttGAACGGTTACTGCCTGGCGTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTAGCTCCGGCTCCGGCCCTACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTGCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTAc                                                                                                                 <  2:181350/241‑1 (MQ=255)                                                                                                                       ccTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAAt                                                                                                              <  1:40595/251‑1 (MQ=255)                                                                                                                          aCCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATc                                                                                                             <  1:395888/249‑1 (MQ=255)                                                                                                                           ccGTGAACGGTTACTGCCGGGCTGCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTTCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCCGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGa                                                                                                          <  2:469846/251‑1 (MQ=255)                                                                                                                             gTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGaaa                                                                                                        <  1:130505/251‑1 (MQ=255)                                                                                                                                    gTTACTGCCGGGCTTCCGTCAGGCGCGGCAGTCGGTTGCGGAAATCGGCTCGTGAGCGAGCGTTAGCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTTGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCTCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATca                                                                                                   <  1:285796/249‑1 (MQ=255)                                                                                                                                                  tCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAg                                                                                    <  1:89776/250‑1 (MQ=255)                                                                                                                                                                   gcgGTTGCGGAAATCGGCGCGGTAGCGGGCGGTACCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAAc                                                                  <  2:344302/251‑1 (MQ=255)                                                                                                                                                                                                 ggTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGGCCGAGGACTGGGAAACTAAATGAAAATCTACAATCTGAAAGATAACAAATAGCAGGGACGCTTTTCGCACACCGTAACCCAAGGGGTGGGGGAAAAA‑CAgggggtggt                                   >  1:68115/1‑246 (MQ=255)                                                                                                                                                                                                      ccccggcgccggcCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTc                              <  2:264586/249‑1 (MQ=255)                                                                                                                                                                                                        ccggctccggcCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTc                              <  1:130608/249‑1 (MQ=255)                                                                                                                                                                                                         cggctccggcCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGtt                                                        <  1:380208/222‑1 (MQ=255)                                                                                                                                                                                                                 gcCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTg                   <  1:124718/251‑1 (MQ=255)                                                                                                                                                                                                                          tGTTCGATTTGTGTG‑ACAAGCCGGATTCCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAAtt           <  1:58581/250‑1 (MQ=255)                                                                                                                                                                                                                              cGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAAt                                                                                                              <  2:115998/148‑1 (MQ=255)                                                                                                                                                                                                                                  gtgtgtg‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCctgac  <  2:237840/250‑1 (MQ=255)                                                                                                                                                                                                                                             aGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGc                  <  1:14997/225‑1 (MQ=255)                                                                                                                                                                                                                                               ccGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTg                                     <  2:264470/204‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            CAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGAC  >  NC_000913/3380‑3855

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: ‑