Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	263,172	A→C	100%	pseudogene (41/63 nt)	ykfN ←	pseudogene, YdiA family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	263,172	0	A	C	96.0%	309.3 / 0.1	92	pseudogene (41/63 nt)	ykfN	pseudogene, YdiA family
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base C (48/40);  total (49/43)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTTCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTAC  >  NC_000913/263075‑263267                                                                                                  |                                                                                                ggCTTCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGc                                                                                               <  2:1338079/100‑1 (MQ=255)  gCTTCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCt                                                                                              <  2:1263414/100‑1 (MQ=255)     tCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCAt                                                                                           >  1:537477/1‑100 (MQ=255)     tCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCAt                                                                                           >  1:898552/1‑100 (MQ=255)      cGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCAtt                                                                                          >  2:1554832/1‑100 (MQ=255)      cGCTCACTATACAACCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCAtt                                                                                          >  2:1699250/1‑100 (MQ=255)       gCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTa                                                                                         >  1:2001253/1‑100 (MQ=255)       gCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTa                                                                                         <  1:1871628/100‑1 (MQ=255)       gCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTa                                                                                         >  2:2053054/1‑100 (MQ=255)        cTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTAt                                                                                        >  2:2270105/1‑100 (MQ=255)            cTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGTTCGAACAAAATCGCTCCTTACGAt                                                                                     >  1:2592966/1‑99 (MQ=255)            cTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATc                                                                                    <  1:779525/100‑1 (MQ=255)             tataCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCg                                                                                   <  1:80331/100‑1 (MQ=255)             tataCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCg                                                                                   <  1:1957582/100‑1 (MQ=255)             tataCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGGGGGGGGTGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCg                                                                                   <  1:773464/100‑1 (MQ=255)              ataCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCgg                                                                                  >  2:381612/1‑100 (MQ=255)               taCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGt                                                                                 >  1:1227367/1‑100 (MQ=255)                aCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGtt                                                                                <  1:1989853/100‑1 (MQ=255)                aCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGtt                                                                                >  1:215333/1‑100 (MQ=255)                    ccTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATa                                                                            <  1:2694774/100‑1 (MQ=255)                    ccTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATa                                                                            >  1:117043/1‑100 (MQ=255)                     cTTTGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTa                                                                               <  2:1058504/96‑1 (MQ=255)                      ttGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACt                                                                          <  1:1470178/100‑1 (MQ=255)                       tGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACtt                                                                         >  2:2340396/1‑100 (MQ=255)                         gCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCa                                                                       <  1:99239/100‑1 (MQ=255)                          cTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCaa                                                                      >  1:1490083/1‑100 (MQ=255)                          cTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCaa                                                                      <  2:461494/100‑1 (MQ=255)                           tgcgagTTTAGTTGTTGGCCTCTCCTTATTTTAATGTGTTGGCAATCTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGATTAATACATCatc                                                                     <  2:1697597/95‑3 (MQ=255)                           tGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAAc                                                                     >  2:1426215/1‑100 (MQ=255)                            gCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAAAAAATTCGCTCATTATGATCGGTTAATACTTCaa                                                                      >  2:1130944/1‑98 (MQ=255)                               ggTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCt                                                                 >  1:141080/1‑100 (MQ=255)                               ggTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTAGCTCATTATGATCGGTTAATACTTCAACTTCt                                                                 <  1:2606187/100‑1 (MQ=255)                                gTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTg                                                                >  1:2597085/1‑100 (MQ=255)                                  ttAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGt                                                              >  2:449393/1‑100 (MQ=255)                                  ttAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGt                                                              <  1:2516309/100‑1 (MQ=255)                                  ttAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGt                                                              >  2:2105167/1‑100 (MQ=255)                                    aGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTg                                                            <  1:2129402/100‑1 (MQ=255)                                     gTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGc                                                           >  2:149728/1‑100 (MQ=255)                                      ttGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCa                                                          <  1:1126408/100‑1 (MQ=255)                                      ttGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCa                                                          <  1:2270538/100‑1 (MQ=255)                                       tGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCAt                                                         >  2:687925/1‑100 (MQ=255)                                       tGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCAt                                                         <  1:633630/100‑1 (MQ=255)                                        gTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATg                                                        <  1:2046222/100‑1 (MQ=255)                                        gTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATg                                                        <  2:2419394/100‑1 (MQ=255)                                           caccacTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAtt                                                     >  1:1510718/1‑100 (MQ=255)                                           caccacTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAtt                                                     >  2:1839965/1‑100 (MQ=255)                                           caccacTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAtt                                                     >  1:215421/1‑100 (MQ=255)                                           caccacTCCTAAATTTAATGT‑‑‑‑‑CAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgtttg                                                <  1:360809/100‑1 (MQ=255)                                            accacTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttg                                                    <  2:2724448/100‑1 (MQ=255)                                            accacTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttg                                                    >  2:1847595/1‑100 (MQ=255)                                               acTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgttt                                                 <  1:2154308/100‑1 (MQ=255)                                               acTCCTAAATTTAATGTGTTGGCAATGTGCTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgttt                                                 <  2:1649916/100‑1 (MQ=255)                                                cTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgtttg                                                >  2:1176420/1‑100 (MQ=255)                                                cTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgtt                                                  <  2:1329743/98‑1 (MQ=255)                                                cTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCGCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgtt                                                  >  1:543865/1‑98 (MQ=255)                                                  ccTAGATTGAATGTGTTGGCAATGTGTTCTATAAACGTCGAACAAATTCGCTCATGATGTGCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTc                                              <  2:2467972/100‑1 (MQ=255)                                                   cTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTcc                                             <  1:1478239/100‑1 (MQ=255)                                                    tAAATGTAATTTGTTGGCAATTTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTc                                              <  2:2604335/98‑1 (MQ=255)                                                      aaTTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTa                                          >  1:2257576/1‑100 (MQ=255)                                                       aTTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTaa                                         >  1:2475564/1‑100 (MQ=255)                                                          tAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTaaaaa                                      >  1:1836614/1‑100 (MQ=255)                                                          tAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTaaaaa                                      <  2:548517/100‑1 (MQ=255)                                                           aaTGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAg                                     <  1:627680/100‑1 (MQ=255)                                                              ttcttcgCAATGTGTTCAATAAAGCTCGAACAAATTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATa                                  <  2:1212461/94‑1 (MQ=255)                                                               tgtTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATaa                                 >  2:2173932/1‑100 (MQ=255)                                                                   ggCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAAcgcg                             <  1:1946981/100‑1 (MQ=255)                                                                     cAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGcc                           <  1:593381/100‑1 (MQ=255)                                                                      aaTGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCt                          <  2:483729/100‑1 (MQ=255)                                                                      aaTGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGATAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACACGCCt                          >  2:504899/1‑100 (MQ=255)                                                                       aTGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTg                         >  1:1048271/1‑100 (MQ=255)                                                                       aTGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTAGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTg                         <  1:1935748/100‑1 (MQ=255)                                                                        tgtgTTCAATAAAGCTCGAACAAATTCGCTCATAATGATCGGTTAATACTTCAACTTCTGGTTGCATCATTGTTTGTCCGTAAAAAGATAACGCGCCTGc                        >  2:1815978/1‑100 (MQ=255)                                                                          tgtTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCg                      >  1:1820531/1‑100 (MQ=255)                                                                          tgtTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCg                      >  1:1613978/1‑100 (MQ=255)                                                                           gtTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCgg                     <  1:942953/100‑1 (MQ=255)                                                                             tCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgt                   >  1:1220159/1‑100 (MQ=255)                                                                              cAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgta                  >  1:584948/1‑100 (MQ=255)                                                                                 tAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgtagta               >  1:2318213/1‑100 (MQ=255)                                                                                  aaaGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGtagtag              <  2:1096192/100‑1 (MQ=255)                                                                                   aaGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGc             >  1:386714/1‑100 (MQ=255)                                                                                     gCTCGAACAAATTCGCTCATTATGTTCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAg           >  2:599217/1‑100 (MQ=255)                                                                                      cTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAgg          >  1:1875504/1‑100 (MQ=255)                                                                                      cTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAgg          <  2:148881/100‑1 (MQ=255)                                                                                      cTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAgg          >  1:2336965/1‑100 (MQ=255)                                                                                      cTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAg           >  2:2184226/1‑99 (MQ=255)                                                                                          aaCAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCTGGTAGTAGCAGGCGCa      >  2:1610557/1‑100 (MQ=255)                                                                                           aCAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCAt     <  2:100864/100‑1 (MQ=255)                                                                                           aCAAATTCGCTCATTATGATCGGTTAATACTTAAACTTCTGGTTGCATGATTGTTTGTACGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCAt     >  2:67137/1‑100 (MQ=255)                                                                                            cAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCAtt    >  2:1013570/1‑100 (MQ=255)                                                                                              aaTTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTAc  >  2:205811/1‑100 (MQ=255)                                                                                              aaTTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTAc  <  1:280641/100‑1 (MQ=255)                                                                                              aaTTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTAc  <  2:386714/100‑1 (MQ=255)                                                                                                  |                                                                                                GGCTTCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTAC  >  NC_000913/263075‑263267

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