Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,849,104	A→G	14.2%	intergenic (+23/‑248)	rlmI → / → ilvD	ribosomal RNA large subunit methyltransferase I/dihydroxy‑acid dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,849,104	0	A	G	14.2%	113.8 / 11.0	51	intergenic (+23/‑248)	rlmI/ilvD	ribosomal RNA large subunit methyltransferase I/dihydroxy‑acid dehydratase
Reads supporting (aligned to +/- strand):  ref base A (21/21);  new base G (2/5);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

TGGACCGCAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCACGCGGCGATCCCGCGCG  >  NC_002947/5848960‑5849249                                                                                                                                                   |                                                                                                                                                  tGGACCGCAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGAATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGACGCCCAGCGGc                                                                                                                                              >  2:364196/1‑150 (MQ=255)     ccGCAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAATCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCTGGGCGGCCCtt                                                                                                                                          <  1:252890/150‑1 (MQ=255)       gCAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGTTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTCCTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCTGGGCGGCCCtttg                                                                                                                                        <  2:601605/150‑2 (MQ=255)       gCAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCTGGGCGGCCCtttt                                                                                                                                        <  2:605059/150‑2 (MQ=255)       gCAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCTGGGCGGCCCtttt                                                                                                                                        <  1:503132/150‑2 (MQ=255)                    gctCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATccc                                                                                                                           <  1:465349/150‑1 (MQ=255)                            gagGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCaa                                                                                                                   <  2:163768/150‑1 (MQ=255)                             agGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAAc                                                                                                                  <  1:658918/150‑1 (MQ=255)                                             gATCACCCGGTACACCCGGCCATCGCCGAAACCCGCGACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAATCGCCTCCCGGCGGCCCTTTTCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCAtt                                                                                                  >  1:529568/1‑150 (MQ=255)                                                  cccGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTc                                                                                             >  2:579632/1‑150 (MQ=255)                                                   ccGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGc                                                                                                                                              <  2:52208/100‑1 (MQ=255)                                                   ccGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGc                                                                                                                                              >  1:52208/1‑100 (MQ=255)                                                   ccGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCa                                                                                            >  1:784020/1‑150 (MQ=255)                                                       tACACCCGGCCAGCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACcgc                                                                                        <  1:579632/150‑1 (MQ=255)                                                          acCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGcc                                                                                     >  2:588815/1‑150 (MQ=255)                                                          acCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGcc                                                                                     >  1:451073/1‑150 (MQ=255)                                                          acCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGcc                                                                                     >  1:760039/1‑150 (MQ=255)                                                            ccGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAg                                                                                   >  1:732123/1‑150 (MQ=255)                                                                ccATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGgt                                                                               <  2:451073/150‑1 (MQ=255)                                                                   tCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTa                                                                            >  1:588831/1‑150 (MQ=255)                                                                    cGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAg                                                                           <  2:760039/150‑1 (MQ=255)                                                                     gCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTTAT‑‑ATCTGCCAAAAAGGGCCGCCCCGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTATAGa                                                                          >  2:412036/1‑150 (MQ=255)                                                                     gCCGAAACCCGCTACATCAAAAGAATAACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCGGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGa                                                                          <  2:529568/150‑1 (MQ=255)                                                                         aaaCCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACTCCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCg                                                                      <  1:484893/150‑1 (MQ=255)                                                                              cGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTa                                                                 >  1:524741/1‑150 (MQ=255)                                                                               ccTACATCAAAACCATCACCTGCCGGTTACTGTCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTAt                                                                <  1:328651/149‑1 (MQ=255)                                                                                cTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTAtt                                                               >  1:677413/1‑150 (MQ=255)                                                                                     tCAAAAGCATCACCTGCCGGTTACTGCCAAACCGCTCATCAATGTGCCACAA‑‑CCATTCCTATCGGCCCTCTCCGTCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCg                                                          >  1:216554/1‑150 (MQ=255)                                                                                     tCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCg                                                                                 >  2:424198/1‑127 (MQ=255)                                                                                     tCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCg                                                                                 <  1:424198/127‑1 (MQ=255)                                                                                       aaaaGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTa                                                                 >  1:273923/1‑141 (MQ=255)                                                                                       aaaaGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTa                                                                 <  2:273923/141‑1 (MQ=255)                                                                                        aaaGCATCCCCTGCCGGTTACGGCCAAACAGCTAAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCa                                                       <  2:792132/150‑1 (MQ=255)                                                                                             aTCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCAt                                                  >  2:266467/1‑150 (MQ=255)                                                                                                     ccGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCAcc                                                                                          >  2:259829/1‑102 (MQ=255)                                                                                                     ccGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCAcc                                                                                          <  1:259829/102‑1 (MQ=255)                                                                                                      cgatctcTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCAGCACCATCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCcggcgg                                         <  1:48571/144‑1 (MQ=255)                                                                                                            cTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTAt                                   <  2:216554/150‑1 (MQ=255)                                                                                                             tGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGACCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCAATCCCCCAGCCAGCGTTGTAGAATCGTCCTATTCAGCGCCCGTCATCCCCGGCgggagatcg                                  >  2:48571/1‑144 (MQ=255)                                                                                                              tccAAACAGCTCAT‑‑CTCTGCCAAAAAGGGCCGCCCCGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGa                                 <  1:627908/149‑1 (MQ=255)                                                                                                               ccAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTAt                                   <  2:697682/147‑1 (MQ=255)                                                                                                               ccAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTAt                                   >  1:697682/1‑147 (MQ=255)                                                                                                               ccAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTAt                                   <  2:147229/147‑1 (MQ=255)                                                                                                               ccAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTAt                                   >  1:147229/1‑147 (MQ=255)                                                                                                                      gCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGt                         <  1:588815/150‑1 (MQ=255)                                                                                                                           t‑‑atCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGc                    <  1:266467/150‑1 (MQ=255)                                                                                                                                 tGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGctc                             >  2:197297/1‑137 (MQ=255)                                                                                                                                 tGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGctc                             <  1:197297/137‑1 (MQ=255)                                                                                                                                      gaaaGGGCCGCTGGGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCACGCGGCGa           >  1:549251/2‑150 (MQ=255)                                                                                                                                          gggCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCACGCGGCGATccc       <  2:524741/150‑1 (MQ=255)                                                                                                                                               gccCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCACGCGGCGATCCcgcgcg  >  1:104578/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TGGACCGCAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGAT‑‑ATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCACGCGGCGATCCCGCGCG  >  NC_002947/5848960‑5849249

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