Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,017,296	A→C	100%	I335S (ATC→AGC)	ycbZ ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,017,296	0	A	C	97.8%	315.5 / ‑1.7	91	I335S (ATC→AGC)	ycbZ	putative peptidase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (36/53);  total (38/53)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

AAAGCGCGTGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGATCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGCCATCAC  >  NC_000913/1017200‑1017392                                                                                                 |                                                                                                 aaaGCGCGTGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGctct                                                                                               <  1:1089569/100‑1 (MQ=255)     cgcgTGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATc                                                                                           >  2:2025513/1‑100 (MQ=255)       cttGGATGACGCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTGCATTCAGGATTTGCTCCTGAAGGCTCTCATCCt                                                                                         <  2:1684076/98‑1 (MQ=255)       cgTGGATGACCCGGAAATTCAATGACCGAAAGGTCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGCTTTGCTCCTGAAGGCTCTCATCCt                                                                                         >  1:397134/1‑100 (MQ=255)       cgTGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCt                                                                                         >  1:1229580/1‑100 (MQ=255)        gTGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTg                                                                                        <  1:2026165/100‑1 (MQ=255)         tGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCtag                                                                                       >  1:1684088/1‑98 (MQ=255)         tGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGc                                                                                       >  1:1275754/1‑100 (MQ=255)          ggATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCa                                                                                      <  1:566289/100‑1 (MQ=255)           gATGACCCGGAAATTGAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCAt                                                                                     <  1:2519687/100‑1 (MQ=255)            aTGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATa                                                                                    <  1:2541417/100‑1 (MQ=255)             tGACCCGGAAAATCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATa                                                                                    <  1:2300996/99‑1 (MQ=255)              gACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACg                                                                                  <  1:1251148/100‑1 (MQ=255)                cccGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATAAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGtt                                                                                <  2:1364679/100‑1 (MQ=255)                  cGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCa                                                                              >  2:2207502/1‑100 (MQ=255)                  cGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCa                                                                              <  2:315275/100‑1 (MQ=255)                   ggAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAg                                                                             <  2:2356580/100‑1 (MQ=255)                   ggAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAg                                                                             >  2:926750/1‑100 (MQ=255)                    gAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGc                                                                            <  1:635955/100‑1 (MQ=255)                    gAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGc                                                                            <  1:232348/100‑1 (MQ=255)                    gAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGc                                                                            <  1:1276974/100‑1 (MQ=255)                       aTTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAg                                                                         <  2:1148367/100‑1 (MQ=255)                       aTTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAg                                                                         >  1:735415/1‑100 (MQ=255)                        ttCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAgg                                                                        <  2:1937891/100‑1 (MQ=255)                        ttCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAgg                                                                        <  1:1852864/100‑1 (MQ=255)                        ttCAATGACCGAAAGGGCGTTAATTGGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAgg                                                                        <  2:2279820/100‑1 (MQ=255)                         tCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGa                                                                       <  1:2548833/100‑1 (MQ=255)                           aaTGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGaaa                                                                     <  2:874662/100‑1 (MQ=255)                            aTGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAAc                                                                    >  1:1979776/1‑100 (MQ=255)                             tGACCGGAAGGTCGTTAATTTGCCTGACTCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCAGCCTGCATACGTTCAGCGAGGAAAcc                                                                   <  2:2349798/100‑1 (MQ=255)                             tGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAAcc                                                                   <  2:778889/100‑1 (MQ=255)                             tGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAAcc                                                                   <  2:2430235/100‑1 (MQ=255)                               acggAAAGGGCGTTAATTGGCCGGTTGTGTGCGTATTAGGTTTGAATTAGGATTTGCTCCTGAAGGCTATCATCCTGCATACGTTCAGCGAGTAAGCCtt                                                                 <  2:74700/97‑1 (MQ=255)                                ccGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCtt                                                                 <  1:573402/99‑1 (MQ=255)                                ccGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTc                                                                >  2:1470211/1‑100 (MQ=255)                                 tccgatctcGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTcg                                                               <  2:647158/92‑1 (MQ=255)                                   aaaGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTcgcg                                                             >  2:789801/1‑100 (MQ=255)                                   aaaGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTcgcg                                                             <  1:2207517/100‑1 (MQ=255)                                   aaaGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTcgcg                                                             <  1:2479944/100‑1 (MQ=255)                                     aGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGcc                                                           >  1:1335572/1‑100 (MQ=255)                                        gCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCAt                                                         <  2:121237/99‑1 (MQ=255)                                         cGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTcgagatcgga                                                       >  1:647165/1‑92 (MQ=255)                                         cGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTc                                                       >  1:1183996/1‑100 (MQ=255)                                         cGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTc                                                       >  1:2469507/1‑100 (MQ=255)                                            tAATTTGGCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCAtt                                                        <  2:1176048/96‑1 (MQ=255)                                                ttGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACgctgc                                                <  1:923654/100‑1 (MQ=255)                                                  gCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGctgctg                                              >  2:1438528/1‑100 (MQ=255)                                                  gCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGctgctg                                              >  1:2502344/1‑100 (MQ=255)                                                   cccGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCACCATTCACGCtgctgc                                             <  1:345934/100‑1 (MQ=255)                                                    ccGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTgca                                            >  2:1916343/1‑100 (MQ=255)                                                     cGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCtgctgcgg                                           <  1:1680357/100‑3 (MQ=255)                                                     cGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTgcag                                           <  1:2347980/100‑1 (MQ=255)                                                      gATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGCAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCAGCTgca                                            <  2:2121298/98‑1 (MQ=255)                                                      gATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTgcagc                                          >  2:208828/1‑100 (MQ=255)                                                         atgagagCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTgcagca                                         <  2:302299/92‑1 (MQ=255)                                                          cGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTa                                      <  1:982172/100‑1 (MQ=255)                                                            ttCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAg                                    <  2:1229572/100‑1 (MQ=255)                                                            cgatctttCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAg                                    <  2:2261730/94‑1 (MQ=255)                                                             tCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGt                                   <  1:660691/100‑1 (MQ=255)                                                              cGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGtt                                  >  1:65716/1‑100 (MQ=255)                                                               gCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGttt                                 <  1:649953/100‑1 (MQ=255)                                                                 cTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAg                               >  1:1285163/1‑100 (MQ=255)                                                                  ttCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAgagatcg                              >  1:2261744/1‑94 (MQ=255)                                                                  ttCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAgc                              <  2:1498233/100‑1 (MQ=255)                                                                  ttCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAgc                              >  2:2361033/1‑100 (MQ=255)                                                                  ttCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAgc                              <  1:815163/100‑1 (MQ=255)                                                                    cGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAgctg                            >  2:2255179/1‑100 (MQ=255)                                                                    cGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAgct                             >  2:55506/1‑99 (MQ=255)                                                                    cGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAgccg                            <  1:1734521/100‑3 (MQ=255)                                                                      gTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAgctgct                          <  1:793815/100‑1 (MQ=255)                                                                         tCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGc                       >  2:301364/1‑100 (MQ=255)                                                                           aaTCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGCGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCg                     >  2:756808/1‑100 (MQ=255)                                                                           aaTCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCg                     >  2:619325/1‑100 (MQ=255)                                                                            aTCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCgg                    <  2:407151/100‑1 (MQ=255)                                                                            aTCAGGATTTGCTCCTGAAGGATCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCgg                    >  2:1221572/1‑100 (MQ=255)                                                                            aTCAGGAATTGATCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCGCGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCgg                    >  2:1341603/1‑100 (MQ=255)                                                                             tCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGga                   >  1:2388554/1‑100 (MQ=255)                                                                              cAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGgag                  <  2:1720768/100‑1 (MQ=255)                                                                               aGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGgaga                 <  1:1047933/100‑1 (MQ=255)                                                                               aGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGgag                  >  1:1571755/1‑99 (MQ=255)                                                                                ggATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCCCTAAGTTTAGCCGCTCGcacgagga                >  2:2355059/1‑93 (MQ=255)                                                                                ggATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGaa                <  1:1854946/100‑1 (MQ=255)                                                                                  aTTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAgg              <  2:366283/100‑1 (MQ=255)                                                                                  aTTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAgg              <  2:1477718/100‑1 (MQ=255)                                                                                      gCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGgtggc          <  1:1373078/100‑3 (MQ=255)                                                                                        tCCTGAAGGATCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGc        >  1:1954567/1‑100 (MQ=255)                                                                                         ccTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGcc       >  2:1822685/1‑100 (MQ=255)                                                                                         ccTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCACCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGcc       >  1:1396600/1‑100 (MQ=255)                                                                                          cgtaaGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGCCa      <  2:1244709/97‑1 (MQ=255)                                                                                           tGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGCCAt     >  2:2409736/1‑100 (MQ=255)                                                                                             aaGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGCCATca   <  1:2163307/100‑1 (MQ=255)                                                                                              aGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGTGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGCCATcac  >  2:1815471/1‑100 (MQ=255)                                                                                                 |                                                                                                 AAAGCGCGTGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGATCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGCCATCAC  >  NC_000913/1017200‑1017392

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