Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,182,712	G→A	A55T (GCC→ACC)	zupT →	zinc transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,182,712	0	G	A	96.8%	293.4 / ‑0.1	94	A55T (GCC→ACC)	zupT	zinc transporter
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (42/49);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

GTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCAC  >  NC_000913/3182616‑3182808                                                                                                 |                                                                                                 gTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCg                                                                                               <  1:1760006/100‑1 (MQ=255)  ttCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGc                                                                                              >  1:240653/1‑100 (MQ=255)   tctcGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCa                                                                                             >  1:1514911/1‑100 (MQ=255)   tctcGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCa                                                                                             <  2:1609974/100‑1 (MQ=255)    ctcGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCAc                                                                                            <  1:660745/100‑1 (MQ=255)    ctcGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCAc                                                                                            <  2:2442594/100‑1 (MQ=255)    ctcGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCAc                                                                                            >  1:2262385/1‑100 (MQ=255)        ggCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagc                                                                                        >  1:1602320/1‑100 (MQ=255)        ggCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagc                                                                                        >  1:891757/1‑100 (MQ=255)          cAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcag                                                                                      >  1:2229197/1‑100 (MQ=255)          cAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagcag                                                                                      >  1:2437325/1‑100 (MQ=255)           aaaaaCCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagcagc                                                                                     <  1:375794/100‑1 (MQ=255)            aaaaCCCTCGAACCGCTTACTGGCGTTTTCGTTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCt                                                                                    <  1:1318983/100‑1 (MQ=255)               aCCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGa                                                                                  >  2:227709/1‑99 (MQ=255)                cccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAg                                                                                <  1:794906/100‑1 (MQ=255)                cccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAg                                                                                >  1:92565/1‑100 (MQ=255)                cccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAg                                                                                <  1:1078478/100‑1 (MQ=255)                 ccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAgg                                                                               <  2:2302095/100‑1 (MQ=255)                 ccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAgg                                                                               <  1:1808021/100‑1 (MQ=255)                 ccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAgg                                                                               <  2:545394/100‑1 (MQ=255)                 ccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAgg                                                                               <  2:1514898/100‑1 (MQ=255)                 ccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAg                                                                                >  2:372217/1‑99 (MQ=255)                 ccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCAGCTCATTAATGGAAATGCTTGCTACCGCACTAGCAGCTGAAg                                                                                <  2:483001/99‑1 (MQ=255)                  cTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGa                                                                              <  1:668944/100‑1 (MQ=255)                   tCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGaa                                                                             <  2:1976048/100‑1 (MQ=255)                     gAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATg                                                                           <  2:2337255/100‑1 (MQ=255)                     gAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATg                                                                           <  1:1616548/100‑1 (MQ=255)                       aCCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTc                                                                         >  2:2152501/1‑100 (MQ=255)                        ccGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCAAGTTGCTCATCTCATTAATGCAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCg                                                                        >  2:1660673/1‑100 (MQ=255)                          gCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGACATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGcc                                                                      <  2:2418707/100‑1 (MQ=255)                          gCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGcc                                                                      <  1:2353433/100‑1 (MQ=255)                          gCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGcc                                                                      <  2:1492119/100‑1 (MQ=255)                          gCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGcc                                                                      >  1:1845763/1‑100 (MQ=255)                          gCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGcc                                                                      >  2:1339801/1‑100 (MQ=255)                               cTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTgtgt                                                                 >  1:2137369/1‑100 (MQ=255)                               cTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTgtgt                                                                 >  1:1255026/1‑100 (MQ=255)                                  gCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTgg                                                              <  1:1129888/100‑1 (MQ=255)                                  gCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTgg                                                              <  1:1153303/100‑1 (MQ=255)                                   cGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTggg                                                             <  2:2229177/100‑1 (MQ=255)                                   cGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTggg                                                             <  1:2730420/100‑1 (MQ=255)                                    gTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTTATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTggg                                                             <  2:2089733/99‑1 (MQ=255)                                    gTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGt                                                            <  2:69138/100‑1 (MQ=255)                                         cGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATgg                                                       >  2:2523995/1‑100 (MQ=255)                                         cGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATgg                                                       >  2:2573999/1‑100 (MQ=255)                                         cGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATgg                                                       <  2:2649326/100‑1 (MQ=255)                                         cGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATgg                                                       <  2:120949/100‑1 (MQ=255)                                          gCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGa                                                      >  2:2328461/1‑100 (MQ=255)                                          gCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTAggatcggaag                                                      >  1:2283193/1‑92 (MQ=255)                                           cTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGaa                                                     <  1:1339814/100‑1 (MQ=255)                                           cTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCCGTGCTGGGTTATgg                                                       >  2:1718042/1‑98 (MQ=255)                                            tAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAAt                                                    <  2:2437303/100‑1 (MQ=255)                                             aGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAAGGTCGCCTGTGTTGGGTTATGGAATg                                                   >  2:121466/1‑100 (MQ=255)                                              ggTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGt                                                  <  2:2763655/100‑1 (MQ=255)                                              ggTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGt                                                  <  1:1391246/100‑1 (MQ=255)                                               gTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGtt                                                 >  1:907282/1‑100 (MQ=255)                                                ttttGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGttt                                                <  1:865168/100‑1 (MQ=255)                                                 tttGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTa                                               <  2:669693/100‑1 (MQ=255)                                                  ttGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTa                                               <  2:700487/99‑1 (MQ=255)                                                   tGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATc                                             >  2:2069687/1‑100 (MQ=255)                                                     cggcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCtt                                           <  1:2474096/100‑1 (MQ=255)                                                       gcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCg                                         >  2:1354465/1‑100 (MQ=255)                                                       gcggGGATCATGTTGCTCAACTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCg                                         >  2:1456072/1‑100 (MQ=255)                                                          gggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTc                                      <  1:1768741/100‑1 (MQ=255)                                                           ggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCt                                     <  1:202787/100‑1 (MQ=255)                                                           ggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCt                                     >  1:355620/1‑100 (MQ=255)                                                           ggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCt                                     >  1:780639/1‑100 (MQ=255)                                                                aTGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTg                                <  1:1253987/100‑1 (MQ=255)                                                                aTGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTg                                >  1:1830884/1‑100 (MQ=255)                                                                aTGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTg                                <  1:2152520/100‑1 (MQ=255)                                                                 tGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTgg                               <  1:2528218/100‑1 (MQ=255)                                                                  tttGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGc                              <  1:2977831/99‑1 (MQ=255)                                                                  gTTGCTCATCTCATTAATGGAAATGCATCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGc                              <  2:252035/100‑1 (MQ=255)                                                                    tGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTa                            <  1:2547798/100‑1 (MQ=255)                                                                     gCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTAt                           >  2:2104677/1‑100 (MQ=255)                                                                      ctcatctcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTAtt                          >  1:2998147/1‑100 (MQ=255)                                                                       tcatctcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTAttt                         >  1:950066/1‑100 (MQ=255)                                                                       tcatctcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTAttt                         >  1:2971308/1‑100 (MQ=255)                                                                        catctcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTAtttt                        >  2:1356428/1‑100 (MQ=255)                                                                        catctcatTAATGAAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGCGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTAtttt                        <  2:2533274/100‑1 (MQ=255)                                                                            tccgatctGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGcc                    <  2:608095/93‑1 (MQ=255)                                                                            tcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGcc                    >  1:2482341/1‑100 (MQ=255)                                                                             catTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCt                   >  2:1442284/1‑100 (MQ=255)                                                                             catTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCt                   >  1:694263/1‑100 (MQ=255)                                                                              atTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGAAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTg                  >  2:1865526/1‑100 (MQ=255)                                                                               ttAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTgg                 <  2:2239007/100‑1 (MQ=255)                                                                                tAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGa                >  1:980877/1‑100 (MQ=255)                                                                                    ggAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGccagatcgga            >  1:608102/1‑92 (MQ=255)                                                                                     gcaaTGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCa           <  2:1845747/98‑1 (MQ=255)                                                                                     gAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCa           <  1:754422/100‑1 (MQ=255)                                                                                         tGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGtt       <  2:1215922/100‑1 (MQ=255)                                                                                          gCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTg      <  1:2069706/100‑1 (MQ=255)                                                                                             tCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCca   >  2:2070373/1‑100 (MQ=255)                                                                                              ccTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGtggcaac  >  2:2911382/1‑94 (MQ=255)                                                                                              ccTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCcac  >  2:1392204/1‑100 (MQ=255)                                                                                                 |                                                                                                 GTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCAC  >  NC_000913/3182616‑3182808

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