Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	68,678	T→C	A457A (GCA→GCG)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	68,678	0	T	C	97.3%	114.6 / NA	37	A457A (GCA→GCG)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (15/21);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGCTCGCTGCACGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGAGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑TGCCATCACGTTATTAACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTC  >  NC_000913/68431‑68911                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          gATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGcc                                                                                                                                                                                                                                        >  2:212050/1‑251 (MQ=255)      cggcgATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCAc                                                                                                                                                                                                                                   <  1:113749/251‑1 (MQ=255)             aaaGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCCTCACGttattg                                                                                                                                                                                                                            >  2:8260/1‑250 (MQ=255)               aGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGttattgac                                                                                                                                                                                                                          >  2:128497/1‑248 (MQ=255)                   gTTCAAAGCGTTGTGCCTTTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCg                                                                                                                                                                                                                      <  2:173137/251‑1 (MQ=255)                                     gttCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAg                                                                                                                                                                                                    <  1:212050/249‑1 (MQ=255)                                             cgcGGTTGCAGGGTTTTCTCTACCGCACTGGCCATGTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATGGACGGCGATCCCCTGATCGTTAAAGCACTCCAt                                                                                                                                                                                            <  2:290793/249‑1 (MQ=255)                                               cGGTTGCAGGGTTTTCTCTCCCGCACTGGCCATTTTTTGCTGGGCTGAGGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTGTCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTGTTTCCGCGCGATGCCTCCCAGCTGC‑ATCACGTTATTGTCTGCGGTCGCCTGATCGGTAAAGCATTCCATGa                                                                                                                                                                                          <  2:25979/251‑1 (MQ=255)                                               cGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATg                                                                                                                                                                                           >  1:199847/1‑250 (MQ=255)                                                       gggTTTTCTCTACCGCACGGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGGCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCATGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGt                                                                                                                                                                                  <  2:283656/251‑1 (MQ=255)                                                        ggTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCTCTCCATGATTGCGCGTg                                                                                                                                                                                 <  2:181253/251‑1 (MQ=255)                                                                   ccGCCCTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAACACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATTCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCg                                                                                                                                                                      <  2:215507/251‑1 (MQ=255)                                                                       aCTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCAGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTgg                                                                                                                                                                  >  2:295528/1‑251 (MQ=255)                                                                            ccATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGgcagc                                                                                                                                                              >  1:265120/1‑250 (MQ=255)                                                                                   ttGCTGGGCTGATGGGATGTCTGCGTGCACTTGCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCaaa                                                                                                                                                       <  2:56751/250‑1 (MQ=255)                                                                                        gggctgaTGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAACTCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCACCAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAAccgccg                                                                                                                                                 >  1:316702/3‑251 (MQ=255)                                                                                         ggctgaTGGGATGTCTGCGTGCCCTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCCCTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGGCCTGGTTTTTCCGCGCGATGCCCCCCAG‑CGCCATCACTTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGa                                                                                                                                                <  1:139084/250‑1 (MQ=255)                                                                                                  gATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGACGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCCATCAAACCGCCGCACAGCAGCg                                                                                                                                       >  1:126075/1‑251 (MQ=255)                                                                                                     gTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGa                                                                                                                                     <  1:263938/250‑1 (MQ=255)                                                                                                                               cagcaAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCCATCCAACCGCCGAACAGCCGCGGAGCGTCGGTAGCGAGGGTAAGATCGGt                                                                                                           >  2:235316/1‑250 (MQ=255)                                                                                                                                   aaaaaTCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCa                                                                                                      <  1:146946/251‑1 (MQ=255)                                                                                                                                   aaaaaTCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCa                                                                                                      <  1:222985/251‑1 (MQ=255)                                                                                                                                   aaaaaTCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTGCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCa                                                                                                      <  2:199611/251‑1 (MQ=255)                                                                                                                                    aaaaTCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCAc                                                                                                     <  1:155160/251‑1 (MQ=255)                                                                                                                                             gcACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAg                                                                                            <  1:268340/251‑1 (MQ=255)                                                                                                                                                ccGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCg                                                                                         <  2:218500/251‑1 (MQ=255)                                                                                                                                                  gAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGGACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGtt                                                                                       <  1:150085/251‑1 (MQ=255)                                                                                                                                                                   cAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGtttggcg                                                                        >  1:83145/1‑244 (MQ=255)                                                                                                                                                                       ggCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGtttggcgagcgac                                                                  <  1:33465/251‑12 (MQ=255)                                                                                                                                                                                    aGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATCATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGGTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAAcc                                                      >  2:25868/1‑250 (MQ=255)                                                                                                                                                                                       gggCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAAcc                                                                                                                                                     <  2:53640/152‑1 (MQ=255)                                                                                                                                                                                       gggCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAAcc                                                                                                                                                     >  1:53640/1‑152 (MQ=255)                                                                                                                                                                                                           gcagGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGAACGTTAAACCAGTCGAGCACCACCGGCAGGTg                                >  2:162309/1‑249 (MQ=255)                                                                                                                                                                                                             agGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCCTGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGCACAGCAGCGGAGCGGCGGTTGCGAGGGTAAGATCGGTAATCACCCCTTTAAGGCGTTGGTTAGCGATTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATcc                            >  2:284529/1‑251 (MQ=255)                                                                                                                                                                                                                                  ttttCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGAAGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGAGATCCAGAGACGGATTTTTGGCCCAt       >  1:309705/1‑251 (MQ=255)                                                                                                                                                                                                                                   tttCCGCGCGATGCCGCCCAG‑CGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCAt       <  2:198046/250‑1 (MQ=255)                                                                                                                                                                                                                                       cgcgcgATGCCGCCCAG‑CGCCATCACGTTATGGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTc  <  2:299943/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          GATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGCTCGCTGCACGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGAGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAG‑TGCCATCACGTTATTAACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTC  >  NC_000913/68431‑68911

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