Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 274,788 A→T 30.2% L387Q (CTG→CAG)  PP_0222 ← DszA family monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947274,7880AT30.2% 32.1 / 13.5 23L387Q (CTG→CAG) PP_0222DszA family monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (8/8);  new base T (4/3);  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 = 9.93e-01

TCGGCCTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGC  >  NC_002947/274644‑274925
                                                                                                                                                |                                                                                                                                         
tCGGCCTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCg                                                                                                                                      >  1:401119/1‑150 (MQ=255)
tCGGCCTGGAACAGTTTTTAGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCg                                                                                                                                      >  2:539105/1‑150 (MQ=255)
    ccTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGGCCTGGACGGttt                                                                                                                                  >  2:394746/1‑150 (MQ=255)
        gAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGGCCTGGACGGTTTCaac                                                                                                                              <  1:394746/150‑3 (MQ=255)
                tttCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCAcg                                                                                                                      <  2:401119/150‑1 (MQ=255)
                         gggTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACt                                                                                                             >  1:276940/1‑150 (MQ=255)
                          ggTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTg                                                                                                            >  2:185322/1‑150 (MQ=255)
                            taccgtCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGATCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTAAGGTTGAAAAGATACAGGCCGGTTTCATTGATCTAAGAGATAAACTggt                                                                                                          >  2:267352/6‑148 (MQ=255)
                               cTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCgg                                                                                                       <  1:397900/150‑1 (MQ=255)
                                ttCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGc                                                                                                      >  2:323567/1‑150 (MQ=255)
                                                        gCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAg                                                                              >  2:136183/1‑150 (MQ=255)
                                                         ccACGCTGCTTAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGTGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAgg                                                                             <  1:539105/150‑1 (MQ=255)
                                                         ccACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAgg                                                                             <  1:185322/150‑1 (MQ=255)
                                                                   gAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGc                                                                   <  2:350993/150‑1 (MQ=255)
                                                                              aTGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAg                                                        <  1:323567/150‑1 (MQ=255)
                                                                                        cGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCATCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGtt                                              <  1:306803/150‑1 (MQ=255)
                                                                                                tCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAgg                                      >  2:309010/1‑150 (MQ=255)
                                                                                                 ctttgAATCTCTCCTGTGTGCCGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCCATCTACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGgc                                     <  1:267352/146‑1 (MQ=255)
                                                                                                                                    atGAAACCGGCCTTGACGGTTTCAGCGATCCACGCGATCGACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCACCTGCCGACCTGCTGCCACAGGCGCCGGCTGGTCCAGGCGTTGTCTTGGAACTTgcgc  >  2:438464/2‑150 (MQ=255)
                                                                                                                                    atGAAACCGGCCTGGACGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgc                                                                                         <  1:525781/62‑1 (MQ=37)
                                                                                                                                    atGAAACCGGCCTGGACGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGgc                                                                                         >  2:525781/2‑63 (MQ=25)
                                                                                                                                     tGAAACCGGCCTGGACGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGtt                  <  2:113944/133‑1 (MQ=255)
                                                                                                                                     tGAAACCGGCCTGGACGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGtt                  >  1:113944/1‑133 (MQ=255)
                                                                                                                                                |                                                                                                                                         
TCGGCCTGGAACAGTTTTTCGCGCAGGGTGCCTTCGGTGTAGGCAGTCTTGTAGCGGCCACGCTGCTGAAGTTGCGGGATGACCAGCTCGATGAAGTCCTCGAAGCTCTCCGGGGTGACGGTGCGGGTCAGGTTGAAACCGTCCAGGCCGGTTTCATCGATCCACGCGATCAACTGCTCGGCCACCTGTTCAGGCGCGCCGACCAGGGTTACGTAGCGGCCGCCCAAGGCGTGCTGTTGCAACAGGCGCCGGCGGGTCCAGGCGTTGTCTTGCAACTGGCGC  >  NC_002947/274644‑274925

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