Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,937,909	G→A	A206T (GCG→ACG)	rbsK →	ribokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,937,909	0	G	A	98.0%	169.1 / NA	49	A206T (GCG→ACG)	rbsK	ribokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/25);  minor base C (1/0);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGGTGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAAC  >  NC_000913/3937665‑3938153                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     ggtggAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTAGTCGTGTTGAAAATGATGAAGATGCAACGaa                                                                                                                                                                                                                                                  <  1:293163/249‑1 (MQ=255)         cgcAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAg                                                                                                                                                                                                                                        <  1:48399/251‑1 (MQ=255)               gTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTg                                                                                                                                                                                                                                  >  1:169975/1‑251 (MQ=255)                            cAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTa                                                                                                                                                                                                                     <  1:252242/251‑1 (MQ=255)                                    cAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCgtact                                                                                                                                                                                                             <  2:303215/251‑1 (MQ=255)                                     aGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCgtactg                                                                                                                                                                                                            >  1:9043/1‑251 (MQ=255)                                              aaTGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGAtt                                                                                                                                                                                                    <  1:64044/250‑1 (MQ=255)                                                tGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATAGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGATGGTGGACATTATTACGCCCAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGAGCAGGTACTGCATGAAAAAGGAATCCGTACTGTACTGATTACt                                                                                                                                                                                                 >  1:173296/1‑251 (MQ=255)                                                      tGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTTCTTTAgg                                                                                                                                                                                            >  2:298261/1‑250 (MQ=255)                                                       ggAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGTAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCCTGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGa                                                                                                                                                                                           >  1:24876/1‑250 (MQ=255)                                                                    gAAAGTGTGATGGCAGCGGCGAACATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCCGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGGGTATg                                                                                                                                                                             >  1:98354/1‑251 (MQ=255)                                                                                  aGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATg                                                                                                                                                               >  1:175981/1‑251 (MQ=255)                                                                                    cggcgAAAATCGCCCGTCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATgat                                                                                                                                                             <  1:365074/251‑3 (MQ=255)                                                                                       cgAAAATCGCCCATCAAAATAAGACTACCGTTGCGCTTAACCCGGCTCCGGCTCGCTAACTTCCTGACGAACTTCTGGCTCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTTTTGAAAATGATGAAGATGCTACGAAGGCGTCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTTATTACTTTAGTAAGTCTTGGTGTATTGGCTAGCGTGAATGATGaa                                                                                                                                                          <  2:316518/251‑1 (MQ=255)                                                                                         aaaaTCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAgg                                                                                                                                                        <  1:393958/251‑1 (MQ=255)                                                                                                        aaTAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTgg                                                                                                                                         >  1:255391/1‑251 (MQ=255)                                                                                                         aTAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGaa                                                                                                                                                          <  1:288431/233‑1 (MQ=255)                                                                                                         aTAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGAATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGaaagatcggaagggcgtcgg                                                                                                                                        >  2:288431/1‑233 (MQ=255)                                                                                                                   cGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGGCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGc                                                                                                                              <  2:153817/251‑1 (MQ=255)                                                                                                                        cgcTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCt                                                                                                                         <  2:77493/251‑1 (MQ=255)                                                                                                                             aaCCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGCAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGAGCAGGCTGTCg                                                                                                                     >  1:357084/1‑250 (MQ=255)                                                                                                                                ccggctccggctCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATGCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGAGTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATAc                                                                                                                 <  1:38060/251‑1 (MQ=255)                                                                                                                                ccggcgccggctCTCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTTAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGAttc                                                                                                                 <  2:175981/251‑3 (MQ=255)                                                                                                                                            cgcgAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGAGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCgg                                                                                                     >  2:336018/1‑251 (MQ=255)                                                                                                                                              cgAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGgag                                                                                                   <  1:208903/251‑1 (MQ=255)                                                                                                                                                  cTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATAc                                                                                               >  2:274999/1‑251 (MQ=255)                                                                                                                                                         aCGAACTGCGAGCGCTGGTGGTCATTATGACGCCAAACTAAAAGGAAGCAGAAAAGCGCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGATGGCGGCGCTGGGACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGACTCCGGGTGCAGGCTGTCGAGACCATTGCTGCCGgagttcccattaa                                                                                         <  2:105867/250‑11 (MQ=255)                                                                                                                                                           gAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGAAGCAAAGCTCACCGGTATTCGTGTTGACAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACgg                                                                                      >  1:267261/1‑251 (MQ=255)                                                                                                                                                                   ggCGCTGGTGGACATTATTACGCCAAACGACACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTa                                                                               >  2:84627/1‑250 (MQ=255)                                                                                                                                                                   ggCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTa                                                                               >  1:234661/1‑250 (MQ=255)                                                                                                                                                                     cgcTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATc                                                                            <  1:392258/251‑1 (MQ=255)                                                                                                                                                                          gtggACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCTATACCATTGCTGCCGGAGATTCCTTTAACGGTGCGTTAATCACGGc                                                                       <  2:73181/251‑1 (MQ=255)                                                                                                                                                                           tggACATTATTACGCCACACGAAACGGAAGCAGAAAAGCTCACCGGCATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGAACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGcc                                                                      >  2:308588/1‑250 (MQ=255)                                                                                                                                                                            ggACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCa                                                                      <  2:124088/250‑1 (MQ=255)                                                                                                                                                                             gACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCAtt                                                                    <  2:131779/251‑1 (MQ=255)                                                                                                                                                                                 ttattaCGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCt                                                                                                                                                                         >  1:217898/1‑146 (MQ=255)                                                                                                                                                                                 ttattaCGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCt                                                                                                                                                                         <  2:217898/146‑1 (MQ=255)                                                                                                                                                                                               gaaacgGAAGCCGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTCTCCGTACTGTACTGATTACTTTAGGCAGTCGTGGAGTATGGGCTAGCGTGAATGATGACGGTAAGCGCGTTCCTGGCTTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAACCCAtt                                                  >  2:270682/1‑251 (MQ=255)                                                                                                                                                                                                    ggAAGAAGACAAGCTCACCGGTATTCGTGTTGAATATGATGAAGATGCAACGACGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTAATGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGACTGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGcc                                               >  2:260440/1‑249 (MQ=255)                                                                                                                                                                                                     gAAGCAGAAAAGCTCCCCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCAt                                                   >  1:303473/1‑244 (MQ=255)                                                                                                                                                                                                     gAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCAt                                                   <  2:303473/244‑1 (MQ=255)                                                                                                                                                                                                      aaGCAGGAAAGCTCTCCGTGATTCGTGTTGAAAATGATGAAGAGGCAACGAATGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTGCTGTACTGATTACTTTAGGAAGTCTGGGTGGATGGGCGAGCTTTAATGATGAAGGTCAGCGCGTTCCTGGATTCCAGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTATAACGGTGCGTTAATCACGGCAGTGCTGGAAGAAAAACCATTTCCagag                                           <  2:128066/251‑1 (MQ=255)                                                                                                                                                                                                            aaaaACTCAACGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGAGGTGTATGGACTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGCAGAAAAACCATGGCCAGAGACGAtt                                     >  2:302092/1‑251 (MQ=255)                                                                                                                                                                                                                 cTCACCGGTATTCGTGTTGAAAATGATGCAGATGCAACGAAGACGGCGCAGGTACTGAATGAAAAAGGAATCCGTACTGTACTGATTACTTTAGGAAGTAGTGGTGTATGGGCTAGCGGGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGACGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGACAGAGGCGATTCGttt                                >  1:105187/1‑251 (MQ=255)                                                                                                                                                                                                                    cccGGTATTCGTGTTGAAAATGATGAAGAGGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCGGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGc                             <  2:274678/250‑1 (MQ=255)                                                                                                                                                                                                                                  tGAAAATGATGAAGATGCAACGAAGGCGGCGCAGTTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTTAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCAGTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCt                <  2:195269/250‑1 (MQ=255)                                                                                                                                                                                                                                            gaAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGCGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGt     >  2:184306/1‑251 (MQ=255)                                                                                                                                                                                                                                              aGATGCACCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTa    >  1:312591/1‑250 (MQ=255)                                                                                                                                                                                                                                               gATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAac  <  1:394081/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     GGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGGTGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAAC  >  NC_000913/3937665‑3938153

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