Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,092,469	T→A	10.1%	intergenic (+72/+19)	deaD → / ← ygiD	ATP‑dependent DEAD‑box RNA helicase DeaD/4,5‑DOPA‑dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,092,469	0	T	A	10.1%	136.9 / 8.5	50	intergenic (+72/+19)	deaD/ygiD	ATP‑dependent DEAD‑box RNA helicase DeaD/4,5‑DOPA‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (22/21);  new base A (2/3);  total (25/25)
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.41e-01

CTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATGCTCATCACTGGGATGGTTACGCA  >  NC_002947/2092326‑2092618                                                                                                                                                 |                                                                                                                                                      ctgctgTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCcg                                                                                                                                                 >  1:1488217/1‑150 (MQ=255)       tCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAaga                                                                                                                                           >  1:1010578/1‑150 (MQ=255)        ccAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAagag                                                                                                                                          <  1:1428329/150‑1 (MQ=255)         cAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAgg                                                                                                                                         >  1:941088/1‑150 (MQ=255)         cAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAgg                                                                                                                                         >  1:378149/1‑150 (MQ=255)              tCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCtt                                                                                                                                    >  2:248116/1‑150 (MQ=255)                      acTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGaa                                                                                                                            >  2:554331/1‑150 (MQ=255)                        tCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCAAGCCCGCGaaagg                                                                                                                                          >  2:1062968/1‑131 (MQ=255)                        tCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCAAGCCCGCGaaagg                                                                                                                                          <  1:1062968/134‑4 (MQ=255)                          gcgTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGa                                                                                                                        <  1:1639390/150‑1 (MQ=255)                          gcgTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGa                                                                                                                        <  1:1158336/150‑1 (MQ=255)                           cgTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGAt                                                                                                                       <  1:1362329/150‑1 (MQ=255)                             tggctggcAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAg                                                                                                                     <  2:378149/150‑1 (MQ=255)                                 tggcAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGt                                                                                                                 <  1:554331/150‑1 (MQ=255)                                   gcAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTcc                                                                                                               <  2:1010578/150‑1 (MQ=255)                                       gcagcTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATcc                                                                                                           <  1:248116/150‑1 (MQ=255)                                          gcTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCa                                                                                                        <  2:1236654/150‑1 (MQ=255)                                            tGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAgc                                                                                                      >  1:1595324/1‑150 (MQ=255)                                            tGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCCTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAg                                                                                                                     <  1:618642/135‑1 (MQ=255)                                            tGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCCTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAg                                                                                                                     >  2:618642/1‑135 (MQ=255)                                             gAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATAACAGGTCTGGCTCTGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAgcg                                                                                                     >  2:1379594/1‑150 (MQ=255)                                                                     ttGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGca                                                                             <  2:1340513/150‑1 (MQ=255)                                                                         tCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAg                                                                                              <  1:873050/129‑1 (MQ=255)                                                                         tCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAg                                                                                              >  2:873050/1‑129 (MQ=255)                                                                         tCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTAAGTCGCCAACGCGCCCAGGGTAAACCCCTGATGcaccac                                                                         >  1:84703/1‑150 (MQ=255)                                                                                  aaCGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACtt                                                                <  2:1488217/150‑1 (MQ=255)                                                                                         aaaTCCCCGACCTAAGCCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGaa                                                                   <  2:243392/140‑1 (MQ=255)                                                                                         aaaTCCCCGACCTAAGCCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGaa                                                                   >  1:243392/1‑140 (MQ=255)                                                                                           aTCCCCGACCTAAGCCGGGGATTTTTTTTACCTGTACCGGCCCTCTTCGCGGGCAAGCCCGCGaaa                                                                                                                                            <  1:625159/66‑2 (MQ=9)                                                                                                  gCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGATGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAgg                                                >  1:866420/1‑150 (MQ=255)                                                                                                    ttAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAgtgg                                              >  1:651553/1‑147 (MQ=255)                                                                                                       ggTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCaaa                                           <  2:941088/150‑1 (MQ=255)                                                                                                             ggATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCAGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAg                                     >  1:1574134/1‑150 (MQ=255)                                                                                                                ttttttttACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCgg                                  >  1:1202756/1‑150 (MQ=255)                                                                                                                   tttttACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGgcag                               <  2:1028856/150‑1 (MQ=255)                                                                                                                       tACCTGTACCGGCCCT‑TTCGCGGGCTAGCCCGCGAAGAGGCTCTTCAGTCGAACCGATAGATGTCCATCCCAAGCGCGCCCAGGGTAAACCCATGATGCACCACGCCGCACTTGCTCACCGCCCCAAGGGCACAGCACAGCGGCAGCaaa                           >  2:678933/1‑150 (MQ=255)                                                                                                                         ccTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTaa                                                                                         <  1:808539/86‑1 (MQ=255)                                                                                                                         ccTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTaa                                                                                         >  2:808539/1‑86 (MQ=255)                                                                                                                         ccTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATg                         <  2:1595324/150‑1 (MQ=255)                                                                                                                          cTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATGc                        >  1:1168976/1‑150 (MQ=255)                                                                                                                               ccGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAggg                                               <  1:377178/122‑1 (MQ=255)                                                                                                                               ccGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAggg                                               >  2:377178/1‑122 (MQ=255)                                                                                                                                  gCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCATTCAGTCGAACCGATAGATGt                                                                                                                 <  2:209875/53‑1 (MQ=255)                                                                                                                                  gCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCATTCAGTCGAACCGATAGATGt                                                                                                                 >  1:209875/1‑53 (MQ=255)                                                                                                                                    ccT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATGCtc                      <  2:1046426/142‑1 (MQ=255)                                                                                                                                    ccT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATGCtc                      >  1:1046426/1‑142 (MQ=255)                                                                                                                                     ctcttCGCGGGCAAGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAg                                                 >  1:828943/4‑115 (MQ=255)                                                                                                                                     ctcttCGCGGGCAAGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAg                                                 <  2:828943/112‑1 (MQ=255)                                                                                                                                     cT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATGCTCATCACTggg             >  1:1030982/1‑150 (MQ=255)                                                                                                                                             gggCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATGCTCATCACTGGGATGGTTa      <  2:866420/150‑1 (MQ=255)                                                                                                                                                 ttGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATGCTCATCACTGGGATGGTTACGCa  >  2:49685/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                      CTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCT‑TTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATGCTCATCACTGGGATGGTTACGCA  >  NC_002947/2092326‑2092618

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