Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	306,291	T→A	7.4%	I69N (ATT→AAT)	hslO →	molecular chaperone HslO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	306,291	0	T	A	7.4%	160.6 / 3.2	54	I69N (ATT→AAT)	hslO	molecular chaperone HslO
Reads supporting (aligned to +/- strand):  ref base T (22/28);  new base A (2/2);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

GCGGCGAGTGGGTGTCTCTCGATGACAGCTACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACC  >  NC_002947/306147‑306439                                                                                                                                                 |                                                                                                                                                      gcggcgAGTGGGTGTCTCTCGATGACAGCTACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCTGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTg                                                                                                                                                   <  1:152828/149‑1 (MQ=255) gcggcgAGTGGGTGTCTCTCGATGACAGCTACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCAGAATcag                                                                                                                                                   <  7:260036/149‑3 (MQ=255)   ggcgAGTGGGTGTCTCTCGATGACAGCTACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCa                                                                                                                                                 <  1:270239/149‑1 (MQ=255)         tGGGTGTCTCTCGATGACAGCTACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCAGAATCAGCAGgccat                                                                                                                                           <  6:134377/149‑4 (MQ=255)           ggTGTCTCTCTATGAAAGCTACGCCGCGTTGTTGGCGCGTCATGAGTACACCCAGCCGGTCCAGGTGCTGCTCTGCGAGTTGATGTCCGCCACAGCATTGCTGGTGGGTGCAATGAAGTTATATGGCCTGCTGATTATGCAGGCGAGtt                                                                                                                                         >  8:8895/1‑149 (MQ=255)               tctctcGATGACAGCTACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCggc                                                                                                                                     <  4:167036/149‑1 (MQ=255)                ctctcGATGACAGCTACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCggcc                                                                                                                                    >  7:45879/1‑149 (MQ=255)                   tcGATGACAGCTACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCggccggg                                                                                                                                 >  6:187430/1‑149 (MQ=255)                   tcGATGACAGCTACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCggccggg                                                                                                                                 >  5:92766/1‑149 (MQ=255)                     gATGACAGCTACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCggccgggcc                                                                                                                               >  3:154755/1‑149 (MQ=255)                       tGACAGCTACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGa                                                                                                                             <  4:237780/149‑1 (MQ=255)                          cAGCTACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTc                                                                                                                          <  3:104112/149‑1 (MQ=255)                              tACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGGCTGCTGATTCTGCa                                                                                                                                                 >  5:281473/1‑122 (MQ=255)                              tACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCa                                                                                                                                                 <  6:281473/122‑1 (MQ=255)                               aCGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCgctg                                                                                                                     >  6:287193/1‑149 (MQ=255)                                cgccgcGGTTTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTTCTGGTTCTGCAGGCGCGTTCGGCCGGGCCGATTCCgctg                                                                                                                     <  2:327673/148‑1 (MQ=255)                                cgccgcGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCgctgc                                                                                                                    >  3:248511/1‑149 (MQ=255)                                       gtgtTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATCGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGAtggt                                                                                                             >  5:266556/1‑149 (MQ=255)                                                cgTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTc                                                                                                    <  4:248511/149‑1 (MQ=255)                                                  tCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCa                                                                                                  <  5:187430/149‑1 (MQ=255)                                                    aTGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGt                                                                                                <  8:99584/149‑1 (MQ=255)                                                      gAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGt                                                                                                >  4:146610/1‑147 (MQ=255)                                                           cccccAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAgcgc                                                                                          >  8:184279/1‑148 (MQ=255)                                                           cccccAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGcgcg                                                                                         <  7:184279/149‑1 (MQ=255)                                                            ccccAGCCGGTCCAGGGGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCga                                                                                        <  4:312198/149‑1 (MQ=255)                                                                   cGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGt                                                                                 <  1:174163/149‑1 (MQ=255)                                                                    ggTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTg                                                                                <  3:9855/149‑1 (MQ=255)                                                                      gcctggTGCTGGTGGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGACTCGTGGc                                                                              <  6:266556/145‑1 (MQ=255)                                                                        cAGGTGCTGCTCTGCGAGTTGATGGCAGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGTCCGATTCCGCTGCTGATGGTTGAATGTTCCAGTGAGCGCGAGATTCGTGGCa                                                                             <  7:8895/148‑1 (MQ=255)                                                                         aGGTGCTGCTCGGCGAGTTGATGGCCGCCACCACCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATg                                                                           >  2:82068/1‑149 (MQ=255)                                                                          ggTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATgg                                                                          <  5:150223/149‑1 (MQ=255)                                                                             gctgctCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGccc                                                                       <  5:193990/149‑1 (MQ=255)                                                                              ctgctCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCg                                                                      <  5:287193/149‑1 (MQ=255)                                                                                                 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    >  5:108466/1‑148 (MQ=255)                                                                                                                                                aTTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCAcctga    <  6:108466/148‑1 (MQ=255)                                                                                                                                                 ttCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCAcctgacc  <  6:183267/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                      GCGGCGAGTGGGTGTCTCTCGATGACAGCTACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCACCTGACC  >  NC_002947/306147‑306439

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