Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,079,177 T→A 21.0% intergenic (‑374/+23) mgtE ← / ← PP_t56 magnesium transporter/tRNA‑Arg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,079,1770TA21.0% 37.7 / 6.7 19intergenic (‑374/+23)mgtE/PP_t56magnesium transporter/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base T (5/10);  new base A (2/2);  total (7/12)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

CTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACT  >  NC_002947/5079028‑5079314
                                                                                                                                                     |                                                                                                                                          
cTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTAtt                                                                                                                                            <  1:520269/150‑1 (MQ=255)
      gTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCgg                                                                                                                                      <  1:216422/150‑1 (MQ=255)
                 gttgtGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATg                                                                                                                           <  1:575997/150‑1 (MQ=255)
                       gATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCTTTTTg                                                                                                                              <  1:167058/141‑1 (MQ=255)
                       gATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCTTTTTg                                                                                                                              >  2:167058/1‑141 (MQ=255)
                             aaaGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGa                                                                                                               <  1:465878/150‑1 (MQ=255)
                             aaaGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGa                                                                                                               <  2:597893/150‑1 (MQ=255)
                                  acaTTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCTGACTcag                                                                                                          <  2:392767/150‑1 (MQ=255)
                                                                        cacCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGAGCACCACCAAAAAGCCCCCTCTATTGAGCTCGCTATTGTTTAGGGAAGGGCTGAGTCTGCAGGAGTTCGAACCTGCGACCGCTCGGTTCGTAGCCGa                                                                    >  2:531545/1‑150 (MQ=255)
                                                                           cGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTAc                                                                >  1:268800/1‑150 (MQ=255)
                                                                                           gagCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCt                                                >  2:205224/1‑150 (MQ=255)
                                                                                                  gATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCttttt                                                                                                                               <  1:254134/65‑1 (MQ=25)
                                                                                                  gATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCttttt                                                                                                                               >  2:254134/1‑65 (MQ=25)
                                                                                                       aCAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTAATGGGGTATGGCGGACTAAGTAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGa                                    <  1:200471/150‑1 (MQ=255)
                                                                                                                       aTTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCa                    <  2:268800/150‑1 (MQ=255)
                                                                                                                               acaaacaaAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAg                   >  1:30437/1‑143 (MQ=255)
                                                                                                                               acaaacaaAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAg                   <  2:30437/143‑1 (MQ=255)
                                                                                                                                          gCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCGGCAGGATTTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACt  <  1:531545/150‑1 (MQ=255)
                                                                                                                                               ctctATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTc                                                                                                 <  2:106582/49‑1 (MQ=255)
                                                                                                                                               ctctATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTc                                                                                                 >  1:106582/1‑49 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
CTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACT  >  NC_002947/5079028‑5079314

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: