Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,606,449 A→T 17.4% intergenic (‑157/+54) PP_1407 ← / ← PP_t64 sulfate transporter/tRNA‑Gly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,606,4490AT17.4% 50.5 / 6.3 23intergenic (‑157/+54)PP_1407/PP_t64sulfate transporter/tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base A (10/9);  new base T (2/2);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

TATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTC  >  NC_002947/1606332‑1606586
                                                                                                                     |                                                                                                                                         
tatCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCTTGGCTCCtt                                                                                                                                  >  2:294296/1‑127 (MQ=255)
tatCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCTTGGCTCCtt                                                                                                                                  <  1:294296/127‑1 (MQ=255)
              cGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCCTGGGCCAGGCCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGa                                                                                                                              >  1:607493/1‑117 (MQ=255)
              cGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGa                                                                                                                              <  2:607493/117‑1 (MQ=255)
                          gcacgcaAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAAc                                                                                                         <  2:402585/126‑1 (MQ=255)
                          gcacgcaAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAAc                                                                                                         >  1:402585/1‑126 (MQ=255)
                          gcacgcaAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAg                                                                                 <  1:570914/150‑1 (MQ=255)
                              gcaAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCggg                                                                             >  2:309963/1‑150 (MQ=255)
                                     cccGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGAtgtg                                                                                    >  1:177966/1‑136 (MQ=255)
                                     cccGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGAtgtg                                                                                    <  2:177966/136‑1 (MQ=255)
                                         aTGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATc                                                                  <  2:212518/150‑1 (MQ=255)
                                               aTTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAAccc                                                            >  1:413740/1‑150 (MQ=255)
                                               aTTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAAccc                                                            >  1:213128/1‑150 (MQ=255)
                                                        gATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTaaa                                                   >  1:3744/1‑150 (MQ=255)
                                                                    aTTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCtt                                               >  2:224333/1‑142 (MQ=255)
                                                                    aTTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCtt                                               <  1:224333/142‑1 (MQ=255)
                                                                            tGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGtttt                               <  2:413740/150‑1 (MQ=255)
                                                                                    gTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTaa                       >  2:279953/1‑150 (MQ=255)
                                                                                             ggATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATAccc              <  1:309963/150‑1 (MQ=255)
                                                                                              gATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCg             >  1:27904/1‑150 (MQ=255)
                                                                                                         aaaaaaGGAGCCAAGGCTCCTTTTGACTGCGATGACGATGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTc  <  1:523189/150‑1 (MQ=255)
                                                                                                            aaaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCggg                                                                             >  1:541417/1‑72 (MQ=255)
                                                                                                            aaaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCggg                                                                             <  2:541417/72‑1 (MQ=39)
                                                                                                                     |                                                                                                                                         
TATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTC  >  NC_002947/1606332‑1606586

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