Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,017,296	A→C	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	98.4%	234.1 / ‑4.0	64	I335S (ATC→AGC)	ycbZ	putative peptidase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (34/29);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

AAAAGCGCGTGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGATCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGCCATCAC  >  NC_000913/1017199‑1017392                                                                                                  |                                                                                                 aaaaGCGCGTGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCtc                                                                                                <  1:89840/100‑1 (MQ=255)  aaaGCGCGTGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGctct                                                                                               <  2:1216253/100‑1 (MQ=255)  aaaGCGCGTGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGctct                                                                                               >  2:614900/1‑100 (MQ=255)  aaaGCGCGTGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGctct                                                                                               <  1:1631673/100‑1 (MQ=255)     gcgcGTGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCAt                                                                                            >  2:494202/1‑100 (MQ=255)       gcgTGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATc                                                                                           <  1:900700/99‑1 (MQ=255)        cgTGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCt                                                                                         >  1:672866/1‑100 (MQ=255)        cgTGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTCTCAATCAGGATTTGCTCCTGAAGGCTCTCATcc                                                                                          <  1:511471/99‑1 (MQ=255)          tGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGc                                                                                       <  2:1256023/100‑1 (MQ=255)          tGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGc                                                                                       <  1:215120/100‑1 (MQ=255)            gATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCTCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCAt                                                                                     <  1:576035/100‑1 (MQ=255)            gATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCAt                                                                                     <  1:1372898/100‑1 (MQ=255)                aCCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGt                                                                                 >  2:251853/1‑100 (MQ=255)                 cccGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGAGTTCGGCTTCGGTGTCACTCAGGATTTGCCCCTGTAGGCTCGCATCGTGCAt                                                                                     >  1:1462577/1‑95 (MQ=255)                   cGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCa                                                                              >  1:5621/1‑100 (MQ=255)                   cGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCa                                                                              >  2:1505565/1‑100 (MQ=255)                    ggAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAg                                                                             <  2:1700712/100‑1 (MQ=255)                        aTTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAg                                                                         <  2:1613177/100‑1 (MQ=255)                         ttCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAgg                                                                        <  2:448734/100‑1 (MQ=255)                           cAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGaa                                                                      >  1:11609/1‑100 (MQ=255)                           cAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGaa                                                                      >  2:1195773/1‑100 (MQ=255)                           cAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGaa                                                                      >  2:1195652/1‑100 (MQ=255)                           cAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGaa                                                                      <  1:1633656/100‑1 (MQ=255)                                aCCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCtt                                                                 >  2:729187/1‑100 (MQ=255)                                  cGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTcg                                                               <  2:118346/100‑1 (MQ=255)                                         gCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCAtt                                                        <  2:672867/100‑1 (MQ=255)                                          cGTTAATTTGCTCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCACCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTc                                                       >  1:57578/1‑100 (MQ=255)                                          cGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTc                                                       <  2:11609/100‑1 (MQ=255)                                          cGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTc                                                       >  2:137019/1‑100 (MQ=255)                                           gTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCa                                                      >  2:1525937/1‑100 (MQ=255)                                           gTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCa                                                      >  2:928267/1‑100 (MQ=255)                                            ttAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCAc                                                     >  2:1609602/1‑100 (MQ=255)                                             tAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACg                                                    <  2:1335523/100‑1 (MQ=255)                                                 ttGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACgctgc                                                <  2:929540/100‑1 (MQ=255)                                                   gCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGctgctg                                              <  1:1585161/100‑1 (MQ=255)                                                      cGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTgcag                                           <  2:1740632/100‑1 (MQ=255)                                                      cGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTgcag                                           <  2:1716446/100‑1 (MQ=255)                                                       gATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTgcagc                                          >  1:1419276/1‑100 (MQ=255)                                                       gATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTgcagc                                          <  1:61970/100‑1 (MQ=255)                                                        aTGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTgcagca                                         <  2:725938/100‑1 (MQ=255)                                                         tGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCAt                                        >  2:664438/1‑100 (MQ=255)                                                           cGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTa                                      <  2:1088681/100‑1 (MQ=255)                                                              tCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGt                                   <  2:1168148/100‑1 (MQ=255)                                                              tCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGATCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGt                                   <  1:144405/100‑1 (MQ=255)                                                                gCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGttt                                 <  1:1493298/100‑1 (MQ=255)                                                                gCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGttt                                 <  2:5621/100‑1 (MQ=255)                                                                gCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGttt                                 >  2:1466437/1‑100 (MQ=255)                                                                       gTTTCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAgctgct                          >  2:94853/1‑100 (MQ=255)                                                                         ttCAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCg                        >  2:437496/1‑100 (MQ=255)                                                                           cAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGcc                      <  2:542380/100‑1 (MQ=255)                                                                           cAATCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGcc                      >  2:886408/1‑100 (MQ=255)                                                                            aaTCAGGATTTGCTACTGAAGGCTCTCATCCTGCATACGTTCAGCGTGGAAAACTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCg                     >  1:439927/1‑100 (MQ=255)                                                                             aTCAGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCgg                    <  1:1459146/100‑1 (MQ=255)                                                                                aGGATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGgaga                 >  2:1279646/1‑100 (MQ=255)                                                                                 ggATTTGCTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGaa                >  2:1365442/1‑100 (MQ=255)                                                                                        cTCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCg         >  1:912846/1‑100 (MQ=255)                                                                                         tCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGc        <  2:189557/100‑1 (MQ=255)                                                                                         tCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGc        >  1:86610/1‑100 (MQ=255)                                                                                         tCCTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTACTCGCCGGAGAAGGTGTCGc        >  2:1435577/1‑100 (MQ=255)                                                                                          ccTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGcc       >  1:1893100/1‑100 (MQ=255)                                                                                          ccTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGcc       >  2:83900/1‑100 (MQ=255)                                                                                          ccTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGcc       >  2:84164/1‑100 (MQ=255)                                                                                           cTGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGTTCGCCGGAGAAGGTGTCGCCa      >  2:1725555/1‑100 (MQ=255)                                                                                            tGAAGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGCCAt     >  2:1244245/1‑100 (MQ=255)                                                                                              aaGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGCCATca   <  2:1893113/100‑1 (MQ=255)                                                                                               aGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGCCATcac  >  2:272884/1‑100 (MQ=255)                                                                                               aGGCTCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGCCATcac  >  2:1895038/1‑100 (MQ=255)                                                                                                  |                                                                                                 AAAAGCGCGTGGATGACCCGGAAATTCAATGACCGAAAGGGCGTTAATTTGCCCGATGCGTTCGCCTTCGGTTTCAATCAGGATTTGCTCCTGAAGGATCTCATCCTGCATACGTTCAGCGAGGAAACCTTCGCGCCATTCACGCTGCTGCAGCATTAAGTTTAGCTGCTCGCCGGAGAAGGTGTCGCCATCAC  >  NC_000913/1017199‑1017392

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