Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,182,712	G→A	71.9%	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	71.9%	138.5 / 60.7	92	A55T (GCC→ACC)	zupT	zinc transporter
Reads supporting (aligned to +/- strand):  ref base G (12/14);  new base A (34/32);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCA  >  NC_000913/3182623‑3182807                                                                                          |                                                                                                ggCAAACACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagc                                                                                       <  2:2196122/100‑1 (MQ=255) ggCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagc                                                                                       >  1:22903/1‑100 (MQ=255) ggCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGAGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagc                                                                                       >  1:700911/1‑100 (MQ=255)   cAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcag                                                                                     >  1:2322019/1‑100 (MQ=255)   cAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCCCTagcag                                                                                     >  2:14752/1‑100 (MQ=255)    aaaaaCCCTCGAACCGCTTAGTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCCAACTa                                                                                         >  2:966288/1‑95 (MQ=255)    aaaaaCCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagcagc                                                                                    <  1:1481014/100‑1 (MQ=255)    aaaaaCCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGAGAATCATGTTACACAACTCATTAATGGAAATGCTTCCTGCCGCACTagcagc                                                                                    >  1:786873/1‑100 (MQ=255)     aaaaCCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCt                                                                                   <  1:1710964/100‑1 (MQ=255)     aaaaCCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCt                                                                                   <  2:1310908/100‑1 (MQ=255)     aaaaCCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCt                                                                                   >  2:1432559/1‑100 (MQ=255)      aaaCCCTCTAACCGCTTACTGGCGTTTCCGCTAGGTTTTGCGGCGGTGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTg                                                                                  <  2:251397/100‑1 (MQ=255)        aCCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGaa                                                                                >  1:280059/1‑100 (MQ=255)         cccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCTGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAg                                                                               <  1:777899/100‑1 (MQ=255)          ccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAgg                                                                              <  1:1888346/100‑1 (MQ=255)            tCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGaa                                                                            <  2:884902/100‑1 (MQ=255)             cGAACCGCTTACTGGCGTTTTCGCCAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAAt                                                                           >  2:1330561/1‑100 (MQ=255)              gAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATg                                                                          >  1:805913/1‑100 (MQ=255)                 ccGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCg                                                                       >  2:1813141/1‑100 (MQ=255)                 ccGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCg                                                                       >  2:891353/1‑100 (MQ=255)                 ccGCTTACTGGCGTTTTCGCTAGGTTTTGCGGAGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAg                                                                               <  2:1901805/92‑1 (MQ=255)                  cGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGc                                                                      <  2:2469683/100‑1 (MQ=255)                   gCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGcc                                                                     >  2:1935377/1‑100 (MQ=255)                   gCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGcc                                                                     <  1:687080/100‑1 (MQ=255)                    cTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGcc                                                                     <  2:1159065/99‑1 (MQ=255)                      tACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCtgt                                                                  <  1:1118853/100‑1 (MQ=255)                      tACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCtgt                                                                  >  2:2195792/1‑100 (MQ=255)                        cTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCTCACTAGCAGCTGAAGGAATGTCGACTgtgt                                                                >  1:601152/1‑100 (MQ=255)                        cTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTgtgt                                                                >  2:894813/1‑100 (MQ=255)                          ggCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTg                                                              >  1:2185188/1‑100 (MQ=255)                          ggCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTg                                                              >  2:57385/1‑100 (MQ=255)                          ggCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTg                                                              >  1:1654830/1‑100 (MQ=255)                           gCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTg                                                              <  2:853209/99‑1 (MQ=255)                            cGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTggg                                                            >  1:625222/1‑100 (MQ=255)                            cGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTggg                                                            <  1:2227862/100‑1 (MQ=255)                             gTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGt                                                           <  1:1978076/100‑1 (MQ=255)                             gTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGt                                                           <  2:890191/100‑1 (MQ=255)                             gTTTTCGCCAGGTTTTGCGGCGGGGATCATGTTGCTCATCTTATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGt                                                           <  1:1330574/100‑1 (MQ=255)                              ttttCGCTAGGGTTTGGGGCGGGGATCATGTTGCTCGTCACTTTAATGTAAATGCTTCCTACCGCGCTAGCAGCTGAAGGAATGTCGCCTGTGTTggg                                                            <  2:2459473/98‑1 (MQ=255)                               tttCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTa                                                         >  2:2374710/1‑100 (MQ=255)                                ttCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTAt                                                        >  2:759002/1‑100 (MQ=255)                                ttCGCTAGGTTTTGCAGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTgg                                                             <  2:2394040/95‑1 (MQ=255)                                 tCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATAATTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATg                                                       >  1:718523/1‑100 (MQ=255)                                  cGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATgg                                                      <  2:2321990/100‑1 (MQ=255)                                  cGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATgg                                                      >  2:323716/1‑100 (MQ=255)                                   gCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGa                                                     >  2:2244237/1‑100 (MQ=255)                                          tttGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTa                                              <  1:2441673/100‑1 (MQ=255)                                          tttGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGATGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTa                                              <  2:1488285/100‑1 (MQ=255)                                          tttGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTa                                              <  2:2514255/100‑1 (MQ=255)                                              cggcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCtt                                          <  1:891362/100‑1 (MQ=255)                                              cggcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCtt                                          >  2:88171/1‑100 (MQ=255)                                               ggcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTc                                         >  1:1438271/1‑100 (MQ=255)                                                gcggGGATCATGTTGCTCATCTCATTAATGTAAATGCTTCCTGCCGCACTAGTAGCTGAAGGAATGTCGCCTGTGTTTGGTTATGGAATGTTTATCtt                                          <  2:2343271/98‑1 (MQ=255)                                                gcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCg                                        >  2:2390437/1‑100 (MQ=255)                                                gcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCg                                        <  2:1255298/100‑1 (MQ=255)                                                 cggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCgg                                       <  1:1597432/100‑1 (MQ=255)                                                  ggggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGTCTGTGTTGGGTTATGGAATGTTTATCTTCgg                                       <  2:1218675/99‑1 (MQ=255)                                                  ggggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGt                                      >  2:2018914/1‑100 (MQ=255)                                                   gggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCTCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTc                                     <  2:2107818/100‑1 (MQ=255)                                                   gggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTc                                     <  1:323718/100‑1 (MQ=255)                                                      aTCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGc                                  >  2:15141/1‑100 (MQ=255)                                                       tCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTCTATCTTCGGTCTGCt                                 >  1:2618754/1‑100 (MQ=255)                                                          tGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTgg                              <  2:575426/100‑1 (MQ=255)                                                             tGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCt                            >  2:2570287/1‑99 (MQ=255)                                                                 catctcatTCATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGGCTGTGCTGGGTTATGGCATGTTTATTTTCGGTCTGCTTGCCTAtttt                       >  2:299642/1‑100 (MQ=255)                                                                  atctcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTg                      <  1:898650/100‑1 (MQ=255)                                                                   tctcatTAATGGAAATGCTTCCTGCCGCACTAGCAGTTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTgg                     <  1:265849/100‑1 (MQ=255)                                                                   tctcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTgg                     <  2:470491/100‑1 (MQ=255)                                                                     tcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGcc                   <  2:884816/100‑1 (MQ=255)                                                                     tcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGcc                   <  1:1663117/100‑1 (MQ=255)                                                                      catTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCt                  >  1:1437883/1‑100 (MQ=255)                                                                      catTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCt                  >  1:1226612/1‑100 (MQ=255)                                                                      catTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTGATCTTCGGTCTGCTTGGCTATTTTGGCCt                  >  2:2250621/1‑100 (MQ=255)                                                                       atTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGTAATGTCGCCTGTGTTGGGTGATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTg                 <  2:1442304/100‑1 (MQ=255)                                                                        ttAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTgg                <  1:14752/100‑1 (MQ=255)                                                                         tAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGa               >  2:940104/1‑100 (MQ=255)                                                                          aaTGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGAc              >  1:1781909/1‑100 (MQ=255)                                                                           aTGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGAcc             >  2:1913405/1‑100 (MQ=255)                                                                           aTGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGAcc             >  1:2582652/1‑100 (MQ=255)                                                                           aTGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGAcc             <  1:1720471/100‑1 (MQ=255)                                                                             ngAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGc           >  1:143351/2‑100 (MQ=255)                                                                             ggAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGc           <  2:124679/100‑1 (MQ=255)                                                                             ggAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGc           >  1:2135025/1‑100 (MQ=255)                                                                             ggAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGc           >  1:2095048/1‑100 (MQ=255)                                                                              gAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCa          <  1:124850/100‑1 (MQ=255)                                                                              gAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCTTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCa          <  1:57385/100‑1 (MQ=255)                                                                                aaTGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATg        >  2:1115793/1‑100 (MQ=255)                                                                                   gCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTg     <  2:2050486/100‑1 (MQ=255)                                                                                   gCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTg     <  2:757167/100‑1 (MQ=255)                                                                                   gCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTg     <  1:1432573/100‑1 (MQ=255)                                                                                     ttCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGcc   >  1:985400/1‑100 (MQ=255)                                                                                      tCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGCACCGCATGTTGCca  >  2:1802833/1‑100 (MQ=255)                                                                                          |                                                                                                GGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCA  >  NC_000913/3182623‑3182807

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