Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,904,699 C→G 15.2% L60F (TTG→TTC)  PP_4312 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,904,6990CG15.2% 77.5 / 9.7 33L60F (TTG→TTC) PP_4312hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (16/12);  new base G (3/2);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-01

GCTGTAGCTTTCGTCCCGGCTGCCTTCCTCGACATGCCGCGTGGGGTTGGCGGCAAACACGTAGTCGGCCTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCTTTCGACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGAC  >  NC_002947/4904567‑4904847
                                                                                                                                    |                                                                                                                                                    
gctgtagctTTCGTCCCGGCTGCCTTCCTCGACATGCCGCGTGGGGTTGGCGGCAAACACGTAGTCGGCCTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTc                                                                                                                                     <  1:224326/150‑1 (MQ=255)
       ctTTCGTCCCGGCTGCCTTCCTCGACATGCCGCGTGGGGTTGGCGGCAAACACGTAGTAGGCCTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGc                                                                                                                              <  1:319601/150‑1 (MQ=255)
        tttCGTCCCGGCTGCCTTCCTCGACATGCCGCGTGGGGTTGGCGGCAAACACGTAGTCGGCCTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCt                                                                                                                             >  1:437014/1‑150 (MQ=255)
        tttCGTCCCGGCTGCCTTCCTCGACATGCCGCGTGGGGTTGGCGGCAAACACGTAGTCGGCCTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCt                                                                                                                             <  1:459731/150‑1 (MQ=255)
             tCCCGGCTGCCTTCCTCGACATGCCGCGTGGGGTTGGCGGCAAACACGTAGTCGGCCTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCg                                                                                                                        >  2:279664/1‑150 (MQ=255)
                       cTTCCTCGACATGCCGCGTGGGGTTGGCGGCAAACACGTAGTCGGCCTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTGAACTGccaag                                                                                                                                            >  2:357094/1‑117 (MQ=255)
                       cTTCCTCGACATGCCGCGTGGGGTTGGCGGCAAACACGTAGTCGGCCTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTGAACTGccaag                                                                                                                                            >  2:357021/1‑117 (MQ=255)
                       cTTCCTCGACATGCCGCGTGGGGTTGGCGGCAAACACGTAGTCGGCCTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTGAACTGccaag                                                                                                                                            <  1:357021/120‑4 (MQ=255)
                       cTTCCTCGACATGCCGCGTGGGGTTGGCGGCAAACACGTAGTCGGCCTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTGAACTGccaag                                                                                                                                            <  1:357094/120‑4 (MQ=255)
                       cTTCCTCGACATGCCGCGTGGGGTTGGCGGCAAACACGTAGTCGGCCTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCTCCCGGCTTCTCGACCTGACGCa                                                                                                              >  2:97561/1‑150 (MQ=255)
                        ttCCTCGACATGCCGCGTGGGGTTGGCGGCAAACACGTAGTCGGCCTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAg                                                                                                             >  1:102278/1‑150 (MQ=255)
                                    ccGCGTGGGGTTGGCGGCAAACACGTAGTCGGCCTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGt                                                                                                 >  2:277443/1‑150 (MQ=255)
                                                                    ccTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGc                                                                                                                              <  2:348160/89‑1 (MQ=255)
                                                                    ccTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGc                                                                                                                              >  1:348160/1‑89 (MQ=255)
                                                                    ccTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCttt                                                                 >  2:160964/1‑150 (MQ=255)
                                                                           aaaGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCTTTCGACAGc                                                          >  2:56159/1‑150 (MQ=255)
                                                                                  tCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCttt                                                                 <  2:338693/136‑1 (MQ=255)
                                                                                  tCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCttt                                                                 >  1:338693/1‑136 (MQ=255)
                                                                                            cacTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCaa                                                                                                                                                    <  1:110416/43‑1 (MQ=255)
                                                                                            cacTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCaa                                                                                                                                                    >  2:110416/1‑43 (MQ=255)
                                                                                              cTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCTTTCGACAGCCTCCAGCGCCAAcggcgcg                                       >  1:153041/1‑150 (MQ=255)
                                                                                                        ccAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCTTTCGACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAg                             <  1:277443/150‑1 (MQ=255)
                                                                                                         cAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCTTTCGACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGt                            <  1:279664/150‑1 (MQ=255)
                                                                                                          aGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGTGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCTTTCGACAGCCTCCAGCGCCAACGGCGCGGCGATg                                 <  2:169446/144‑1 (MQ=255)
                                                                                                          aGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCTTTCGACAGCCTCCAGCGCCAACGGCGCGGCGATg                                 >  1:169446/1‑144 (MQ=255)
                                                                                                          aGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCTTTCGACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGTg                           <  2:153041/150‑1 (MQ=255)
                                                                                                                      agaCCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCTTTCGACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTc               >  1:84091/1‑150 (MQ=255)
                                                                                                                        aCCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCTTTCGACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAg                   >  1:445501/1‑144 (MQ=255)
                                                                                                                        aCCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCTTTCGACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAg                   <  2:445501/144‑1 (MQ=255)
                                                                                                                         ctcgggCAGTTGAACTGCCCGAGCCCTTCGCCCGGCTTCTCGGCATGACGCAGAAAAACGAGGGTCTGGGTGCCATCCACCGGCTGTGCCGGGCTTTCGGCAGCCTCCAGCGCCATCAGCGCGGAGATGGCAGTGAGCAACAGGCTcagc            >  2:223728/5‑150 (MQ=255)
                                                                                                                          cTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCTTTCGACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTcagca           >  1:276994/1‑150 (MQ=255)
                                                                                                                                  ttCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCTTTCAACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTga   >  2:502233/1‑150 (MQ=255)
                                                                                                                                   tCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAATACCAGGGTCTGGGTGCCAGCCACCGGTTGTGCCCGGCTTTCGACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTgac  <  2:68084/150‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                                    
GCTGTAGCTTTCGTCCCGGCTGCCTTCCTCGACATGCCGCGTGGGGTTGGCGGCAAACACGTAGTCGGCCTTGCCAAAGCGCTCGGGCAGTACACTGGCCAGGTCCAGCGCACGGTTGAGACCTTGGCAGTTCAACTGCCCGAGCCCTTCGCCCGGCTTCTCGCCATGACGCAGAAACACCAGGGTCTGGGTGCCATCCACCGGTTGTGCCCGGCTTTCGACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGAC  >  NC_002947/4904567‑4904847

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: ‑