Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,763	C→T	N10N (AAC→AAT)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,763	0	C	T	91.7%	90.4 / ‑5.6	38	N10N (AAC→AAT)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (9/24);  minor base A (3/0);  total (12/26)
Fisher's exact test for biased strand distribution p-value = 3.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGAT  >  NC_000913/3530‑4009                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                           caggcgcggcaggcgGTTGCGTAAATCGGCTCGTTAGCGGGCGGTATCTCCTGCTCCGGCCCGACCTTGTTCGCGCTGTGTTACAAGCCGGATACCGCCCAGCGCGTTGCCGTCTGGTTGGGTAAGATCTACCTGCAAAATCAGGAAGGTTTTGTGCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAACGAATCTTT‑CCAATCTGAAAGATCGCAATGAGCAgggaatctgtgc                                                                                                                                                                                                                                           <  2:281530/251‑11 (MQ=255)                  tgcgGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAg                                                                                                                                                                                                                          <  1:31102/249‑1 (MQ=255)                  ggggcaaaTCTGCGCGGGAGCGAGCGGTATCTCCGGCTCCGGGTCGGCCTTTGGCGCTCGGTGTGACAAGCCGGGTACTGTCCAGGGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAAGCAGGAAGGTTGTGTTCACATGTGCCGGCTGGATACGGCGGGCACACGAGTACTGGAAAACTGAATGAAACTCT‑ACAATCGGAAAGATCACAATGAGCAGGTCAGCTTGGCGTAAGCCGTAA‑CCCAgg                                                                                                                                                                                                                         <  2:352715/246‑1 (MQ=255)                           cggcgctggagtgagcgGTATCTCCGGTCCCGGCCCGATCTTGTTCGGTCTTTGTGACAGGCGGTATCCCGCCCGGTGCGGTGCCTGCGGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTGTTGTTCATACTTGCCGGCTGGATACGGCGTGCGCACGAGTACTGGAAAACTAAATGAAACACT‑ACAATCTGAAAGAGCGCAAGGAGCAGGGCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCa                                                                                                                                                                                                               <  2:447401/239‑1 (MQ=255)                                 ggTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATc                                                                                                                                                                                                         <  1:637205/251‑1 (MQ=255)                                            gTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑ttt                                                                                                                                                                                             <  1:315958/251‑1 (MQ=255)                                                  tctgcaccggcCCGTCCTTGTTCGCTCTTTGGGACAATCCGGATTCCGCCCAGCGCGTTGCCGACTGGTTGTGTAAGAACTACCTGCAAAAGCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAACCTCT‑GCAATCTTAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑ttttt                                                                                                                                                                                           <  2:557738/241‑1 (MQ=255)                                                                     gttggctatgtgtgACAAGCCGGATACAGCCTATCGTGGTGCCGACTGGGTGGGTAAGAACTACTTGCAAAACCAGGAAGGTTTTGATCATATTTGCCGGCTGGATACGGCGGGCGCAAGATTACGGGAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAATGAGCAGGCCAGCGTTGCGCAAGCCGGAA‑CCCAGGGG‑GTGGGCAAAAATCATGGGCTG‑TGTTTTCCGCACGACCTGCCTGAATTc                                                                                                                                                                     <  2:190210/242‑1 (MQ=255)                                                                            ctGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCctg                                                                                                                                                               <  1:399093/249‑1 (MQ=255)                                                                                           gataCCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGGGGACTGGAGAAATAAATGAAACTCT‑ACAATCTGAAAGATCACAAAGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAAGGG‑GTGGGGAAAAAACAGGGGGGG‑GTTTTTTCGCCCGACCCCCCGGGATTCCGCCTGCCTGGAA‑ttagtgag                                                                                                                                               >  1:343832/4‑243 (MQ=255)                                                                                                 gCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCT‑ACAAGCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGTCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGaa                                                                                                                                       <  1:175156/251‑1 (MQ=255)                                                                                                    catcgcgTTCCCGACTGGTAGGGAAATATATACTTGTAAAAACAGGAAGGGTTTGTGCATATTCACCGGGTGGATACGGCGGGCGCACGAGTACTGGAAAACTAATTGAAGCACG‑ACAATGAGAAAGATCACAAGAAGAATGTCAGCTTTGCTCAAGCCGTAA‑CGCAGGGG‑TTGGGCAAAAATCAGGGACTG‑TATTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGca                                                                                                                                    <  2:197617/248‑2 (MQ=255)                                                                                                       cgcgTTGCCGCCTGGTTGGGTAAGAACTACCTGCCAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGATTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGa                                                                                                                                 <  1:161999/251‑1 (MQ=255)                                                                                                                                  tACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGGGGACTGGAAAAATAAATAAAAATCTCACAA‑CTTAAAGATCACAATGAGCAGGGCAGCTTTGTGCCAGCCGTAA‑ACCAGGGGGTGGGGAAAAAAACAGGGG‑TG‑TTTTTTTTCCACGACCTTCCCGGATTCAACCTGACTGAAAATTG‑TGAGAAGGTGAAGCTGGAgttttttaccccccgggcggaggaccc                                                                                                      >  1:429231/1‑228 (MQ=255)                                                                                                                                   aCCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATGTTGTCACCCGCAGTGCGAAGATCct                                                                                                      <  1:387441/250‑1 (MQ=255)                                                                                                                                       gCAAAATCAGGAAGTTGTGGTTCATATTGGGCGGCTGGATACGGCGGGCGCACGAGTACTGGAAACCTAAATGAAACTCT‑ACAATCTTTAAGATAACAATGAGCAGTTCAGCTTTGCGCATGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTGTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGTTGAAGCTGGATTTTGGCTCCCGCAGTGCGAAGATCCTCTCg                                                                                                  <  2:608267/250‑1 (MQ=255)                                                                                                                                       gCAAAAGCAGGAAGGGTTTGGTCATAGTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAAGCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCtctc                                                                                                   <  1:615779/249‑1 (MQ=255)                                                                                                                                           aataaGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACGGGAAAACTAAATGAAACTCT‑ACAAGCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGTG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGtt                                                                                             <  2:259656/247‑1 (MQ=255)                                                                                                                                                 gAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTgg                                                                                       <  1:500011/251‑1 (MQ=255)                                                                                                                                                      ttttGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGACTGGGAAACTAAATTAAACTCT‑ACAAACTGTAAAATAACAAAGAGCAGGACAGCTTTTCGCAAGCCGTAAACCCGGGGG‑TGGGGCAAAA‑TAAGGGGGGGGTTTTTTCC‑CCCGACCTGCCCGGGTTTAGCCCGCCTGAAA‑ATGATGAGAGGCTGAAACTTGGTTTTTTCACCCCCCGGGGGAAGATCCCCTCCGCGTTTTTTGGtgttg                                                                                  >  1:333529/1‑248 (MQ=255)                                                                                                                                                          gtttatatTTGCCGGCGGGATACCGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCG‑CCAATCTGAAAGATCCCAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGTGCTG‑TTTTTTCCGCACGTCCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGAAGCTGAAGCTTGATTTTGTCTCCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGaaa                                                                               <  1:558208/246‑1 (MQ=255)                                                                                                                                                                       ggCTGGATACGGCGGGCGCACGATTACTGGAAAACTAAAGGAAACTCT‑ACAATCTGAAAGAGCACAATGAGCAGGGCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCGG‑TTTTTTCCGCACGACCTGCCGGAATGCAGCCTGACTGAAA‑GTGATTAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGAGCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATc                                                                 <  2:568463/251‑1 (MQ=255)                                                                                                                                                                                 ggcgggcgCACGAGTACTGGAAAACTAAATGATACTCT‑ACAATCTGAAAGATCACAAAGAGCCGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGAACGACCTGCCGGAATTCAGCCTGAATGAAA‑TTGATGAGATGCTGAAGATGGATTTTGTCACCCGCAGGGCGAAGATCCTCTCGGCGTTTATGGGTGATCAAATCCCGCAGGCAATCATGGacaggg                                                       >  2:452767/1‑246 (MQ=255)                                                                                                                                                                                 ggcgggcgCACGAGTACTGGAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAgc                                                       >  1:36548/1‑251 (MQ=255)                                                                                                                                                                                    gggcgCACGAGTACTGGAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGcgcg                                                    <  1:134337/251‑1 (MQ=255)                                                                                                                                                                                              gTACTGGAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGACGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCgcggc                                          >  1:469803/1‑251 (MQ=255)                                                                                                                                                                                                aCTGGAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCctgactg                                                                                                                                                           <  1:584822/137‑1 (MQ=255)                                                                                                                                                                                                aCTGGAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGACAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCctgactg                                                                                                                                                           >  2:584822/1‑137 (MQ=255)                                                                                                                                                                                                   ggAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGtt                                       <  1:442421/249‑1 (MQ=255)                                                                                                                                                                                                                aaaCTCT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTcc                        >  1:120943/1‑251 (MQ=255)                                                                                                                                                                                                                aaaCTCT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTCTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTcc                        <  1:582599/251‑1 (MQ=255)                                                                                                                                                                                                                aaaCTAT‑ACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGCTCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTcc                        >  2:149845/1‑251 (MQ=255)                                                                                                                                                                                                                               aaaaGATCACAATGAGCAGGTCAGCGTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAAGCAGGGGCTG‑TTTTTGCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCGGGAGTTTGTCACCCGCAGTGCGAAGATCCGCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTg          <  1:289172/250‑1 (MQ=255)                                                                                                                                                                                                                                aaaGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCAGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAGGATCCTCTCGGCGTTTATTGGCGATGAAATCCCCCAGGAAATCCTGGAAGCGCGCGTGCGCGCGGCGTTTGCCTTCCAGGCTCCGGTCGCCAATGTTg          >  2:238623/1‑250 (MQ=255)                                                                                                                                                                                                                                      cacaATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTg          >  1:591959/1‑244 (MQ=255)                                                                                                                                                                                                                                      cacaATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTg          <  2:591959/244‑1 (MQ=255)                                                                                                                                                                                                                                      cacaATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCg    >  2:286167/1‑250 (MQ=255)                                                                                                                                                                                                                                       acaATGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGAt  <  1:173170/251‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                           CAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCT‑ACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGG‑TTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAA‑TTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGAT  >  NC_000913/3530‑4009

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