Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,069,743	C→T	A33T (GCA→ACA)	rutG ←	pyrimidine:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,069,743	0	C	T	97.8%	132.2 / NA	46	A33T (GCA→ACA)	rutG	pyrimidine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (31/13);  minor base G (0/1);  total (31/15)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGCCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGG  >  NC_000913/1069513‑1069989                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        gATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATc                                                                                                                                                                                                                                    >  1:239994/1‑251 (MQ=255)         ggtaaagccccTGACCGTAAAACCCGGTCGCGGCAATCCCCACGCCGAAAAAAGCCGCGCGGGAGCCGAGATAACGGGGGAGGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATcc                                                                                                                                                                                                                            <  1:145902/243‑1 (MQ=255)               tGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGGCGCGCCCGCCGGTAATAAAGAACAACAGTAACGTGCCGATCCCCGGCATTACAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGAATTAACACCGTCGCGCCAAACATCGCCACCGCGGGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGgcggg                                                                                                                                                                                                                      >  2:405885/1‑249 (MQ=255)                                                                               aTAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATc                                                                                                                                                       >  1:155224/1‑250 (MQ=255)                                                                                     ggggcggcgcCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGtc                                                                                                                                                <  1:170535/245‑1 (MQ=255)                                                                                     ggggACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGCCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATAAGCTGtc                                                                                                                                                >  1:43824/1‑251 (MQ=255)                                                                                        gccgcgCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTcc                                                                                                                                             <  1:144040/249‑1 (MQ=255)                                                                                        gACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTcc                                                                                                                                             >  1:267782/1‑251 (MQ=255)                                                                                            cgcccgccGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTa                                                                                                                                          >  1:424385/1‑250 (MQ=255)                                                                                                 gccgGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAaggagg                                                                                                                                    >  1:196031/1‑251 (MQ=255)                                                                                                   cggaaataaAGAAAAACAGTAACGTGCCAACCCCGGACATTAAAATGGAAAGATTGGGATCCGGCCCCATCAGAATCGGCATTAACGCCGTCGCGCCAAACATCGCTACCGCGTGT‑TGAACCCCCAGTACTGTCGTCTGAGCAAACGG‑CAATCGTTCAGCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGtt                                                                                                                                  <  1:312102/247‑1 (MQ=255)                                                                                                            aGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACAACGCCGCTCTCTGTAGAGGTCGCTTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGc                                                                                                                           >  1:301759/1‑249 (MQ=255)                                                                                                              aaaaacattcccgttccGATCCCCGACATTCAAATGGAACGATTGGGATCCAGCCCCATCAGCATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGCGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAggg                                                                                                                       >  1:392790/16‑251 (MQ=255)                                                                                                                      tAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATc                                                                                                                 >  1:499067/1‑249 (MQ=255)                                                                                                                          gTGCCGATCCCCGAAATTACAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCAACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAACCCGAACATTGCCATCCGCTGTCTCCTTAAGGAGGTTATCAAGCAGGGCGCATCAGcgcg                                                                                                           >  2:213632/1‑251 (MQ=255)                                                                                                                              cGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CGATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGAt                                                                                                       >  1:149020/1‑251 (MQ=255)                                                                                                                               gATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCAAAACATCGCCACCGCGTGT‑TGACCGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGCCGCGACCACGCCGCTCTATGTCGAGGACGATTTTAACTGCCAGTGCGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGAt                                                                                                       >  2:182573/1‑250 (MQ=255)                                                                                                                                     cGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCAATAAAACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACAGAACATTGAAATAAGCTGTCTCCTTAAGGAGGTTAAAAAGCAGGGCGCAGCAGCGCGAGATAACTGCg                                                                                                >  1:311625/1‑251 (MQ=255)                                                                                                                                         aTTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCg                                                                                            >  1:50611/1‑251 (MQ=255)                                                                                                                                             aaaTGGAAAGATGGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCCCCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGCCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAAcc                                                                                        <  1:275121/251‑1 (MQ=255)                                                                                                                                              aaTGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAAcc                                                                                        >  2:178819/1‑250 (MQ=255)                                                                                                                                                               ccAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCCCCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGt                                                                            <  1:207286/245‑1 (MQ=255)                                                                                                                                                               ccAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtg                                                                      >  1:175375/1‑251 (MQ=255)                                                                                                                                                               ccAGCCACATCAGAATCGGCATTCACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCCAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGt                                                                            >  2:207286/1‑245 (MQ=255)                                                                                                                                                                 aGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtg                                                                    >  2:160003/1‑251 (MQ=255)                                                                                                                                                                   ccccATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtggt                                                                  <  1:178819/251‑1 (MQ=255)                                                                                                                                                                   ccccATCAAAATCGGCATTACCACCGTCGCTCCAAAACTCGCCACCGCGTGTGTGA‑CGCCCATGACTGTCGTCGTAGCAAAAGT‑CAATCGTGCAGCCGGCGCTACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtggt                                                                  <  2:422467/251‑1 (MQ=255)                                                                                                                                                                    cccATCAGAACGGGCATTAACACCGTCGCGCCAAATAGCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGGGCAAACGG‑CAATCGTTCATCCGGCGCGTCCACGCCGCTCTCTGTAGAGGTCGAGTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtggtg                                                                 <  1:475629/251‑1 (MQ=255)                                                                                                                                                                     ccATCAGAATCGGCGTTAACACAGTCGCGCCAAACATCGCCCCCGCGTGT‑TGACCGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCCGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtggtgt                                                                >  2:222565/1‑251 (MQ=255)                                                                                                                                                                      cATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCTCGTGT‑GGACCGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGgtgtggtgtg                                                               <  1:469465/251‑1 (MQ=255)                                                                                                                                                                      cATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGAGATCACTGCGATCGAACCACACCAGCCCGTAGGgtgtggtgtg                                                               >  2:247628/1‑251 (MQ=255)                                                                                                                                                                       aTCAGAAGCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGa                                                              <  1:307192/251‑1 (MQ=255)                                                                                                                                                                        tCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGGGTGGTGTGAc                                                             >  1:411855/1‑251 (MQ=255)                                                                                                                                                                                  cATTAACACCGTAGCGCCAAACATCGCCACAGCGTGT‑TGACCGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCCGCGCGTGATAACTGAGATCGAAACACACCAGCCCGTAGGGTGCGGTGTGACGATGAATCGc                                                   >  2:351191/1‑251 (MQ=255)                                                                                                                                                                                    ttAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAgg                                                                                                                                                                                                <  1:481359/108‑1 (MQ=255)                                                                                                                                                                                    ttAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCACATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAgg                                                                                                                                                                                                >  2:481359/1‑108 (MQ=255)                                                                                                                                                                                     tAACAACGTCGCGCCAAACATGGCCTCCGCGTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACAACGCCGCTCTCTGTAGAGGTAGATTTTACGTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGCTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGAACGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTc                                                >  2:131376/1‑251 (MQ=255)                                                                                                                                                                                            gccccgcCAAACAACGCCACCGCGGAT‑TGAACGCCCATTACGTCCGTCTGAGCAAACGG‑CAATCTTTCACCCGGCGCGCCCACGCCGCTCTCTGGAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGTTGTGGTGTGACGATGAATCGCTTCGATGgcg                                         <  2:174309/247‑1 (MQ=255)                                                                                                                                                                                                  ccaaacataacgaccgagggt‑TGAAAGCCCATGACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGGGACCACGCAGCACTCTGTAGAGGTCGATATTAAATGACCGGGAGGAAAACCGAACATTGCCATCAGATGACTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGGGATAACTGCGAGCGAACCAAACCAGCCCGGAGGGGGCGGTGTGCCGAGGAGTCGCTTCGATAGAGCaaaa                                    >  2:248584/20‑250 (MQ=255)                                                                                                                                                                                                           gCCACCACGTGT‑TGACCGCCCATTCCTGTCGCCTGAGCAAACGG‑AAATCGTTCATCCGGCGCGACCACGCCACTCTATGTAGAGAACGATTTTAACTGCCAGTGCGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGGGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTCGGCTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAAAAGacagtcg                          >  2:47402/1‑245 (MQ=255)                                                                                                                                                                                                                 gcgTGT‑TGAACGCCCATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGAGTGGTGTGACGCTGAATCGCTTCGATGGCGCAAAACGGAATGTCGTggggg                    >  2:368670/1‑249 (MQ=255)                                                                                                                                                                                                                   gtgt‑TGAACGCCCATTACTTTCGGCGGAGCAAACGGTCAA‑CGTTCACCCGGCGCGCCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCCTCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGAGCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGcc                  <  2:325429/251‑1 (MQ=255)                                                                                                                                                                                                                           aCGCCCATTACTGGCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTc           <  1:349865/251‑1 (MQ=255)                                                                                                                                                                                                                               ccATTACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTcacc        >  1:373847/1‑250 (MQ=255)                                                                                                                                                                                                                                   tACTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTg   >  1:280029/1‑251 (MQ=255)                                                                                                                                                                                                                                    aCTGTCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTgg  >  1:117284/1‑251 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        GATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGT‑TGAACGCCCATTACTGCCGTCTGAGCAAACGG‑CAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGG  >  NC_000913/1069513‑1069989

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