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	98.5%	179.5 / NA	67	E274D (GAA→GAT)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (24/42);  minor base G (0/1);  total (24/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGG‑TCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAG  >  NC_000913/3380‑3868                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                           cAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGATGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGCACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑GTGTGGAACTCGGCGCGGTAGCCAGCGGTATCTCCTGCTCCGGCCCCACCTTGTTTGCTCTGTGTGACAAGCCGGATACCGcccc                                                                                                                                                                                                                                                    >  1:353812/1‑250 (MQ=255)        gccgccAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTTCGGGAATCGGCGCGGTTGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTTTTCGCTCTGTGTGACAAGCCGGaatcccccccgcccgtt                                                                                                                                                                                                                                             >  2:615222/1‑236 (MQ=255)           gccAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTCGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGcc                                                                                                                                                                                                                                         >  1:453054/1‑251 (MQ=255)               ggATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGGAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACt                                                                                                                                                                                                                                     >  1:511198/1‑251 (MQ=255)                gcttgaattgcGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGGAAACGGCGCGGTAGCGAGAGGGATCTCCGGCTCCGGCCCGACCTTGTTTGCCCTTTGTGACAAACCGGGTAACGCCCAG‑GCGCTTGCCGACTg                                                                                                                                                                                                                                    >  2:408943/7‑251 (MQ=255)                gATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTg                                                                                                                                                                                                                                    <  2:285957/251‑1 (MQ=255)                 attgcattgcGCACGGGCGACCTCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGCCCTTGTTCGCTCTGTGTGCCAAGCCGGATACCGCCCCGCGCG‑TTGCCGACTgg                                                                                                                                                                                                                                   >  1:207131/1‑251 (MQ=255)                     cattgcGCACGGGCGACATCGTGCTGGCTTCATTCACGCCTGACATTCCCGTCAGCCTGAGCTTGCCGCTAAGCTGATGAAAGATGTTAGCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCACGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTgg                                                                                                                                                                                                                               <  2:165440/251‑1 (MQ=255)                      attgcGCACGGGCGACTGCTGGCTGGCTTCATTCACGCCTGTTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGAGGAAAGATGTTATCGCTGAGCCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCAGCTCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCA‑TTGCCGACTGGTTggg                                                                                                                                                                                                                              <  2:459118/251‑1 (MQ=255)                       ttgcGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCCGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGt                                                                                                                                                                                                                             >  1:297839/1‑251 (MQ=255)                                ggCGACATCTTGCTGGCTTCATTCACTCCTGCTATTCCCGGCATCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTGCTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATTCCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTAc                                                                                                                                                                                                                    <  2:503356/251‑1 (MQ=255)                                ggCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGGTCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTAc                                                                                                                                                                                                                    <  1:60601/251‑1 (MQ=255)                                   gACATCTGGTTGGCTTCATTCACGCCTGCTATTCCCTGCAGCCTGAGCGTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTTACAAGCCGGATACCGCCCAGCGCG‑TTGTCGACTGGTTGGGTAAGAACTACCTg                                                                                                                                                                                                                 <  2:580093/251‑1 (MQ=255)                                     cATCGGGCTGGCTTCATTCACGCCTGCTAGTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACGGTTTGGGTAAGAACTACCTGCa                                                                                                                                                                                                               <  2:207131/251‑1 (MQ=255)                                      aTCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGGTGGGTAAGAACTACCTGCa                                                                                                                                                                                                               >  1:493145/1‑250 (MQ=255)                                          ggctggCTTCATTCACGCCTGCTATTCCCGGCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCaaa                                                                                                                                                                                                             <  2:300670/244‑1 (MQ=255)                                              ggCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGa                                                                                                                                                                                                      <  1:164701/251‑1 (MQ=255)                                                      tcaCGCTTGCTATTCACGTCGGCCTGAGATTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTCCCGTGAACGGTTACTGCCTGGCCTCCGGCAGGCGCGGCAGGCGGGTC‑AGGAAACCGGCGCGGTCGCGCGCGGTACCTCCGGCCCCGGCCCGACCTTGTTCGCTCTGTGTGACACGGCGGATTCCGCCCCCCGCG‑TTGCCCACTAGTTGGGGCAGCACTACCTTCaacattcggcagggtttg                                                                                                                                                                                              >  2:230302/1‑235 (MQ=255)                                                        aCGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCCAAATCAGGAAGGTTTTGtt                                                                                                                                                                                            >  1:365319/1‑251 (MQ=255)                                                         cGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGGATCTCCGGCTCCGGCCCGGCCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCCGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTTCAAAAACAGGAAGGGTTTGttt                                                                                                                                                                                           >  2:321548/1‑250 (MQ=255)                                                          gCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTACGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCTCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGGAAGAACTCACTGCAAAATCAGGAAGGTTTTGTTc                                                                                                                                                                                           >  1:583889/1‑250 (MQ=255)                                                           ccTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTc                                                                                                                                                                                           <  2:529923/249‑1 (MQ=255)                                                              gCTATTCCCGTAAGGTTGGGCTTGCCGCGAAGCTGAGGAAAGATGTGATCGCTTAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCTGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCATGAAGGTTGTGTTCata                                                                                                                                                                                        <  2:303466/249‑1 (MQ=255)                                                               cTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACCAGCCGGGTACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTTAGAACTACCTGCAAAATCCGGAAGGTTTTGTTCATtttt                                                                                                                                                                                     >  1:8070/1‑251 (MQ=255)                                                                tATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TAGCGGAAATCGGCGCGGGAGCGAGCGGTATCTCCGGCTCCGGACCGACCTTGTTCGCTCTGTATGACAAGCCGGATAACGCCCAGCGCG‑TTGCCGACTGGATGGGAAAGAACCACCTGGAAAATCAGGATGGGTTTGTTCATATTTg                                                                                                                                                                                    >  1:377353/1‑251 (MQ=255)                                                                 aTTCCCGTCTGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCCTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTg                                                                                                                                                                                    <  1:83417/250‑1 (MQ=255)                                                                    cccGTCAGCCTGAGCTTGGCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTTAACGGTTACTGCCTGGCTTCCGGCAGTCGCGGCAGGCGG‑TGTCGGAAATCGGCGCTGTAGCGAGCGGTAGCTCCGGCTCCGGCGCGAGCTTGTTCGCTCTGGGGGACAAGCCGGATACCGCTCAGCGCG‑TTGCCGACTGGTTGGGTAAGAATTACCTGCAAAATCAGGAAGGTTTTGTTCATACTTGCCg                                                                                                                                                                                 <  2:230513/250‑1 (MQ=255)                                                                    cccGTCAGCCTGAGCTGGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCGTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGGAGCGAGCGGTATCTCCGGCTACGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGACCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGAAAAATCAGGAAGGGTTTGTTCATATTTGCCgg                                                                                                                                                                                <  2:250858/251‑1 (MQ=255)                                                                       gTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTg                                                                                                                                                                             <  2:267448/251‑1 (MQ=255)                                                                         cAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGGAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGGAAGAACTACCTGCAAAATCAGGAAGGGTTTGTTCATATTTGCCGGCTggg                                                                                                                                                                           >  1:593496/1‑250 (MQ=255)                                                                             cTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGGTAcg                                                                                                                                                                       >  1:219377/1‑251 (MQ=255)                                                                                aGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGAATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATAcggcg                                                                                                                                                                    >  1:535405/1‑251 (MQ=255)                                                                                    tGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACggcggg                                                                                                                                                                  <  1:273032/249‑1 (MQ=255)                                                                                        gcgAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACg                                                                                                                                                            <  2:197491/251‑1 (MQ=255)                                                                                         cgAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGa                                                                                                                                                           <  2:439557/251‑1 (MQ=255)                                                                                         cgAAGCTGATGAAAGATGTTAACGCTGAACCCTACCGCGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGAGGTTGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCCGCGCG‑TTGCCGACTGGTTGGGGAAGAACTCCCTGCAAAAACAGGAAGGTTTTGTTCATATTTGCCGGCCGGATACggagggggcgcgg                                                                                                                                                           >  1:145025/1‑244 (MQ=255)                                                                                           aaGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGAGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGGTGGGGAAGAAATACCTGCAAAATAAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACggggggcgggcgggg                                                                                                                                                         >  1:443224/1‑244 (MQ=255)                                                                                                gatgaAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑gg                                                                                                                                                   <  1:479853/251‑1 (MQ=255)                                                                                                 atgaAAGATGTTATCGCTGAACCCTACCGTGAACGGTTAATGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCCCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCAACCTTGTTCGCTCTGGGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGGCTGGTTGGGTGAGAACTACCTGCAAAATCAGGACGGATTTGTTTATCTTTGCCAGCTGGTTACGGGGGGGGCAcaacgatgggg                                                                                                                                                   >  1:356223/1‑241 (MQ=255)                                                                                                        aTGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTaa                                                                                                                                           <  1:165209/251‑1 (MQ=255)                                                                                                          gTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTAGCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAAt                                                                                                                                         <  2:575267/251‑1 (MQ=255)                                                                                                           ttATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGGAGCGAACGGTATCCCCGGGTCCGGCCCGACCTTGGTCGCTCTGTGTGACAAGCCGGATAACGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTTCCGGCTGGCTACggcggggggacacggtgagtgaaaaactaa                                                                                                                                          >  2:593507/1‑226 (MQ=255)                                                                                                                     aCCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACa                                                                                                                             <  1:473555/251‑1 (MQ=255)                                                                                                                      cccTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACaa                                                                                                                            <  1:7247/251‑1 (MQ=255)                                                                                                                             gTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGGATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGGAAGAACTACCTGCAAAATCAGGAAGGGTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGGCCGCGGACTTGGAAAA‑TAAATGAAA‑ATCTACaaacttca                                                                                                                      >  1:320473/1‑244 (MQ=255)                                                                                                                               gAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGa                                                                                                                   <  1:321548/251‑1 (MQ=255)                                                                                                                                aaCGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGAt                                                                                                                  <  1:117452/251‑1 (MQ=255)                                                                                                                                    gTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGATcaca                                                                                                              <  2:343246/251‑1 (MQ=255)                                                                                                                                    gTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGATcaca                                                                                                              <  2:365319/251‑1 (MQ=255)                                                                                                                                    gTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGATcaca                                                                                                              <  1:348670/251‑1 (MQ=255)                                                                                                                                                  tCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGATCACAATGAGCAGGTCa                                                                                                  <  1:246053/249‑1 (MQ=255)                                                                                                                                                                 aggcgG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTa                                                                                 <  1:621000/251‑1 (MQ=255)                                                                                                                                                                    cgG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGGCCGAGGACT‑GGGAAACTAAAAAAAA‑ATCTTAAATATGAAAAATCACAATGAGCAgggacacctttcgcaaaccgtaaac                                                                              >  1:101763/1‑228 (MQ=255)                                                                                                                                                                    cgG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAAcc                                                                              <  1:499070/251‑1 (MQ=255)                                                                                                                                                                       gttgcgGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCC‑a                                                                           <  2:493145/248‑1 (MQ=255)                                                                                                                                                                          gcgGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCC‑Agg                                                                         <  1:145299/250‑1 (MQ=255)                                                                                                                                                                                     cgcgGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CGCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCC‑AGGGGTTGGGCaaaa                                                             <  1:615222/251‑1 (MQ=255)                                                                                                                                                                                      gcgttagcgtggggtttctcCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAGACCGCCCAGCGTG‑TTGACGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGCTCATATTTGCCGGCTGGATGCGTCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACTATCTGAAAGACCGCAATGAGCAGGTCAGCTTTGCGCAAGCCGTATCAC‑AGGGGGTGGGCaaaaa                                                            <  2:528377/235‑1 (MQ=255)                                                                                                                                                                                                   gTATCTCCGGCTCCGGCCCGACTTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGGAACCC‑AGGTGTTGGGCAAAAATCAGGGGCTGttt                                               <  2:593496/251‑1 (MQ=255)                                                                                                                                                                                                         ccggctccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTcc                                          <  2:156231/250‑1 (MQ=255)                                                                                                                                                                                                            gctccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGGGGAAT‑GGAAAAATAAATGAAA‑CACTACAATCTTGAAGATCACAAAAAGCAGAGCAGCTTTGCGCAAACCCTAACCC‑AGGGGGTGGGGAAAAATCAAGGGATGGTTTTTTCGCAc                                      >  2:410505/1‑251 (MQ=255)                                                                                                                                                                                                               ccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACg                                     <  1:224215/249‑1 (MQ=255)                                                                                                                                                                                                               cccgcacttccTTGGTCGTTCTGTGGGCCAAGCCGGATACCGCCGAGCGGG‑TTGTCGACTGGTTGGGTAAGAACTACCTGCAAAATCGGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CCCTACAATCTGAAAGAGCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCC‑AGGGTTTGGGCAAAAATCAGGGGATGTTTTTTCCGCACGAc                                   <  2:461903/242‑1 (MQ=255)                                                                                                                                                                                                                   cccGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGc                               <  1:119443/251‑1 (MQ=255)                                                                                                                                                                                                                        ccTTGTTCGCTCTGTGTGACAAGCCGGATACCGCC‑AGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAACATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCCCGGGTACT‑GGAAAAATAAATGAAAACACT‑CAATCTTAAAGATCACAATCAGCCGGTCAGATTTCCGCAAGCCGGAACCCCGGGGGTTGGGCACAA‑TCAGGGGCTGGTTTTTTCGCACGACCTGCCggcat                         >  1:291432/1‑248 (MQ=255)                                                                                                                                                                                                                                gCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCC‑AGGGGTTTGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCct                  <  1:596710/251‑1 (MQ=255)                                                                                                                                                                                                                                        ggACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CGCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGaaa          <  2:453054/250‑1 (MQ=255)                                                                                                                                                                                                                                                cGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGTGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATgag  <  2:587062/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                           CAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGG‑TCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACT‑GGAAAACTAAATGAAA‑CTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAG  >  NC_000913/3380‑3868

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