Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,732,082	G→A	R698C (CGT→TGT)	clpB ←	ClpB chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,732,082	0	G	A	98.0%	169.0 / NA	51	R698C (CGT→TGT)	clpB	ClpB chaperone
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (28/22);  minor base C (0/1);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAGGTAGCCACCTTCTTCATAACCGAC  >  NC_000913/2731947‑2732214                                                                                                                                        |                                                                                                                                     cGAGCACCAGCTCTTTCATGTGTGCATAAGCCAGTTCCCCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCca                                                                                                                                  <  1:401742/140‑1 (MQ=255)        cAGCTCTTTCATGTGCGCATAATCCAGTTCACCGACGCGCTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCa                                                                                                                           <  2:313011/140‑1 (MQ=255)               ttCATGTGCGCATAATCGAGTTCACCGAAGCTTTCGTGATTCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTg                                                                                                                    <  2:406636/140‑1 (MQ=255)               ttCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTg                                                                                                                    <  1:890628/140‑1 (MQ=255)                      gcgcATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGa                                                                                                             >  2:525269/1‑140 (MQ=255)                       cgcATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACATCCGTCATCCTGTACCTGCAACAGaa                                                                                                            >  1:481285/1‑140 (MQ=255)                         cATAATCCAGTTCACCGAAGCGATCCTGAATGAGATCGGTACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATg                                                                                                          <  2:384300/140‑1 (MQ=255)                           tAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGtt                                                                                                        <  1:796364/140‑1 (MQ=255)                                         gAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCGGACACCCATCATCCAGTAACTGTAAGAGAATGATGAAGACATCCGGAt                                                                                          >  1:345666/1‑140 (MQ=255)                                         gAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCCTCCAGTACCTGAAAGAGAATGTTGAAGACATCCGGAt                                                                                          >  1:358260/1‑140 (MQ=255)                                           aGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATgc                                                                                        <  1:780250/140‑1 (MQ=255)                                               ttCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGGACCTGCAACAGTATATTGAAGTCATCCGGATGCGCtt                                                                                    >  1:312873/1‑140 (MQ=255)                                                 ccTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGGACCCGCAACAGAATGTCGAAGAAATCCGGATGCGCtttt                                                                                  >  2:452526/1‑140 (MQ=255)                                                     aaTCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCa                                                                              <  1:1017339/140‑1 (MQ=255)                                                       tCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCATCAGAATGTTGAGGAGATCCGGTTGCGCTTTTTCCACt                                                                            >  1:296389/1‑140 (MQ=255)                                                           aTCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCAt                                                                        <  1:525269/140‑1 (MQ=255)                                                            tCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATc                                                                       >  1:31464/1‑140 (MQ=255)                                                            tCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATc                                                                       >  2:35936/1‑140 (MQ=255)                                                                    gAGGTTAGAGGTCATAATGACGCCCATATCACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGa                                                               <  1:283649/140‑1 (MQ=255)                                                                    gAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGa                                                               <  1:20540/140‑1 (MQ=255)                                                                     aGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGAt                                                              >  2:935354/1‑140 (MQ=255)                                                                     aGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGAt                                                              >  1:531742/1‑140 (MQ=255)                                                                     aGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGAt                                                              >  1:1020496/1‑140 (MQ=255)                                                                                  tAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGga                                                 <  1:227125/140‑1 (MQ=255)                                                                                        cgacCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGgacgacga                                           >  2:617670/1‑140 (MQ=255)                                                                                        cgacCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGgacgacga                                           >  2:565953/1‑140 (MQ=255)                                                                                              tATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCGTTTTCCACTTCAGCCATCAGGATGAAGGAATACGGACGACGACGCAcc                                     >  2:362749/1‑140 (MQ=255)                                                                                                  aCGGAAGTCGACCGTTCACCCTAGCCCGTCATTCATACCACCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCtt                                 <  2:342578/140‑1 (MQ=255)                                                                                                      aaGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGt                             >  2:788635/1‑140 (MQ=255)                                                                                                      aaGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGt                             >  2:767918/1‑140 (MQ=255)                                                                                                      aaGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGt                             >  1:24849/1‑140 (MQ=255)                                                                                                       aGTCGACCGTTCTCCTTTGCCCGTCAGTCAGACACCCGTCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTc                            <  1:452526/140‑1 (MQ=255)                                                                                                       aGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAAGATGACGGAATACGGACGACGACGCACCGCTTCGGTc                            >  1:340200/1‑140 (MQ=255)                                                                                                        gTCGACCGTTCTCCCTTGCTCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCa                           >  2:587474/1‑140 (MQ=255)                                                                                                        gTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCa                           >  1:609886/1‑140 (MQ=255)                                                                                                        gTCGACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCa                           >  1:697775/1‑140 (MQ=255)                                                                                                           gACCGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAGGt                        >  1:681147/1‑140 (MQ=255)                                                                                                             cctttCTCCCTTGCCCTTAGGTCAGACACCAACAAGCCGATCCCTGCATAAGAATGTGGAGGACTTCCGG‑TGCCCCTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAGGTAGc                     <  2:481285/137‑1 (MQ=255)                                                                                                              cGTTCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAGGTAGc                     <  2:527828/140‑1 (MQ=255)                                                                                                                ttCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGTATACGGAAGACGACACACCGCGTCGGTCAGGTAGCCa                   >  1:300043/1‑140 (MQ=255)                                                                                                                ttCTCCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGAGGACAACACACCGCATCGGTCAGGTAGCCa                   >  1:311985/1‑140 (MQ=255)                                                                                                                   tcCCTTGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAGGTAGCCACct                <  1:4174/140‑1 (MQ=255)                                                                                                                      cTTGCCCGTCAGTCAGACACCCCTCATCCAGCACCGGCAGCAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAGGTAGCCACcttct             <  2:312873/140‑1 (MQ=255)                                                                                                                       ttGCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAGGTAGCCACcttctt            >  1:922334/1‑140 (MQ=255)                                                                                                                         gCCCGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAGGTAGCCACCTTCTTCa          <  2:697775/140‑1 (MQ=255)                                                                                                                          cccGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAgg                         <  1:141157/124‑1 (MQ=255)                                                                                                                          cccGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAgg                         >  2:141157/1‑124 (MQ=255)                                                                                                                          cccGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAGGTAGCCACCTTCTTCAt         >  2:751875/1‑140 (MQ=255)                                                                                                                            cGTCAGTCAGACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAGGTAGCCACCTTCTTCATaa       <  2:24849/140‑1 (MQ=255)                                                                                                                               cagtcagACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAGGTAGCCACCTTCTTCATAACCg    >  1:890413/1‑140 (MQ=255)                                                                                                                                agtcagACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAGGTAGCCACCTTCTTCATAACCGa   <  1:35936/140‑1 (MQ=255)                                                                                                                                agtcagACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAGGTAGCCACCTTCTTCATAACCGa   <  2:890413/140‑1 (MQ=255)                                                                                                                                 gtcagACACCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAGGTAGCCACCTTCTTCATAACCGAc  <  2:609886/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     CGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGACGCACCGCTTCGGTCAGGTAGCCACCTTCTTCATAACCGAC  >  NC_000913/2731947‑2732214

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