Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,420	T→A	A229A (GCT→GCA)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,420	0	T	A	97.2%	104.9 / NA	37	A229A (GCT→GCA)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (19/16);  minor base G (0/1);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGT‑TGCGCAGCTGCCGCAGGAGA‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTG‑AACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGC  >  NC_000913/4210‑4664                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        aGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTg                                                                                                                                                                                                               >  1:13772/1‑251 (MQ=255)      aaaTCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCg                                                                                                                                                                                                          >  1:78451/1‑251 (MQ=255)                     aaGAAAAGCTGTTCTGGACATTGGGCAGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGGCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGGGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACCTCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCTGGAAC‑CGAGCCACCAGCCGGTTGTCTCGGTGCCAAGCTGAAAc                                                                                                                                                                                           >  2:75739/1‑251 (MQ=255)                          aaaCTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAg                                                                                                                                                                                                                                      <  2:279666/204‑1 (MQ=255)                          aaaCTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAg                                                                                                                                                                                                                                      >  1:279666/1‑204 (MQ=255)                          aaaCTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCgg                                                                                                                                                                                      >  1:441408/1‑251 (MQ=255)                           aaCTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGCTTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGc                                                                                                                                                                                     >  1:455692/1‑251 (MQ=255)                                        atggggCGGCAATATCGAAACTGTTGCCATCGACGGCGATTGCGATGCCTGTCAGGCGCTGGGGAAGCATGCGTTTGATGATGAAGAGCGTAAAGGGGCGCTGGGGTTAAACTCTGGTAACGCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTCGAAGCAGT‑TGCTCAGCTGGCGCAGGAAG‑TGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGAcggcga                                                                                                                                                                        <  2:285996/248‑2 (MQ=255)                                            ggcggcAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCgggggg                                                                                                                                                                    >  2:119311/1‑248 (MQ=255)                                              cggcAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTctgc                                                                                                                                                                  >  1:499323/1‑251 (MQ=255)                                               ggcAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTctgct                                                                                                                                                                 >  1:283346/1‑251 (MQ=255)                                                 cAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGtgggcttg                                                                                                                                                               >  2:312730/1‑244 (MQ=255)                                                     atCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGaa                                                                                                                                                           <  1:402086/251‑1 (MQ=255)                                                       cGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGt                                                                                                                                                         <  1:130415/251‑1 (MQ=255)                                                         aaaCTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGACTGCTGGGGAAGtac                                                                                                                                                       >  1:145136/1‑249 (MQ=255)                                                                 gCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGTGGTTAAACTCAGCTAACTCGATTACCATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTTACGTGGGGTCTGCTGGCGAAGTCACTCGGTCt                                                                                                                                               <  1:366352/251‑1 (MQ=255)                                                                   cATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGTGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCATCCAGCTGGTTGTCTCGGTGCCAAGCgg                                                                                                                                                                                               <  1:438621/201‑1 (MQ=255)                                                                   cATCGACGGCGATTTCGATGCCTGACAGGCGATGGTGAAACAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCCACCAGCTGGTTGTCTCGGTGCCAAGCgg                                                                                                                                                                                               >  2:438621/1‑201 (MQ=255)                                                                    aTCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGTGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTGTGAGGGCGAAGTGACTCGGttttgc                                                                                                                                            >  1:6450/1‑246 (MQ=255)                                                                          ggCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAACCTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGAGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTAGGTGCCAAGCGGCAACTTCGGCGATTTGACGGCGGGTGTGCTGGCgggaggaatcggtctggacgggg                                                                                                                                       >  1:73114/1‑228 (MQ=255)                                                                          ggCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCGGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGaa                                                                                                                                      <  2:132463/251‑1 (MQ=255)                                                                                            agggcgtgggagcaGCAGGCTATTGATGATGAAGAGCTGAAAGTGGTGCTGGGGTTAAACGCAGCTAACTCGATGAACATCGGCCGCTTGCTGTCGATGATGTGCTACTACGTTGAAGCAGG‑GGCGCAGTGGACGTAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTGCGGGGAGTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGAc                                                                                                                    <  2:172250/238‑1 (MQ=255)                                                                                                cgcTGGTGAAGCAGGCGGTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAAc                                                                                                                <  1:31195/251‑1 (MQ=255)                                                                                                                         ggaagagCTGAAAGTGGCGCTGGGTTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTG‑AACGATACCGTGCCACGTTTcc                                                                                      <  2:455692/250‑1 (MQ=255)                                                                                                                              ggcTGAAAGTGGCGCTGGGGTTAAACCCAGTTAACGCGATTAACATCGGCCGTTTGCGGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCGGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATGGCTGCGACCAACGTG‑AACGATACCGTGCCACGTTTCCTGCac                                                                                 <  2:435718/249‑1 (MQ=255)                                                                                                                                             gggggtttacctcagcttaCTCGATTAACGTTGGCCGTTTGCTGGCTCAGATTTGCTACTACTTTGAAGCGGT‑GGCGCAGCTGCCGCAGG‑GAGCGGGCAACCAGCTGTGTGTCTCGGTGCCAAGCGGAAACTGCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGCCTGCCGGTGAAACGTTTTATTGCGGCGGCCAACGTG‑AACTATACCGTGCCACGTTTACTGCACGACGGTCAGTGGTCa                                                                  <  2:474427/233‑1 (MQ=255)                                                                                                                                                  ttAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTTCTGGGGGAGGCACTCGGTCTGCCGGGGGAACGGTTTATTTCTGCGACCCACGTG‑GACGATACCGTGCCCCGTTTCCTGCACGACGGCCCGGGGTCAcccc                                                              >  1:332394/1‑249 (MQ=255)                                                                                                                                                     aaCTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCCGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCTCTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTG‑AACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGc                                                          <  1:206459/251‑1 (MQ=255)                                                                                                                                                               aCTCGATTAACATCAGCCGTTTGCTGGCGCTGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTG‑AACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCg                                                <  1:127871/251‑1 (MQ=255)                                                                                                                                                                   gtttAACAACAGCCGTTGGCTGGCGCAGATTTGTTACTACTTTGAAGC‑GTCTGCGCGTCGGCCGCAGGAAG‑CGCTCACCCAGCTGGTGGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGTCGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTGATTGCTGCGACCAACGTG‑AACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGa                                               <  2:295475/246‑1 (MQ=255)                                                                                                                                                                   gATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGGGAAGTCACTCGGTATGCCGGCGAAAAGTTTTTTTTATTCGAACAACGTGCAAC‑CTAACGTGCCACGGTTTCTGCCCGGCGGGCAATGGTCCACCAAAGCGAGTCAggagg                                               >  1:295475/1‑245 (MQ=255)                                                                                                                                                                                   ttGCTGGCGCAGATTTGCTACTACTTTGAATCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAGCCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTG‑AACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATgg                             <  2:447411/250‑1 (MQ=255)                                                                                                                                                                                        ggCGCAGATTTGCTACTACTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTGCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTG‑AACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGa                       <  1:130441/251‑1 (MQ=255)                                                                                                                                                                                               aTTTGCTACTTCTTGGATGCAGT‑TGCGCAGAGGCCGCAGGAAG‑CGCGCAACCAGCGGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTG‑AACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGcc                <  2:78451/251‑1 (MQ=255)                                                                                                                                                                                                    ctactaCTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAATCACTCGGTCTGCCGGGGGAACGGTTTATTGTTGCGACCAACGTG‑AAAGCTACCGGGCCCCGCTTCCTTCACGCCGGTCAGTGGTCACCCCAAACGACTCAGGCGACGTTCTCCAACGCGAGGGGCGTGAGTCAGCCGaac            >  1:278424/1‑250 (MQ=255)                                                                                                                                                                                                     tactacTTTGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGGAGTCACTCGGTCTGCCGGTGAAAAGGTTTTTTTCTGCGACCAACGGG‑AAAGATACCGTGCCACGTTTCCTGCACGACGGGCAGTGGTCACCCCAAACGACTCAGGCGACGATATCCAACGCGAGGGGCGTGGATCAGCCGagac           >  1:278361/1‑247 (MQ=255)                                                                                                                                                                                                             tGAAGCAGT‑TGCGCAGCTGCCGCAGGAAG‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCCCTCGGTCTGCCGGTGAAAAGGTTTTTTTCTGCGACCCACGTG‑AACGATACCGTGCCACGGTTCCTGCACGACCGTCAGTGGTCACCCAAAGCGACTCCGGCGGCGTTATCCAACGCGATGGACGGGGGGCAGCCCAACAACAGGCcgc  >  1:250062/1‑251 (MQ=255)                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        AGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGT‑TGCGCAGCTGCCGCAGGAGA‑CGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTG‑AACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGC  >  NC_000913/4210‑4664

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