Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,867,164	C→T	A130T (GCG→ACG)	rpoS ←	RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,867,164	0	C	T	93.0%	217.0 / 6.3	71	A130T (GCG→ACG)	rpoS	RNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (27/39);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

TTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTT  >  NC_000913/2867068‑2867260                                                                                                 |                                                                                                 tttGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACcgtgcg                                                                                               <  1:787700/100‑1 (MQ=255) tttGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACcgtgcg                                                                                               <  1:334549/100‑1 (MQ=255)    ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGACGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGAt                                                                                            >  2:1554416/1‑100 (MQ=255)    ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGAt                                                                                            <  2:534356/100‑1 (MQ=255)       tCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAg                                                                                         <  2:611351/100‑1 (MQ=255)        cATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGc                                                                                        >  2:863406/1‑100 (MQ=255)          tAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGccc                                                                                      >  2:459715/1‑100 (MQ=255)          tAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGccc                                                                                      >  2:459870/1‑100 (MQ=255)            aTCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCa                                                                                    <  1:176533/100‑1 (MQ=255)            aTCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCa                                                                                    >  2:708148/1‑100 (MQ=255)              cGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAgg                                                                                  <  2:343394/100‑1 (MQ=255)              cGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAgg                                                                                  <  1:1740085/100‑1 (MQ=255)               gCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGt                                                                                 <  1:1001240/100‑1 (MQ=255)               gCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGt                                                                                 >  2:982809/1‑100 (MQ=255)               gCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGt                                                                                 <  2:989375/100‑1 (MQ=255)                cccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGtt                                                                                <  1:908009/100‑1 (MQ=255)                cccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGtt                                                                                >  1:1289865/1‑100 (MQ=255)                  cGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGc                                                                              <  2:1668916/100‑1 (MQ=255)                  cGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGc                                                                              <  1:1818119/100‑1 (MQ=255)                   gTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGcc                                                                             >  2:862090/1‑100 (MQ=255)                     tCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCct                                                                           >  2:216615/1‑100 (MQ=255)                      cAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCctc                                                                          >  1:622428/1‑100 (MQ=255)                         tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCctct                                                                         <  1:1600036/98‑1 (MQ=255)                         tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTc                                                                       <  1:1411704/100‑1 (MQ=255)                         tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTc                                                                       <  2:816533/100‑1 (MQ=255)                          cGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCg                                                                      <  1:459713/100‑1 (MQ=255)                          cGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCg                                                                      <  1:459868/100‑1 (MQ=255)                          cGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCg                                                                      <  1:708147/100‑1 (MQ=255)                           gTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGa                                                                     <  2:1204490/100‑1 (MQ=255)                              tGGCGAATCCACCAGGTTGCGTATGTTGCGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATaa                                                                  <  2:715381/100‑1 (MQ=255)                               ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAg                                                                 <  2:323835/100‑1 (MQ=255)                                 cGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGt                                                               >  1:1387623/1‑100 (MQ=255)                                    aTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCa                                                            <  1:1108488/100‑1 (MQ=255)                                      ccaccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGc                                                          >  1:1083995/1‑100 (MQ=255)                                          caGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACg                                                      <  1:636338/100‑1 (MQ=255)                                           aGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGc                                                     <  2:571056/100‑1 (MQ=255)                                                 cGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGAcc                                               >  2:249401/1‑100 (MQ=255)                                                 cGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGAcc                                               >  2:1845982/1‑100 (MQ=255)                                                    aTGTTGAGAGGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACg                                            >  1:1473439/1‑100 (MQ=255)                                                     tGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                           >  2:698917/1‑100 (MQ=255)                                                      gTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGAt                                          <  1:1112555/100‑1 (MQ=255)                                                      gTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGAt                                          <  2:1199404/100‑1 (MQ=255)                                                      gTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGAt                                          >  2:1479943/1‑100 (MQ=255)                                                      gTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGAt                                          >  2:1478391/1‑100 (MQ=255)                                                        tGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGAtt                                         >  1:385172/1‑99 (MQ=255)                                                        tGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTg                                        <  1:20423/100‑1 (MQ=255)                                                         gagaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGc                                       <  1:703331/100‑1 (MQ=255)                                                           gaAGCGGAGACCACGTTCCTGGTCAAAGTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGcc                                      <  2:357524/99‑1 (MQ=255)                                                           gaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCa                                     <  2:1904123/100‑1 (MQ=255)                                                           gaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCa                                     <  1:541320/100‑1 (MQ=255)                                                           gaAGCGGAAACCACGTTCCGGGTCAAACTTCACTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCa                                     >  2:1307573/1‑100 (MQ=255)                                                             aGCGGAAACCGCGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATa                                   <  1:910334/100‑1 (MQ=255)                                                                ggAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcg                                >  1:543657/1‑100 (MQ=255)                                                                  aaaCCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggc                              >  2:1375881/1‑100 (MQ=255)                                                                  aaaCCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggc                              <  1:167081/100‑1 (MQ=255)                                                                    aCCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcgg                            >  2:1007255/1‑100 (MQ=255)                                                                    aCCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcgg                            <  1:1478383/100‑1 (MQ=255)                                                                    aCCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcgg                            <  1:1479935/100‑1 (MQ=255)                                                                     ccACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCggg                           >  1:1119540/1‑100 (MQ=255)                                                                     ccACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGGCCATCAACGCCAGACCACGATTTCCATAAcggcg                             >  1:1620929/1‑98 (MQ=255)                                                                              gggTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa                  <  1:1381931/100‑1 (MQ=255)                                                                               ggTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTac                 >  2:137004/1‑100 (MQ=255)                                                                                gTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacc                <  1:944074/100‑1 (MQ=255)                                                                                gTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacc                <  2:781692/100‑1 (MQ=255)                                                                                gTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacc                <  2:1083996/100‑1 (MQ=255)                                                                                  cAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccac              <  2:961673/100‑1 (MQ=255)                                                                                  cAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccac              >  1:754098/1‑100 (MQ=255)                                                                                      cTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACTCCAGACCACGATTGACCTAAAGGCGGGCAATTTTTAacaaca            >  1:182892/1‑98 (MQ=255)                                                                                           ctACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCa      <  1:195901/99‑1 (MQ=255)                                                                                            tACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAg    <  1:884987/100‑1 (MQ=255)                                                                                            tACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACAACCAGACGCAAg    <  2:1276267/100‑1 (MQ=255)                                                                                             aCCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGt   <  1:702981/100‑1 (MQ=255)                                                                                              ccGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGtt  >  2:582192/1‑100 (MQ=255)                                                                                                 |                                                                                                 TTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTT  >  NC_000913/2867068‑2867260

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