Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,674,716 A→G 13.8% intergenic (+27/+114) rhlE‑II → / ← yceI ATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,674,7160AG13.8% 69.3 / 3.7 29intergenic (+27/+114)rhlE‑II/yceIATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor
Reads supporting (aligned to +/- strand):  ref base A (11/14);  new base G (2/2);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

ATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGA‑‑‑TCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAG  >  NC_002947/5674579‑5674865
                                                                                                                                         |                                                                                                                                                        
aTCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGCGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGa                                                                                                                                              <  1:64080/150‑1 (MQ=255)
aTCATTCGAAAGGAGTCCAAGCTGGACCCTTTCCCTACGCCTCATCAACTGGATGACCTGCCAAGCCGCCCACGCGTTGAGCGCCCAGCGCTGGTGACCCGCAAGGGCTGTTAAAATTCCCGAGCGGGCTGCTTCGCAGCCCTATCACGa                                                                                                                                              <  1:662032/150‑1 (MQ=255)
                   aGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCag                                                                                                                           <  2:240392/150‑1 (MQ=255)
                            gTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCt                                                                                                                  >  1:443772/1‑150 (MQ=255)
                             tttCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCtt                                                                                                                 <  2:23242/150‑1 (MQ=255)
                                              aaCTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGc                                                                                                                                      >  2:185355/1‑112 (MQ=255)
                                              aaCTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGc                                                                                                                                      <  1:185355/112‑1 (MQ=255)
                                               aCTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCGGCCCTATCACGACACAAgg                                                                                                                                       >  2:691492/1‑110 (MQ=255)
                                               aCTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCGGCCCTATCACGACACAAgg                                                                                                                                       <  1:691492/110‑1 (MQ=255)
                                               aCTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGA‑‑‑TCTCTACAGGAGCGGc                                                                                            <  2:327553/150‑1 (MQ=255)
                                                  ggATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGA‑‑‑TCTCTACAGGAGCGGCCtt                                                                                         <  2:581673/150‑1 (MQ=255)
                                                                     ccACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGA‑‑‑TCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCa                                                                      >  2:193812/1‑150 (MQ=255)
                                                                            gTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGA‑‑‑TCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGc                                                               >  1:314244/1‑150 (MQ=255)
                                                                                gcgcCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCaga                                                                                                                          <  2:269489/90‑1 (MQ=255)
                                                                                gcgcCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCaga                                                                                                                          >  1:269489/1‑90 (MQ=255)
                                                                                   cccAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCt                                                                                                                  >  2:126500/1‑95 (MQ=255)
                                                                                   cccAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCt                                                                                                                  <  1:126500/95‑1 (MQ=255)
                                                                                    ccAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGA‑‑‑TCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGccc                                                             >  1:740529/1‑144 (MQ=255)
                                                                                    ccAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGCGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGA‑‑‑TCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGccc                                                             <  2:740529/144‑1 (MQ=255)
                                                                                           tgctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGA‑‑‑TCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAg                                                <  1:108885/150‑1 (MQ=255)
                                                                                                                       ccGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTCGA‑‑‑TCTCTACAGGAGCGGCCtt                                                                                         <  1:398838/81‑1 (MQ=255)
                                                                                                                       ccGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTCGA‑‑‑TCTCTACAGGAGCGGCCtt                                                                                         >  2:398838/1‑81 (MQ=255)
                                                                                                                              ggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGA‑‑‑TCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGaca             >  2:169307/1‑150 (MQ=255)
                                                                                                                               gctgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGA‑‑‑TCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGacac            <  1:193812/150‑1 (MQ=255)
                                                                                                                                 tgctTCGCAGCCCTATCACGACACAAGGTCGCTCCTGCAGAGGATTTCTTGA‑‑‑TCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAACTTACTTCTGCTTCACACCTTCGACACTg          >  1:173852/1‑150 (MQ=255)
                                                                                                                                 cgctTTGCGGCCCTATCACGACACAAGGCCGCTCCTGTAGAGA‑‑‑TCAAGAAATTCTCTGCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACt           >  1:707824/2‑149 (MQ=11)
                                                                                                                                 cgctTTGCGGCCCTATCACGACACAAGGCCGCTCCTGTAGAGA‑‑‑TCAAGAAATTCTCTGCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACt           <  2:707824/148‑1 (MQ=11)
                                                                                                                                     tCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGA‑‑‑TCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGa               <  1:624777/141‑1 (MQ=255)
                                                                                                                                     tCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGA‑‑‑TCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGa               >  2:624777/1‑141 (MQ=255)
                                                                                                                                         aGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGA‑‑‑TCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccag  >  2:118977/1‑150 (MQ=255)
                                                                                                                                         |                                                                                                                                                        
ATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGA‑‑‑TCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAG  >  NC_002947/5674579‑5674865

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: