Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	306,378	A→T	8.9%	E98V (GAA→GTA)	hslO →	molecular chaperone HslO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	306,378	0	A	T	8.9%	123.8 / 3.4	45	E98V (GAA→GTA)	hslO	molecular chaperone HslO
Reads supporting (aligned to +/- strand):  ref base A (22/19);  new base T (2/2);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01

CCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCAGGGGACCGTCGACCTCGATGGCGCCAACCTGTCTGAATGTTTCACCAACTAC  >  NC_002947/306243‑306526                                                                                                                                        |                                                                                                                                                      ccgccACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGAtt                                                                                                                                         >  2:269133/1‑149 (MQ=255) ccgccACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGCCATGGCCCGCTACGAAGC‑GGACCAGAtt                                                                                                                                         >  5:270178/1‑149 (MQ=255)     cACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCgg                                                                                                                                     >  5:183267/1‑149 (MQ=255)      aCCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGc                                                                                                                                    <  3:146610/149‑1 (MQ=255)       ccGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCt                                                                                                                                   <  4:154755/149‑1 (MQ=255)         gCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTgg                                                                                                                                 <  6:92766/149‑1 (MQ=255)               cTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGGCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCa                                                                                                                            <  1:269133/148‑1 (MQ=255)                tggtggGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCAcc                                                                                                                          >  4:128507/1‑149 (MQ=255)                tggtggGTGCAATGAAGTTCGATGGCCTCCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCAcc                                                                                                                          >  1:184295/1‑149 (MQ=255)                  gtgggtgCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCt                                                                                                                        >  3:97189/1‑149 (MQ=255)                     ggtgAAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGCCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGCTTCCGGCTGGTCCCACCCTGTc                                                                                                                     >  6:43192/1‑149 (MQ=255)                         cAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCa                                                                                                                  <  8:45879/148‑1 (MQ=255)                            tGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGTCCGCTTCGTAGC‑gg                                                                                                                                                 >  6:189400/1‑114 (MQ=255)                            tGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGGTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGTCCGCTTCGTAGC‑gg                                                                                                                                                 <  5:189400/114‑1 (MQ=255)                                    acgcgCCTGCAGAATCAGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGa                                                                                                      >  4:10444/5‑149 (MQ=255)                                      gcgCCTGCAGAATCAGCAGGCGCGTTCGGCCCGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTTAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACg                                                                                                    >  8:69382/3‑149 (MQ=255)                                        gCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGc                                                                                                  <  6:270178/149‑1 (MQ=255)                                          ctgctgATTCTGCAGGCGCGTTCGGCCGGGCCGATTCAGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCAACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCa                                                                                                >  5:14822/1‑149 (MQ=255)                                             ctgATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGCTGGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCAc                                                                                               <  1:82068/148‑1 (MQ=255)                                                aTTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCAcctga                                                                                           <  6:108466/148‑1 (MQ=255)                                                aTTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCAcctga                                                                                           >  5:108466/1‑148 (MQ=255)                                                 ttCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCAcctgacc                                                                                         <  6:183267/149‑1 (MQ=255)                                                    tGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCAcctgaccct                                                                                       >  5:318392/1‑148 (MQ=255)                                                       aGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCg                                                                                                                                      >  8:91230/1‑98 (MQ=255)                                                       aGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCg                                                                                                                                      <  7:91230/98‑1 (MQ=255)                                                             gTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGAc                                                                             <  6:318392/149‑1 (MQ=255)                                                              ttCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGAcc                                                                            <  3:10444/149‑1 (MQ=255)                                                                 ggccgggccgATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCgg                                                                         >  2:116905/1‑149 (MQ=255)                                                                          gATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGc                                                                 >  3:137000/1‑148 (MQ=255)                                                                                     tgATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCAg                                                     >  2:118569/1‑149 (MQ=255)                                                                                             gAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCAGGGGACCGt                                             <  4:97189/149‑1 (MQ=255)                                                                                                  ttCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCAGGGGACCGTCGAcc                                        >  4:233759/1‑149 (MQ=255)                                                                                                  ttCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCAGGGGACCGTCGAc                                         <  2:184295/148‑1 (MQ=255)                                                                                                   tCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCAGGGGACCGTCGACCt                                       >  5:161396/1‑149 (MQ=255)                                                                                                        tGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCa                                                               <  2:310578/120‑1 (MQ=255)                                                                                                        tGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCa                                                               >  1:310578/1‑120 (MQ=255)                                                                                                               gagaTTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCAGGGGACCGTCGACCTCGATGGCGCCaa                           >  1:141790/1‑149 (MQ=255)                                                                                                                      gTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCAGGGGACCGTCGACCTCGATGGCGCCAACCTGTCt                    >  4:63268/1‑149 (MQ=255)                                                                                                                       gaatctGGTCCGCTTCGTAGC‑GGGCCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCAGGGGACCGTCGAcc                                        >  2:186513/6‑128 (MQ=255)                                                                                                                       gaatctGGTCCGCTTCGTAGC‑GGGCCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCAGGGGACCGTCGAcc                                        <  1:186513/123‑1 (MQ=255)                                                                                                                         gCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCAGGGGACCGTCGACCTCGATGGCGCCAACCTGTCTGaa                 <  6:161396/149‑1 (MQ=255)                                                                                                                          cctGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCGGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCTTGGGACCGTCGACCTCGATGGCGCCAACCTGTCTGAAt                <  6:14822/147‑1 (MQ=255)                                                                                                                              gCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCAGGGGACCGTCGACCTCGATGGCGCCAACCTGTCTGAATGttt            <  3:63268/149‑1 (MQ=255)                                                                                                                              gCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCAGGGGACCGTCGACCTCGATGGCGCCAACCTGTCTGAATGtt             <  2:141790/148‑1 (MQ=255)                                                                                                                                      cGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCAGGGGACCGTCGACCTCGATGGCGCCAACCTGTCTGAATGTTTCACCAac     >  1:143852/1‑148 (MQ=255)                                                                                                                                        aaGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCAGGGGACCGTCGACCTCGATGGCGCCAACCTGTCTGAATGTTTCACCAactac  >  1:57029/1‑149 (MQ=255)                                                                                                                                        |                                                                                                                                                      CCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACCCTGACCATCGACCCGGTCAAGGGCCAGCGCTACCAGGGGACCGTCGACCTCGATGGCGCCAACCTGTCTGAATGTTTCACCAACTAC  >  NC_002947/306243‑306526

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑