Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,018,260	(A)8→7	100%	coding (315/444 nt)	fliJ →	flagellar protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,018,253	0	A	.	98.5%	290.3 / ‑2.1	65	coding (308/444 nt)	fliJ	flagellar protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (31/33);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

AGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGG  >  NC_000913/2018154‑2018348                                                                                                    |                                                                                                agcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCgaga                                                                                                 >  1:651504/1‑100 (MQ=255) agcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCgaga                                                                                                 >  2:877420/1‑100 (MQ=255)  gcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaa                                                                                                >  2:1718717/1‑100 (MQ=255)    agTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaa                                                                                              >  1:1319225/1‑100 (MQ=255)        tATCCAGACGCTGGAAAAAGCCATTACTCAACATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaac                                                                                          >  1:1344122/1‑99 (MQ=255)         aTCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaaca                                                                                         >  1:22321/1‑98 (MQ=255)          tCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCCGTAGACGCAAAAAGTTGACATTGAACTGAACAGTTGGCGAGaaaaaaacaa                                                                                        >  1:440137/1‑97 (MQ=255)          tCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAACAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaa                                                                                        >  2:1190946/1‑97 (MQ=255)            cAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACg                                                                                     >  2:680239/1‑100 (MQ=255)            cAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACg                                                                                     >  2:1720976/1‑100 (MQ=255)            cAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACg                                                                                     >  1:2116229/1‑100 (MQ=255)                  cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCa                                                                               >  1:1708137/1‑100 (MQ=255)                  cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCa                                                                               >  1:459368/1‑100 (MQ=255)                   tGGAAAAAGCCATTACTGATCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTTAACATTTGGCGAG‑AAAAAAACAACg                                                                                     <  2:850433/93‑1 (MQ=255)                    ggAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAgg                                                                             >  2:125181/1‑100 (MQ=255)                     gAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGc                                                                            <  2:2189285/100‑1 (MQ=255)                      aaaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGcc                                                                           <  1:585879/100‑1 (MQ=255)                        aaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTg                                                                         >  1:411141/1‑100 (MQ=255)                           gCCATTACTCAGCATCGCCAGCATCTTAATCATTGGACGCAGAAAGTTTACATTGCTCTGAACAGTTTGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGc                                                                       <  2:1331522/99‑1 (MQ=255)                           gCCATTACTCAACATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCa                                                                      <  2:1344116/100‑1 (MQ=255)                            ccATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAg                                                                     <  1:1056878/100‑1 (MQ=255)                             cATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGCTGGCGAGAAAAGAAACAACGTTTGCAGGCCTGGCAg                                                                     >  2:1884338/1‑100 (MQ=255)                               ttACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGaca                                                                  <  2:757751/100‑1 (MQ=255)                                tACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGacac                                                                 <  2:1658998/100‑1 (MQ=255)                                  cTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACt                                                                <  1:461579/99‑1 (MQ=255)                                   tCAACATCTCCAGCAACTTAATCAGTGGACGCAGAAAGTTTACATTGCCCTGTACAGTTGGCTAG‑AAAAAAACTACTTTTGCAGGCCTGGCAGACACTGc                                                              <  2:1092639/100‑1 (MQ=255)                                     aGCATCGCCAGCAACATAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAg                                                            <  1:2209115/100‑1 (MQ=255)                                      gCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAgg                                                           <  1:877421/100‑1 (MQ=255)                                           gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACgg                                                      >  2:1651273/1‑100 (MQ=255)                                           gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACgg                                                      <  1:1839221/100‑1 (MQ=255)                                           gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACgg                                                      <  1:1809888/100‑1 (MQ=255)                                           gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACg                                                       <  1:67046/99‑1 (MQ=255)                                            ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGGCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  2:2060627/100‑1 (MQ=255)                                                cAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATc                                                 <  2:173117/100‑1 (MQ=255)                                                 aaCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATcc                                                <  2:1425965/100‑1 (MQ=255)                                                  ccTTAATCAGTGGACGCAGAAAGTTGACATTACCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGAAGACACTGCAGGAACGGCAATCCa                                               >  2:924703/2‑100 (MQ=255)                                                  aCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAAGAACGGCACTCCa                                               >  2:925066/1‑100 (MQ=255)                                                    ttAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcg                                             <  1:479522/100‑1 (MQ=255)                                                          aGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCa                                       >  1:1484793/1‑100 (MQ=255)                                                            tGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACt                                     >  1:84985/1‑100 (MQ=255)                                                              gACGCAGAAAGTTGACATTGCCCTGAACAGTTGGAGAG‑AAAAAAACACCGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGCCGGCACTGc                                   >  1:386455/1‑100 (MQ=255)                                                              gACGCAGAAAGTTGACATTGCCATAAACAGTAGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTg                                    <  1:1110526/99‑1 (MQ=255)                                                               aCGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCt                                  <  1:5515/100‑1 (MQ=255)                                                                 gCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTg                                >  2:2213865/1‑100 (MQ=255)                                                                  cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               <  1:70501/100‑1 (MQ=255)                                                                  cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               >  1:1665731/1‑100 (MQ=255)                                                                  cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               >  1:1665369/1‑100 (MQ=255)                                                                   aGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               <  1:1008554/99‑1 (MQ=255)                                                                   aGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCa                              <  1:2231506/100‑1 (MQ=255)                                                                     aaaGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCCATCCACGGCGGCACTGCTTGCAGa                            >  2:417070/1‑100 (MQ=255)                                                                        gTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaaa                         >  1:664617/1‑100 (MQ=255)                                                                        gTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCAACGGCGACACTACTTGCAGaaaa                         >  1:187814/1‑100 (MQ=255)                                                                         ttGACATTGCCCTGAACAGTTGGAGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACTGCGGCACTGCTTGCAGAAAAc                        <  2:183023/100‑1 (MQ=255)                                                                          tGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAGCAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAcc                       <  2:10649/100‑1 (MQ=255)                                                                              aTTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGcc                    >  1:1460469/1‑99 (MQ=255)                                                                              aTTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCt                   >  2:697027/1‑100 (MQ=255)                                                                               ttGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTc                  >  1:1932872/1‑100 (MQ=255)                                                                               ttGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTc                  <  2:2116215/100‑1 (MQ=255)                                                                               ttGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTc                  >  1:425766/1‑100 (MQ=255)                                                                                 gCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGa                >  1:2106514/1‑100 (MQ=255)                                                                                 gCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGCGAGACACTGCAGGAACGGCAATCCGTGGCGGTACTTCTTGTAGAAAACCGCCGCGa                <  2:303965/100‑1 (MQ=255)                                                                                   ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              >  2:333165/1‑100 (MQ=255)                                                                                   ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  2:1034744/100‑1 (MQ=255)                                                                                    cTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa             <  2:1734300/100‑1 (MQ=255)                                                                                    cTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa             >  2:576215/1‑100 (MQ=255)                                                                                     tGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGAt               <  1:230984/97‑1 (MQ=255)                                                                                      gAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAg            >  2:525431/1‑99 (MQ=255)                                                                                      gAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGACAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGa           >  2:1611325/1‑100 (MQ=255)                                                                                        aCAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaaa         >  2:187672/1‑100 (MQ=255)                                                                                              ggCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatg   <  1:2205993/100‑1 (MQ=255)                                                                                               gCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatgg  <  1:1698587/100‑1 (MQ=255)                                                                                               gCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatgg  >  2:2207683/1‑100 (MQ=255)                                                                                               gCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatgg  >  1:1966754/1‑100 (MQ=255)                                                                                                    |                                                                                                AGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGG  >  NC_000913/2018154‑2018348

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