Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	67,834	G→A	V168V (GTC→GTT)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	67,834	0	G	A	95.8%	148.8 / ‑2.4	48	V168V (GTC→GTT)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (29/17);  total (30/18)
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.56e-01

CGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTG  >  NC_000913/67589‑68080                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       cgattgtcgCCGCGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCAAACCGCCCAGGCATTGCCGGAGAAACCGAACTTGATCTGCGCGGCAACTTTATCACCATCGGGGACCGCCACTTCACGCATGTTATCGCCAAAACGGCCGACTTTCCGATGACGGGTCTCCTGTTTCGAAACCGc                                                                                                                                                                                                                                                   >  2:331528/9‑251 (MQ=255)            gTGGATTTGTGTGGCAGGCGTCTGGGTGTCGCAGTTTTCGTCCTCATCCACCAGCGCTTTAACACCGCCGGCGCTGATGGAGTCCCCCCCCTCCACCAGATCGCCAACCGCCCAGGTATTGCCGGAAAAACCGAACTTGCTCTGTGCGCCACCTTTACCACCATCGTTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATcc                                                                                                                                                                                                                                        <  2:259673/251‑1 (MQ=255)             tGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATcc                                                                                                                                                                                                                                        >  1:325181/1‑250 (MQ=255)                   tgggttgCAGGCTTCATGGTGTAGCAGCTTTCGTCCCCATCGCCCCGCGCTTTACCATCGCGTTCTGGGAGGGAGTCCACCACCGGCTCCAAAGGCCCACCCCCCCAGGTATGGACGGAGAAACCGATTTTGATCTGGGCGGCAACTTTACCACCATCGGGGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGc                                                                                                                                                                                                                                  <  1:238432/244‑1 (MQ=255)                      gttgCAGGCGTCGTGGTGTAGCATCTTCCGTACCCATCGCCCAGCGCTTTAACATCCCCTTCGCTGATGGATTTCCCCACCGGCCCCAGATCGCCAACCGCCCAGGTATGGAGGGAAAAACCGAACTTGTTCTGGGCGGCAACTTTATCCCCACCGGGGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGa                                                                                                                                                                                                                               <  1:220480/247‑1 (MQ=255)                         gCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACg                                                                                                                                                                                                                           >  1:222815/1‑251 (MQ=255)                               gTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATg                                                                                                                                                                                                                     >  2:79140/1‑251 (MQ=255)                                    ggTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCtt                                                                                                                                                                                                                >  1:44535/1‑251 (MQ=255)                                    ggTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCtt                                                                                                                                                                                                                >  1:175584/1‑251 (MQ=255)                                       gtAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCCCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGtt                                                                                                                                                                                                             >  2:14564/1‑251 (MQ=255)                                             gcTTTCGTAATCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCt                                                                                                                                                                                                       >  1:249840/1‑251 (MQ=255)                                                ttCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGcc                                                                                                                                                                                                    >  1:224380/1‑251 (MQ=255)                                                   gTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGCCGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGt                                                                                                                                                                                                 >  2:244009/1‑251 (MQ=255)                                                                                     cTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCTTAcg                                                                                                                                                               >  2:295283/1‑251 (MQ=255)                                                                                     cTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATAcg                                                                                                                                                               >  1:293938/1‑251 (MQ=255)                                                                                              ttCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAAt                                                                                                                                                      >  2:335058/1‑251 (MQ=255)                                                                                                  ccaccTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGa                                                                                                                                                    >  1:98392/1‑249 (MQ=255)                                                                                                          aCCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTg                                                                                                                                                                                                >  2:195056/1‑197 (MQ=255)                                                                                                          aCCAGACCGCCACCCGCCCAGGTATTGACGGAGAAACCGATCTTGATCTGGGCGGCACTTTTATCCCCAGGGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTg                                                                                                                                                                                                <  1:195056/197‑1 (MQ=255)                                                                                                                            cAGGTATTGACGGAAAAACCGAACTTAATCTGTGCGACAATTTTACCCCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGt                                                                                                                        <  1:228358/251‑1 (MQ=255)                                                                                                                               gTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTg                                                                                                                     >  2:92483/1‑251 (MQ=255)                                                                                                                                   tGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttc                                                                                                                 >  1:196914/1‑251 (MQ=255)                                                                                                                                    gagggaaaaaCCGACTTTGAGCTGGGCGGCAACTTTACCACCACCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttca                                                                                                                <  1:270628/244‑1 (MQ=255)                                                                                                                                        gagaAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGCGCCCGACCCTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGAAACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttcaggtt                                                                                                            >  1:367616/1‑251 (MQ=255)                                                                                                                                         agaAACCGAACTTGATCTGTACGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttcaggttc                                                                                                           >  2:366733/1‑251 (MQ=255)                                                                                                                                           aaaCCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCAt                                                                                                         >  1:210499/1‑251 (MQ=255)                                                                                                                                                     tGATCTGGGCGGCAACTTTATCCCCACGGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCat                                                                                                <  1:3728/250‑1 (MQ=255)                                                                                                                                                      gATCTGTGCGGCAACTTTACCCCCAGCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCata                                                                                               <  1:130338/250‑1 (MQ=255)                                                                                                                                                           gtgCGGCATCTTAACCGCCACCGGTGACCGCCTCTTCACGCATGTTATCCCCAAAGCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCg                                                                                            <  1:14564/248‑1 (MQ=255)                                                                                                                                                            tgCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGATATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAc                                                                                        >  2:185572/1‑251 (MQ=255)                                                                                                                                                                                 ggTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGt                                                                   >  1:126962/1‑251 (MQ=255)                                                                                                                                                                                    gcccGCCCCTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCAGGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGGCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGa                                                                <  2:161347/249‑1 (MQ=255)                                                                                                                                                                                     aCCGCCACTTCAAGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGCAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATACAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGa                                                                >  2:175311/1‑250 (MQ=255)                                                                                                                                                                                          cACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTggg                                                          <  1:335897/251‑1 (MQ=255)                                                                                                                                                                                             ttCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGgtgt                                                       >  1:156802/1‑251 (MQ=255)                                                                                                                                                                                              tCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGtgtg                                                      >  1:269250/1‑251 (MQ=255)                                                                                                                                                                                                        tATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTgcagc                                            >  1:372308/1‑251 (MQ=255)                                                                                                                                                                                                         agcGCCAAAACTGCAGACTTTCATATGACGGGTATCCTGGTTAGAAACCGCCGGACGCATCCAGGGGCCGATACGCTCAGGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGTCCGCCATGTGCGGTCTGGTTCAGGTTCATAAAGTCCAGATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTTTAATTgcagc                                            <  2:306699/248‑1 (MQ=255)                                                                                                                                                                                                         aTCGCCAAAACGGCAGACTTTCAGATGACGG‑TATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGGCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTgcagca                                           <  1:242256/250‑1 (MQ=255)                                                                                                                                                                                                           cGCCAAAACGGCAGTCTGTCAGATGACGGGGATCGTGTTTAGAAACCGCCGGTCGCATCAAGGAGCCGATACGCTCATGGGCTTTTTCATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGACCTGGGTGTGGAATTGCAGCAAc                                         <  2:110064/251‑1 (MQ=255)                                                                                                                                                                                                            gCCAAAACGGCAGACTTTCAGATGCCGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTGGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACg                                        <  1:66590/251‑1 (MQ=255)                                                                                                                                                                                                                  actgCAGCCTTTCAGCAGACGGGTATCCTGTTTAGAACCCGCCTGACCCATCACGGCGCAGATAAGCTCATGGGCTTGTTCATCCTGCCAGTGACCGGTACCGCCGGCATGTTGCTGCCGCATACGCGCGCCACTGAAGCCGACCTCGCCACCGCCCTGTGCAGTCTGGTTCAGGCTCATAAAGTCCATATCGATACTGTCCCACGTCAGCGCCGCGTTGGCCTGGGTGTGGACTTGCAGCAATGGTttgt                                  >  2:363726/4‑251 (MQ=255)                                                                                                                                                                                                                       aGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGCTACGCTCATGGGCTTGTTTATCCTGCCAGTGCCCGGTAACGACGGCATGTTGCTGACGCATCCGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCCTATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTTCTGCAACGGTTTGTTGAGc                             >  2:223024/1‑251 (MQ=255)                                                                                                                                                                                                                               aGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTTAGCATGGTCAg                     <  1:83978/251‑1 (MQ=255)                                                                                                                                                                                                                               aGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAg                     <  1:188849/251‑1 (MQ=255)                                                                                                                                                                                                                                  ggACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGcc                  <  2:123889/250‑1 (MQ=255)                                                                                                                                                                                                                                                 ttAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACAtttt   >  2:288013/1‑251 (MQ=255)                                                                                                                                                                                                                                                  tAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTg  >  1:241195/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       CGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTG  >  NC_000913/67589‑68080

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