Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	410,479	A→G	100%	S112S (TCA→TCG)	mak →	manno(fructo)kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	410,479	0	A	G	98.9%	319.3 / ‑4.6	88	S112S (TCA→TCG)	mak	manno(fructo)kinase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (45/42);  total (45/43)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

CAACCTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCG  >  NC_000913/410382‑410575                                                                                                  |                                                                                                 cAACCTGGCTCAACGGTCAGCCATTCGATCAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGGTTCGGa                                                                                                >  2:954233/1‑100 (MQ=255) cAACCTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGa                                                                                                <  2:1033443/100‑1 (MQ=255)   aCCTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAg                                                                                              >  1:196021/1‑100 (MQ=255)   aCCTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAg                                                                                              >  1:2067481/1‑100 (MQ=255)    ccTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGc                                                                                             <  2:1212246/100‑1 (MQ=255)     cTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCa                                                                                            <  1:147793/100‑1 (MQ=255)        gCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTa                                                                                         <  1:1932251/100‑1 (MQ=255)        gCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTa                                                                                         <  1:1326585/100‑1 (MQ=255)        gCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTa                                                                                         <  2:1403302/100‑1 (MQ=255)         cTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAg                                                                                        >  1:1048217/1‑100 (MQ=255)         cTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAg                                                                                        >  1:1053760/1‑100 (MQ=255)            aaCGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATg                                                                                     <  2:754737/100‑1 (MQ=255)              cGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGgc                                                                                   >  2:2331959/1‑100 (MQ=255)                gTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGgcgc                                                                                 >  2:1775268/1‑100 (MQ=255)                  cAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGcgcg                                                                                >  2:1475617/1‑99 (MQ=255)                   aGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGc                                                                              <  1:1096386/100‑1 (MQ=255)                    gCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCa                                                                             >  1:1230299/1‑100 (MQ=255)                     ccATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAg                                                                            >  2:2147154/1‑100 (MQ=255)                     ccATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAg                                                                            >  1:2076055/1‑100 (MQ=255)                      cATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGc                                                                           >  2:183480/1‑100 (MQ=255)                       aTTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCg                                                                          <  2:1435478/100‑1 (MQ=255)                       aTTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCg                                                                          >  1:462813/1‑100 (MQ=255)                       aTTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCg                                                                          >  1:1479514/1‑100 (MQ=255)                           gATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAg                                                                      >  1:1735121/1‑100 (MQ=255)                             tAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAgcg                                                                    <  1:775667/100‑1 (MQ=255)                              aaaGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAgcgc                                                                   <  2:1529237/100‑1 (MQ=255)                              aaaGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAgcgc                                                                   >  2:1062721/1‑100 (MQ=255)                              aaaGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAgcgc                                                                   <  2:1795780/100‑1 (MQ=255)                                aGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAg                                                                 <  2:284377/100‑1 (MQ=255)                                 gACTTAAGCGCGAGGTTGCATCGGGAAGTCCGGCTCGCAAATGACTCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCa                                                                  <  2:1587277/98‑1 (MQ=255)                                  aCTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGAc                                                               <  1:1764104/100‑1 (MQ=255)                                   cTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACg                                                              <  1:809802/100‑1 (MQ=255)                                     tAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGt                                                            >  2:2202037/1‑100 (MQ=255)                                       aGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTAt                                                          >  2:925491/1‑100 (MQ=255)                                        gcgcGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTAtt                                                         <  2:2067465/100‑1 (MQ=255)                                        gcgcGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTAt                                                          <  2:2221764/99‑1 (MQ=255)                                        gcgcGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGAAGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTAtt                                                         >  1:1144484/1‑100 (MQ=255)                                           cgAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGc                                                      >  1:53246/1‑100 (MQ=255)                                           cgAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGc                                                      >  2:1300984/1‑100 (MQ=255)                                            gAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGTCGCTAAGTGACTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGcc                                                     <  2:1651812/100‑1 (MQ=255)                                              ggTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGt                                                   <  1:1869179/100‑1 (MQ=255)                                               gTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTg                                                  >  1:1667487/1‑100 (MQ=255)                                               gTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTg                                                  <  2:1668044/100‑1 (MQ=255)                                                 tGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGAt                                                <  2:1503008/100‑1 (MQ=255)                                                  gCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGAGGGCGCGACAGCGGGAGCGCAGACGGTATTTGCCGTGAtt                                               >  2:1654413/1‑100 (MQ=255)                                                   cAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTa                                              >  2:1769031/1‑100 (MQ=255)                                                    aGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTAt                                             >  2:578985/1‑100 (MQ=255)                                                     gCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATc                                            <  1:825154/100‑1 (MQ=255)                                                     gCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATc                                            >  2:554776/1‑100 (MQ=255)                                                         gAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCa                                        >  2:1002854/1‑100 (MQ=255)                                                         gAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCa                                        >  2:894852/1‑100 (MQ=255)                                                           aGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACg                                      <  2:771181/100‑1 (MQ=255)                                                             tGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACggg                                    <  2:663012/100‑1 (MQ=255)                                                             tGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACggg                                    <  2:835769/100‑1 (MQ=255)                                                              gCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGa                                   >  1:1476693/1‑100 (MQ=255)                                                              gCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGa                                   >  2:1709977/1‑100 (MQ=255)                                                               cGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGAt                                  <  1:1446114/100‑1 (MQ=255)                                                               cGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGAt                                  <  1:244049/100‑1 (MQ=255)                                                                ggCTGGCAAATGACGCTAACTGTCCGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATg                                 >  2:1711760/1‑100 (MQ=255)                                                                    ggCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATgcggc                             >  1:1051689/1‑100 (MQ=255)                                                                    ggCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATgcggc                             >  1:234218/1‑100 (MQ=255)                                                                         aTGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCggg                        <  2:2050540/100‑1 (MQ=255)                                                                         aTGACGCGAACTGTCTGGTGGTTTCGGAAGCAGTAGATGGCGGGGCAGCGGGAGCGCAGACGTTATTTGCCGTGATTATCGGCACGGGATGCGGCGCggn                        <  2:1263254/100‑2 (MQ=255)                                                                           gACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCg                      <  2:1388474/100‑1 (MQ=255)                                                                           gACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCg                      >  2:260522/1‑100 (MQ=255)                                                                            aCGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGt                     <  1:236850/100‑1 (MQ=255)                                                                             cGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTg                    >  1:1180287/1‑100 (MQ=255)                                                                             cGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTg                    <  2:1495551/100‑1 (MQ=255)                                                                              gCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCTGTGATTATCGGCACGGGATGCGGCGCGGGCGTgg                   <  2:118991/100‑1 (MQ=255)                                                                              gCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTgg                   <  1:1300990/100‑1 (MQ=255)                                                                               cTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGctcgtg                  <  1:158711/100‑6 (MQ=255)                                                                                 aaCTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCGGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCa                 >  2:2139165/1‑99 (MQ=255)                                                                                     gTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCaa            <  1:1484000/100‑1 (MQ=255)                                                                                     gTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCaa            >  1:1986895/1‑100 (MQ=255)                                                                                     gTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCaa            >  2:507993/1‑100 (MQ=255)                                                                                      tCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAAt           <  1:470880/100‑1 (MQ=255)                                                                                       cTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATg          <  2:1327818/100‑1 (MQ=255)                                                                                       cTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATg          >  2:1112228/1‑100 (MQ=255)                                                                                       cTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATg          <  2:1818803/100‑1 (MQ=255)                                                                                        tGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATgg         >  2:442974/1‑100 (MQ=255)                                                                                        tGGCGGTGTCGGAAGCAGTAGATGGCCCGGCAGCGGGAGCGCATACGGGATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATg          <  2:1036834/99‑1 (MQ=255)                                                                                         ggCGGTTTCGGAAGCATTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATggg        <  2:1377559/100‑1 (MQ=255)                                                                                         ggCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATggg        >  2:547307/1‑100 (MQ=255)                                                                                         ggCGGTTTCGGAAACAGTAGATGGCGCGGCATCGGGAGCGCAGACGGTACTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATggg        >  1:1581653/1‑100 (MQ=255)                                                                                          gCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATgggc       <  1:529616/100‑1 (MQ=255)                                                                                             gTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATgggcggg    >  2:2228486/1‑100 (MQ=255)                                                                                              tttCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATgggcgggc   >  1:1247875/1‑100 (MQ=255)                                                                                               ttCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGggcgggcg  >  1:225846/1‑100 (MQ=255)                                                                                                  |                                                                                                 CAACCTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCG  >  NC_000913/410382‑410575

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