Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,629,506	A→T	14.0%	P11P (CCT→CCA)	rbbA ←	ribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,629,506	0	A	T	14.0%	121.8 / 11.7	43	P11P (CCT→CCA)	rbbA	ribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
Reads supporting (aligned to +/- strand):  ref base A (22/15);  new base T (3/3);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

AATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGAGTAATTCCGG  >  NC_000913/3629365‑3629654                                                                                                                                               |                                                                                                                                                     aaTCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCTATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGcggga                                                                                                                                               <  1:240444/149‑1 (MQ=255)    cAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCTGTTGCGCAACTgg                                                                                                                                                    <  3:48081/141‑1 (MQ=255)    cAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCTGTTGCGCAACTgg                                                                                                                                                    >  4:48081/1‑141 (MQ=255)    cAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGcgggacgg                                                                                                                                            <  1:399761/149‑1 (MQ=255)    cAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGcgggacg                                                                                                                                             >  3:300934/1‑148 (MQ=255)     aaCGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGcgggacggg                                                                                                                                           >  1:227762/1‑149 (MQ=255)     aaCGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGcgggacggg                                                                                                                                           <  1:98422/149‑1 (MQ=255)        gACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGaa                                                                                                                                         >  4:169661/1‑148 (MQ=255)            acaaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAg                                                                                                                                     >  4:131245/1‑148 (MQ=255)            acaaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGt                                                                                                                                    >  2:201103/1‑149 (MQ=255)             caaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGtt                                                                                                                                   <  1:122395/149‑1 (MQ=255)              aaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGtt                                                                                                                                   >  3:370719/1‑148 (MQ=255)               aGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTc                                                                                                                                  <  3:169661/148‑1 (MQ=255)                           cGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCAt                                                                                                                      >  3:302593/1‑148 (MQ=255)                               gCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTcc                                                                                                                                 >  1:119271/1‑133 (MQ=255)                               gCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTcc                                                                                                                                 <  2:119271/133‑1 (MQ=255)                                ccGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCgg                                                                                                                 >  3:78908/1‑148 (MQ=255)                                   tCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCaa                                                                                                              <  4:302593/148‑1 (MQ=255)                                          cAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGTTCCCGTCCCGCCTGTCGCGCAACTGGcgg                                                                                                                                                 <  2:378376/107‑1 (MQ=25)                                          cAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGTTCCCGTCCCGCCTGTCGCGCAACTGGcgg                                                                                                                                                 >  1:378376/1‑107 (MQ=25)                                              cAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAg                                                                                                  <  2:227762/149‑1 (MQ=255)                                                                 gCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAg                                                                                <  4:300934/148‑1 (MQ=255)                                                                   cGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttc                                                                             >  2:253564/1‑149 (MQ=255)                                                                           cGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCCCGGAAGTTCttcattc                                                                      >  4:103612/1‑148 (MQ=255)                                                                                gTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCg                                                                 >  3:70122/1‑148 (MQ=255)                                                                                 tGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGc                                                                >  3:277390/1‑148 (MQ=255)                                                                                     atTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCAccc                                                            >  3:120822/1‑148 (MQ=255)                                                                                      ttGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCAc                                                          >  2:300149/1‑149 (MQ=255)                                                                                        gTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGc                                                        <  1:109977/149‑1 (MQ=255)                                                                                            aGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACg                                                    >  2:385263/1‑149 (MQ=255)                                                                                               gcAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGccc                                                 <  2:174767/149‑1 (MQ=255)                                                                                               gcAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGcc                                                  <  3:131245/148‑1 (MQ=255)                                                                                                          tCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGt                                      >  1:322851/1‑149 (MQ=255)                                                                                                             aTAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGttt                                    <  4:78097/148‑1 (MQ=255)                                                                                                              tAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtttt                                   >  3:49124/1‑148 (MQ=255)                                                                                                               aaTGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGa                                 <  2:228707/149‑1 (MQ=255)                                                                                                                 tGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACa                               >  1:254121/1‑149 (MQ=255)                                                                                                                          aCG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAg                       <  4:120822/148‑1 (MQ=255)                                                                                                                                cgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACTACGAGGTCGTCAGGCCACGGAAGttct                                                                            <  3:27478/90‑1 (MQ=255)                                                                                                                                cgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttct                                                                            >  4:27478/1‑90 (MQ=255)                                                                                                                                 gccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtgctg                 >  3:137633/1‑148 (MQ=255)                                                                                                                                  ccAGTTGCGCGACAGGCGTGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATgctgct               >  2:405871/1‑149 (MQ=255)                                                                                                                                   cAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATgctgct               >  3:343063/1‑148 (MQ=255)                                                                                                                                               aGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGAGTAATTCCgg  >  1:179353/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                     AATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGAGTAATTCCGG  >  NC_000913/3629365‑3629654

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