Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,307,038	G→A	5.3%	F254F (TTC→TTT)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,307,038	0	G	A	5.3%	241.7 / 2.3	75	F254F (TTC→TTT)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand):  ref base G (35/36);  new base A (2/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

TGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCA  >  NC_000913/4306894‑4307183                                                                                                                                                 |                                                                                                                                                  tgGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGa                                                                                                                                                  <  6:348503/146‑1 (MQ=255)   gCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAgg                                                                                                                                               <  8:342836/147‑1 (MQ=255)   gCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGc                                                                                                                                              <  5:295117/148‑1 (MQ=255)   gCAAAAAACGTTGACGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCTCACGACCTGATGCGGCAGGGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGc                                                                                                                                              <  6:342569/148‑1 (MQ=255)    cAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCg                                                                                                                                             >  3:354953/1‑148 (MQ=255)    cAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCg                                                                                                                                             <  4:354953/148‑1 (MQ=255)       aaaaCGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCggg                                                                                                                                           >  3:333417/1‑147 (MQ=255)       aaaaCGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGc                                                                                                                                          <  6:361362/148‑1 (MQ=255)                cagcggCAATATTGCTGGCCCCTTCGCCCCAATAAAGTCAGAACTTGCCGCCGCAATTAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATc                                                                                                                                 <  4:268373/146‑1 (MQ=255)                aTGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATc                                                                                                                                 <  4:309547/148‑1 (MQ=255)                    gcCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGc                                                                                                                             <  4:176896/148‑1 (MQ=255)                         tatTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGa                                                                                                                                                  >  2:222936/1‑122 (MQ=255)                         tatTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGa                                                                                                                                                  <  1:222936/122‑1 (MQ=255)                            tgctgcACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGcc                                                                                                                      >  8:385093/1‑147 (MQ=255)                                 cACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGAt                                                                                                                <  4:226274/148‑1 (MQ=255)                                     cTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAACGCTCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTGTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCAcc                                                                                                            >  3:250105/1‑148 (MQ=255)                                     cTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCAcc                                                                                                            >  4:424513/1‑148 (MQ=255)                                        aGCGCCATTACAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCa                                                                                                         <  6:336238/148‑1 (MQ=255)                                          cgcCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGGCATGGCAACGAGAGTCTCGCGTGGGAAGGCGTGCATATCCATCACGCCACCGCCCAGGATCACCGCATc                                                                                                       <  6:411268/148‑1 (MQ=255)                                               ttAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCa                                                                                                              >  5:331252/1‑136 (MQ=255)                                               ttAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCa                                                                                                              <  6:331252/136‑1 (MQ=255)                                                tAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGCCCTGATGCGGCAGTGGACGGCGCAGGAACATTTGGGTCATGGCAACGAGAGTCTCGAGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATc                                                                                                 >  8:94236/1‑148 (MQ=255)                                                tAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGTGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATc                                                                                                 >  5:118844/1‑148 (MQ=255)                                                tAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATc                                                                                                 >  6:140078/1‑148 (MQ=255)                                                   tgTCAGAAGATGAGGCGGCAATAATGCGCACGAACTGATGCGGCAGTTGACGGCGCAGGTACTTTTGGTTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGa                                                                                               <  4:250105/146‑1 (MQ=255)                                                          agaTGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGAtt                                                                                        >  1:15925/1‑147 (MQ=255)                                                          agaTGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTa                                                                                       <  4:333417/148‑1 (MQ=255)                                                           gaTGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTaa                                                                                      <  5:140078/148‑1 (MQ=255)                                                             tGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGAATCACCGCATCggg                                                                                                    >  7:393334/1‑132 (MQ=255)                                                             tGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGGCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGAATCACCGCATCggg                                                                                                    <  8:393334/132‑1 (MQ=255)                                                                 gcggcAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCtgg                                                                                <  2:262768/148‑1 (MQ=255)                                                                   ggcAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCCCGCCACCGCCCAGGATCACCGCATCGGGATCGCACAGATTAATGCtggtg                                                                              >  3:61366/1‑148 (MQ=255)                                                                      aaTAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCTGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCa                                                                           <  5:423101/148‑1 (MQ=255)                                                                      aaTAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCa                                                                           <  3:8856/148‑1 (MQ=255)                                                                      aaTAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCa                                                                           <  5:318642/148‑1 (MQ=255)                                                                       aTAAAGCGCACGACCTGATGCGGCAGTGGATGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCaa                                                                          >  5:83440/1‑148 (MQ=255)                                                                         aaaGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGAt                                                                                                  <  4:140312/122‑1 (MQ=255)                                                                         aaaGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGAt                                                                                                  >  3:140312/1‑122 (MQ=255)                                                                          aaGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATgg                                                                       >  2:288135/1‑148 (MQ=255)                                                                                     cTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCcgcg                                                            >  2:4375/1‑148 (MQ=255)                                                                                         tGCGGCAGTTGACGGCGCATGTACTTTAGGGTAATGGCAACGAGAGTCTCGCGTGTGAAGGCGGGCATATCCATCACTCCACCGCCCATGATCACCGCCTCTGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGtttt                                                        <  8:408126/148‑1 (MQ=255)                                                                                          gCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGAt                                                                                                                >  4:92465/1‑91 (MQ=255)                                                                                          gCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGAt                                                                                                                <  3:92465/91‑1 (MQ=255)                                                                                           cGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCa                                                      <  1:288135/148‑1 (MQ=255)                                                                                               aGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGAAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAg                                                    >  1:336626/1‑146 (MQ=255)                                                                                               aGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGAAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAg                                                    >  1:342965/1‑146 (MQ=255)                                                                                                  ggACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGAc                                               >  7:229689/1‑148 (MQ=255)                                                                                                   gACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGAct                                              <  2:15925/148‑1 (MQ=255)                                                                                                       gcgcAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTg                                           >  4:35222/1‑147 (MQ=255)                                                                                                       gcgcAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTg                                           >  5:332705/1‑147 (MQ=255)                                                                                                         gcAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGAc                                        >  1:139627/1‑148 (MQ=255)                                                                                                          cAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGAc                                        >  3:171319/1‑147 (MQ=255)                                                                                                            ggTACTTTTGTGTCATGGCAACGAGAGTCTCGCGTGGGACGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCATACTCTGGACGaa                                     <  1:302957/148‑1 (MQ=255)                                                                                                            ggTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGaa                                     <  8:229689/148‑1 (MQ=255)                                                                                                            ggTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGaa                                     <  6:118844/148‑1 (MQ=255)                                                                                                             gTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGTCATATCCATCACGCCACCGCCCAGGATCACCGCATCTTGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTTGACGaa                                     <  7:94236/147‑1 (MQ=255)                                                                                                             gTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGa                                      >  5:242117/1‑146 (MQ=255)                                                                                                              tACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAg                                   >  1:467316/1‑148 (MQ=255)                                                                                                              tACTTTTGGGTCATGGCAACGAGAGTCTCGAGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCAGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAg                                   >  6:98250/1‑148 (MQ=255)                                                                                                               aCTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGaaa                                    <  4:354509/146‑1 (MQ=255)                                                                                                               aCTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGaaa                                    >  3:354509/1‑146 (MQ=255)                                                                                                                  tttGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGccc                                                                    >  3:227517/1‑111 (MQ=255)                                                                                                                  tttGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGccc                                                                    <  4:227517/111‑1 (MQ=255)                                                                                                                      gggcATGGCAACGAGAGTCTCGCGTGGAAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGtt                            <  2:342965/144‑1 (MQ=255)                                                                                                                      ggTCATGGCAACGAGAGTCTCGCGTGGAAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGtt                            <  2:336626/147‑1 (MQ=255)                                                                                                                            ggCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCaa                                                     <  1:74416/116‑1 (MQ=255)                                                                                                                            ggCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCaa                                                     >  2:74416/1‑116 (MQ=255)                                                                                                                              cAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATg                   >  2:249517/1‑148 (MQ=255)                                                                                                                                aCGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGa                 >  5:21969/1‑148 (MQ=255)                                                                                                                                aCGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGa                 >  3:419699/1‑148 (MQ=255)                                                                                                                                  gagaGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGaa                                     >  2:171226/1‑126 (MQ=255)                                                                                                                                  gagaGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGaa                                     <  1:171226/126‑1 (MQ=255)                                                                                                                                          cgcgTGGGAAGGCGGGCATATCCATCACGCCAACGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATc       >  7:350505/1‑148 (MQ=255)                                                                                                                                           gcgTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATcg      >  4:298179/1‑148 (MQ=255)                                                                                                                                           gcgTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATc       <  1:4375/147‑1 (MQ=255)                                                                                                                                               gggACGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCAGCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCCGACTCTGGACGAAAGGGGCGGTTTCCGCATGGACGAAAAGATCGCGCa  <  4:94355/148‑1 (MQ=255)                                                                                                                                               gggAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCgcgc   >  8:378873/1‑147 (MQ=255)                                                                                                                                                 |                                                                                                                                                  TGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCA  >  NC_000913/4306894‑4307183

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