Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 495,972 G→C 17.4% A22G (GCC→GGC)  PP_0408 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947495,9720GC17.4% 65.8 / 6.4 28A22G (GCC→GGC) PP_0408hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (10/13);  new base C (3/2);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCTGCAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTT  >  NC_002947/495824‑496113
                                                                                                                                                    |                                                                                                                                             
ggCCTGCAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGgc                                                                                                                                              <  2:113070/150‑1 (MQ=255)
    tGCAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCg                                                                                                                                          >  1:528284/1‑150 (MQ=255)
          gTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGa                                                                                                                                    >  1:722654/1‑150 (MQ=255)
           tCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGaa                                                                                                                                   <  1:629120/150‑1 (MQ=255)
            ctgctgCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAg                                                                                                                                  <  2:528284/150‑1 (MQ=255)
             tgctgcTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAgg                                                                                                                                 <  2:415092/150‑1 (MQ=255)
                  ctCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGa                                                                                                                            >  1:370190/1‑150 (MQ=255)
                  ctCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGa                                                                                                                            <  2:189567/150‑1 (MQ=255)
                         ttGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGCCCTGGCCGCCGat                                                                                                                                        <  1:199241/131‑2 (MQ=255)
                         ttGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGCCCTGGCCGCCGat                                                                                                                                        >  2:199241/1‑130 (MQ=255)
                                      gTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGCCCTGGCCGCCGatgc                                                                                                                                      >  1:297350/1‑117 (MQ=255)
                                      gTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGCCCTGGCCGCCGatgc                                                                                                                                      <  2:297350/120‑4 (MQ=255)
                                                  gTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGACGTGCGAt                                                                                            >  2:49102/1‑150 (MQ=255)
                                                     cATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTATCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATaaa                                                                                         <  2:128488/150‑1 (MQ=255)
                                                          tcctGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGtgtg                                                                                    <  2:722654/150‑1 (MQ=255)
                                                            ctGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAgt                                                                                       >  1:1849/1‑145 (MQ=255)
                                                            ctGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAgt                                                                                       <  2:1849/145‑1 (MQ=255)
                                                            ctGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGa                                                                                  >  1:821871/1‑150 (MQ=255)
                                                            ctGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGa                                                                                  <  2:504717/150‑1 (MQ=255)
                                                               gCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACAt                                                                               <  2:118218/150‑1 (MQ=255)
                                                                 tGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGa                                                                             <  2:370190/150‑1 (MQ=255)
                                                                     ttCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCa                                                                         >  2:231788/1‑150 (MQ=255)
                                                                        cAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAAcc                                                                      >  1:725244/1‑150 (MQ=255)
                                                                          gCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCa                                                                    >  2:381965/1‑150 (MQ=255)
                                                                                                     ccGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAAt                                         <  1:49102/150‑1 (MQ=255)
                                                                                                             aTGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGCGTGGTAGGACCGAGGCCAGGCAGCACATTGACGGCGTCGTGCGAGAAAGGTTGGACAGGATTCTACGCAAGGACAGAATGAAACCGGCAGGCTAATAGAAtgtg                                 >  2:97160/1‑150 (MQ=255)
                                                                                                                                      actCGGCGCCCTGGCCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGcc        >  1:829561/3‑150 (MQ=255)
                                                                                                                                           gcggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGTcac   <  2:139896/150‑3 (MQ=255)
                                                                                                                                            cggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCttt  <  1:682618/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GGCCTGCAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTT  >  NC_002947/495824‑496113

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: