Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,306,960	G→A	A280A (GCC→GCT)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,306,960	0	G	A	95.3%	188.9 / 0.1	64	A280A (GCC→GCT)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (32/29);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

ATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTT  >  NC_000913/4306728‑4307204                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     attaGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGAcc                                                                                                                                                                                                                                    >  1:162430/1‑251 (MQ=255)        aTTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTATG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAATCGCACGACCTGATGCg                                                                                                                                                                                                                             <  1:354591/251‑1 (MQ=255)                     agGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGGG‑CCAACGTTGTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGAGGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGgcg                                                                                                                                                                                                                <  2:156668/251‑1 (MQ=255)                           aaTTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGCCCTGATGCGGCAGTGAACGGCGCAgga                                                                                                                                                                                                           >  1:199925/1‑249 (MQ=255)                            atgaCAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTAc                                                                                                                                                                                                         <  1:431404/248‑1 (MQ=255)                                                  aGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACgag                                                                                                                                                                                     >  1:327814/1‑249 (MQ=255)                                                   gAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGagagg                                                                                                                                                                                  >  2:306210/1‑250 (MQ=255)                                                       ataataTGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTcg                                                                                                                                                                              <  1:423302/251‑1 (MQ=255)                                                            atGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGGCATGGCAACGAGAGCCTCGCGTgg                                                                                                                                                                         >  1:57416/1‑251 (MQ=255)                                                            atGAGTCGTTTTATGTGTTTATTATATATATTGGACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGGCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGGCGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTgg                                                                                                                                                                         <  1:428551/251‑1 (MQ=255)                                                             tGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTggg                                                                                                                                                                        >  1:467118/1‑251 (MQ=255)                                                                gTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGa                                                                                                                                                                       >  1:11182/1‑249 (MQ=255)                                                                      ttATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCACATATTGCTGCACCCTGAGAGCAATTAAAGTCAGAAGATAAGACAGCAATAAAGCGCACGACCTGATGCGGCGGTGGACGGGGGAGGGACTTTTGGGTCAAGGCAACgagggctacgcgtgggaggtggggg                                                                                                                                                               >  2:400992/1‑229 (MQ=255)                                                                       tATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGc                                                                                                                                                               <  1:101650/250‑1 (MQ=255)                                                                          gtgtTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCata                                                                                                                                                            <  1:154742/250‑1 (MQ=255)                                                                                    tatatatTGTACCTGACGTTAGCACGGCGGGCCATTGCGCTATTGAA‑ACAACGTTTTACTCATGGGGCTTTAGCACAGACCTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCCCCCTGAGCGCCATTCAAGTCAGAAGATGAAGCAGCTATAAAGCGCACGACCTGATGCGGACGTGGACGGAGCAGGAACTTTTGGGCAATGGAAACGAAAGACTCGCGTGGGCAGGAGGGCatatacgactcgtc                                                                                                                                                 >  2:387487/1‑241 (MQ=255)                                                                                     atatatTGTACCTGGCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGcc                                                                                                                                                 <  1:33034/250‑1 (MQ=255)                                                                                       atatTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCAcc                                                                                                                                              >  1:386768/1‑251 (MQ=255)                                                                                        tatTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAACGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGCGTCTCGCGTGGGAAGGCGGGGTAATCCATCCCGCCACCg                                                                                                                                             >  2:221334/1‑251 (MQ=255)                                                                                        tatTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTAATGGGGCTTTAGCACAGAACAGTGGCAAAAAAAGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAACGTCAGAAGATGAGGCAGCAATAAAGCGCACGCCCTGATGCGGCCGTGGACGGCGCAGGTACTTTTGGATCATGGCCAAGAGAGTCTCGCGTGGGAAGGCGGGCATATCCCACACGCCAccc                                                                                                                                             >  1:467479/1‑250 (MQ=255)                                                                                            gTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGAAAGGCGGGCATATCCATAAAGCCACCGccc                                                                                                                                          >  1:235225/1‑250 (MQ=255)                                                                                             tACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGccc                                                                                                                                          >  1:317088/1‑249 (MQ=255)                                                                                                cTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGCCTCGCGTGGGAAGGAGGGCATATCCATCACGCCACCGCCCAGGAt                                                                                                                                     >  2:315245/1‑251 (MQ=255)                                                                                                   gcGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCAc                                                                                                                                  <  2:370166/250‑1 (MQ=255)                                                                                                      ttGGCACAGCGGTACATTGCGCCAC‑GTGGCCAACGGTTTACTCATGGGTCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAAGATTGCTGCTCCCGGAGCGCCATTAAATGGAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTGTGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCTCGCCACCGCCCAGTATCAcc                                                                                                                                 <  2:467242/249‑1 (MQ=255)                                                                                                      ttGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGCTGCGGCAGt                                                                                                                                                                                                                        >  2:431901/1‑162 (MQ=255)                                                                                                      ttGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGt                                                                                                                                                                                                                        <  1:431901/162‑1 (MQ=255)                                                                                                        ggCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCAt                                                                                                                             <  1:306210/251‑1 (MQ=255)                                                                                                         gCACAGCGGGGCATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATc                                                                                                                            >  1:416941/1‑251 (MQ=255)                                                                                                             aGCGTGACATTGCGCGACTGTG‑TCAACGTTTTACTCATGGGGTTTTAGCACAGAACTGTGGCAAAAAACGGTGATGCGCCTATATTGCTGCACCCTGATCGCCATTAAAGGCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCATGTACTTTGGGGTCATGGCAGCGAGAGTCTCGCGCGGGAAGGCTTGCATATCCAGGACGCCACCGCACAGGATCACCGCATCgggt                                                                                                                        <  2:56866/251‑2 (MQ=255)                                                                                                                  gACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTAATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGa                                                                                                                    >  1:57514/1‑250 (MQ=255)                                                                                                                   aCATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAAc                                                                                                                  <  2:467118/251‑1 (MQ=255)                                                                                                                    cATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCAcc                                                                                                                                 >  1:243284/1‑235 (MQ=255)                                                                                                                    cATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCAcc                                                                                                                                 <  2:243284/235‑1 (MQ=255)                                                                                                                      ttGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAAc                                                                                                                  >  1:326975/1‑248 (MQ=255)                                                                                                                             aCTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATg                                                                                                         <  2:317088/250‑1 (MQ=255)                                                                                                                             aCTGTG‑CCAACGTTTTACTCATGGGACTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAAt                                                                                                          <  2:327814/249‑1 (MQ=255)                                                                                                                                    ccAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCtggtgg                                                                                                  <  1:119192/251‑1 (MQ=255)                                                                                                                                             tACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGcc                                                                                          <  1:202377/250‑1 (MQ=255)                                                                                                                                               cTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGt                                                                                       >  1:404852/1‑251 (MQ=255)                                                                                                                                                tCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCCATGGCCCGTg                                                                                      >  2:359353/1‑251 (MQ=255)                                                                                                                                                    ggggCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAg                                                                                                                                        <  1:312823/197‑1 (MQ=255)                                                                                                                                                        cTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGAAACGAAGGTCTCGCGGGGGAAGGCGGGCATATCCATCACGCCCCCGCCCAGGATGACCGCATCGGGTGCGACAAGTTTCATGCGGGTGGCAATGGcacggtcccggttt                                                                              >  2:12079/1‑238 (MQ=255)                                                                                                                                                          ttAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTAGACGGCGCAGGTACTTTTAGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCAAACAGGTTAATGCAGGAGGCAATGGCCCggtacgctggttt                                                                            >  2:320118/1‑239 (MQ=255)                                                                                                                                                             gCACAGAACGGTGGCAAAAAACGGTGATGAGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCCGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCAGCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAg                                                                         <  2:114245/251‑1 (MQ=255)                                                                                                                                                                                aaCGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAg                                                        >  1:377925/1‑249 (MQ=255)                                                                                                                                                                                  cGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATACAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGTATCGAGCAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGGTTTCAAGCAGACTCTGGACGGAAGAGGc                                                    >  2:166870/1‑251 (MQ=255)                                                                                                                                                                                   gTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCg                                                   >  1:196924/1‑251 (MQ=255)                                                                                                                                                                                      gatacgcCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATAGGGATCGAACAGATTAGTGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAATCAGACTCTGGACGAAAGAGGCGttt                                                >  1:481542/5‑251 (MQ=255)                                                                                                                                                                                                tatTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATgg                                       >  1:200806/1‑250 (MQ=255)                                                                                                                                                                                                 atTGCTGCACCCTGAGCGCCATGAAAGTCAGAAGATGATGCAGCAATAAAGCGCTCGACCTGATGCGGTAGTGGACGGCGCAGGTTCTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATTGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACTATAGATGCGTTTTCTGCATGGa                                      <  2:410051/250‑1 (MQ=255)                                                                                                                                                                                                      ggcACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCGGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGaaaa                                <  2:196924/250‑1 (MQ=255)                                                                                                                                                                                                       gcACCCTGAGCGCCATTAAAGGCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGTCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATTCTGGTGGCAATGTCCCGTGCCGCTTTTTCAAGCAGTCGCTGGACGAAAGAGGCGTTTGCCTCATGGACGaaat                                <  2:240642/250‑2 (MQ=255)                                                                                                                                                                                                        cACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCAAGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTACTGCTGGTGGCAATGGCCCGAGCCGCGTTTTCAAGCAGACTCTGGGCGAAAGAGGCGTTTTCCGCATGGACgcaa                                 >  2:247397/1‑245 (MQ=255)                                                                                                                                                                                                         aCCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCACTACAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACCGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCCTGGCCGAAAAGa                              >  2:410608/1‑250 (MQ=255)                                                                                                                                                                                                                 cgcCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAAt                     <  1:247397/251‑1 (MQ=255)                                                                                                                                                                                                                  cccATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGTAGTGGACGGCGCCGGTACTTTTGGGTCATGGCAACGAGAGTCCCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCATTGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATg                    <  2:416941/250‑1 (MQ=255)                                                                                                                                                                                                                   ccATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACGGTTGGGTCATGGCAACGAGAGGCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAAGGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGATTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATgg                   <  1:170458/251‑1 (MQ=255)                                                                                                                                                                                                                     aTTAAAGTCAGAAGATGAGGCAGCAATGAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGt                 <  2:326975/251‑1 (MQ=255)                                                                                                                                                                                                                       tAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGt                                                                                       <  2:45378/179‑1 (MQ=255)                                                                                                                                                                                                                       tAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGt                                                                                       >  1:45378/1‑179 (MQ=255)                                                                                                                                                                                                                              agaagaTGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTTGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCggg        >  1:182598/1‑251 (MQ=255)                                                                                                                                                                                                                                   aTGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCATGTCCTTTTGGGTCATGGCAACGAGAGTCGCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAAGGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACTAAAAGATCGCGCAATGGGTAATTTCGGGGCtgt   <  1:341984/251‑1 (MQ=255)                                                                                                                                                                                                                                    ggAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCtgtt  <  2:404852/250‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     ATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTG‑CCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTT  >  NC_000913/4306728‑4307204

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