Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	306,288	T→A	8.0%	L68Q (CTG→CAG)	hslO →	molecular chaperone HslO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	306,288	0	T	A	8.0%	143.8 / 3.4	50	L68Q (CTG→CAG)	hslO	molecular chaperone HslO
Reads supporting (aligned to +/- strand):  ref base T (20/26);  new base A (2/2);  total (22/28)
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.54e-01

GCGGCGAGTGGGTGTCTCTCGATGACAGCTACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCAC  >  NC_002947/306147‑306433                                                                                                                                              |                                                                                                                                                   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)                                                                                                                                              |                                                                                                                                                   GCGGCGAGTGGGTGTCTCTCGATGACAGCTACGCCGCGGTGTTGGCGCGTCATGAGTACCCCCAGCCGGTCCAGGTGCTGCTCGGCGAGTTGATGGCCGCCACCGCCTTGCTGGTGGGTGCAATGAAGTTCGATGGCCTGCTGATTCTGCAGGCGCGTTCGGCCGGGCCGATTCCGCTGCTGATGGTGGAATGTTCCAGTGAGCGCGAGATTCGTGGCATGGCCCGCTACGAAGC‑GGACCAGATTCCGGCTGGTGCCACCCTGTCGCAGCTGATGCCCGACGGCCAC  >  NC_002947/306147‑306433

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