Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,041,669 T→C 11.4% D104G (GAC→GGC)  rph ← ribonuclease PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,041,6690TC11.4% 89.7 / 3.3 35D104G (GAC→GGC) rphribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (21/10);  new base C (2/2);  total (23/12)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

GCCGCCGCGCTTCTTGATCACTGCCAAAGCGTCGCACAAGGCAACCATGGCGCCGGTGATGGACGCGGTACGGGTACCGCCATCGGCCTGGATCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTGGTTACGTTCGCCGGTGGAGCGCGGCAGCATGCCATACTCGGCGGTCAGCCAGCCTTGGCCCTG  >  NC_002947/6041530‑6041817
                                                                                                                                           |                                                                                                                                                    
gccgccGCGCTTCTTGATCACTGCCAAAGCGTCGCACAAGGCAACCATGGCGCCGGTGATGGACGCGGTACGGGTACCGCCATCGGCCTGGATCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAgcgcggc                                                                                                                                            <  1:286041/150‑1 (MQ=255)
     cgcgCTTCTTGATCACTGCCAAAGCGTCGCACAAGGCAACCATGGCGCCGGTGATGGACGCGGTACGGGTACCGCCATCGGCCTGGATCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCa                                                                                                                                       >  2:470728/1‑150 (MQ=255)
                 tCACTGCCAAAGCGTCGCACAAGGCAACCATGGCGCCGGTGATGGACGCGGTACGGGTACCGCCATCGGCCTGGATCACGTCGCAATCGACGTACAGGGTGATGTCACTGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGa                                                                                                                           >  1:697135/1‑150 (MQ=255)
                     tGCCAAAGCGTCGCACAAGGCAACCATGGCGCCGGTGATGGACGCGGTACGGGTACCGCCATCGGCCTGGATCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAg                                                                                                                       >  1:588002/1‑150 (MQ=255)
                                               tGGCGCCGGTGATGGACGCGGTACGGGTACCGCCATCGGCCTGGATCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCAc                                                                                             <  2:615097/150‑1 (MQ=255)
                                                      ggTGATGGACCCGGTACGGGTACCGCCATCGGCCTGGATCACGTCGCAATCGACGTACTGATAAACATTACTCAGCTCGATCATGTCCAGCGAGGCGCGCAGCGACCGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCtt                                                                                      <  1:242751/150‑1 (MQ=255)
                                                       gTGATGGACGCGGTACGGGTACCGCCATCGGCCTGGATCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTg                                                                                     >  2:246956/1‑150 (MQ=255)
                                                          aTGGACGCGGTACGGGTACCGCCATCGGCCTGGATCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTACTTGCCa                                                                                  >  2:288369/1‑150 (MQ=255)
                                                              aCGCGGTACGGGTACCGCCGGTGGCCTGGATCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGAAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACggc                                                                              <  1:228398/150‑1 (MQ=255)
                                                                       gggTACCGCCATCGGCCTGGATCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTcgc                                                                     >  1:319051/1‑150 (MQ=255)
                                                                         gTACCGCCATCGGCCTGGATCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCgcgc                                                                   >  1:679543/1‑150 (MQ=255)
                                                                            ccGCCATCGGCCTGGATCACGTCGCAATCGACGTACAGGGTGATGTCACTGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTgg                                                                <  2:697135/150‑1 (MQ=255)
                                                                              gCCATCGGCCTGGATCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTGGtt                                                              >  1:322892/1‑150 (MQ=255)
                                                                                cATCGGCCTGGATCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTGGTTAc                                                            >  2:250223/1‑150 (MQ=255)
                                                                                       cTGGATCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTGGTTACGTTCGcc                                                     >  1:54357/1‑150 (MQ=255)
                                                                                          gATCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCt                                                                  <  2:634320/134‑1 (MQ=255)
                                                                                          gATCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCt                                                                  >  1:634320/1‑134 (MQ=255)
                                                                                           aTCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTGGTTACGTTCGCCGGTg                                                 >  2:691235/1‑150 (MQ=255)
                                                                                            tCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCa                                                                                                         >  2:148996/1‑93 (MQ=255)
                                                                                            tCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCa                                                                                                         <  1:148996/93‑1 (MQ=255)
                                                                                                        cGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTGGTTACGTTCGCCGGTGGAGCGCGgcagca                                    >  2:725419/1‑150 (MQ=255)
                                                                                                            gTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTGGtt                                                              <  1:515769/120‑1 (MQ=255)
                                                                                                            gTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTGGtt                                                              >  2:515769/1‑120 (MQ=255)
                                                                                                             tACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTGGTTACGTTCGCCGGTGGAGCGCGGCAGCATGCCa                               >  1:548198/1‑150 (MQ=255)
                                                                                                                  ggTGATGTCACCGAGCTTGCTCATGCCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCtt                                                                                      >  2:743257/1‑90 (MQ=255)
                                                                                                                  ggTGATGTCACCGAGCTTGCTCATGCCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCtt                                                                                      <  1:743257/90‑1 (MQ=255)
                                                                                                                       tGTCACCGAGCTTGCTCATGCCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAg                                                                                                        >  2:690736/1‑67 (MQ=255)
                                                                                                                       tGTCACCGAGCTTGCTCATGCCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAg                                                                                                        <  1:690736/67‑1 (MQ=255)
                                                                                                                           aCCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTGGTTACGTTCGCCGGTGGAGCGCGGCAGCATGCCATACTcg                         <  2:451564/142‑1 (MQ=255)
                                                                                                                           aCCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTGGTTACGTTCGCCGGTGGAGCGCGGCAGCATGCCATACTcg                         >  1:451564/1‑142 (MQ=255)
                                                                                                                           aCCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTGGTTACGTTCGCCGGTGGAGCGCGGCAGCATGCCATACTCGGCGGTcag                 >  2:676637/1‑150 (MQ=255)
                                                                                                                                  ttGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGAGCTGGATATCCAGCGGGCGGCCACCCTGCTTGCTACGGCTGGCTTCGCGCTGGTTACGTTCGCCGGTGGAGCGCGGCAGCATGCCATACTCGGCGGTCAGCCAGCCt          >  1:655938/1‑150 (MQ=255)
                                                                                                                                      tCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTGGTTACGTTCGCCGGTGGAGCGCGGCAGCATGCCATACTCGGCGGTCAGCCAGCCTTGGc      >  1:431413/1‑150 (MQ=255)
                                                                                                                                          gTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTGGTTACGTTCGCCGGTGGAGCGCGGCAGCATGCCATACTCGGCGGTCAGCCAGCCTTGGCCCTg  <  2:54357/150‑1 (MQ=255)
                                                                                                                                          gTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTGGTTACGTTCGCCGGTGGAGCGCGGCAGCATGCCATACTCGGCGGTCAGCCAGCCTTGGCCCTg  >  1:466080/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
GCCGCCGCGCTTCTTGATCACTGCCAAAGCGTCGCACAAGGCAACCATGGCGCCGGTGATGGACGCGGTACGGGTACCGCCATCGGCCTGGATCACGTCGCAATCGACGTACAGGGTGATGTCACCGAGCTTGCTCATGTCCAGCGCGGCGCGCAGCGAACGGCCGATCAGGCGCTGGATTTCCAGGGTGCGGCCACCCTGCTTGCCACGGCTGGCTTCGCGCTGGTTACGTTCGCCGGTGGAGCGCGGCAGCATGCCATACTCGGCGGTCAGCCAGCCTTGGCCCTG  >  NC_002947/6041530‑6041817

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: