Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,827,626 C→T 8.3% V1203V (GTG→GTA pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,827,6260CT 8.3% 133.3 / 2.6 48V1203V (GTG→GTApvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (23/21);  new base T (2/2);  total (25/23)
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.85e-01

CAGGTGCAACTGGTCGAGGGCAATGGCGCTGACACCGTCGACCTGCGGCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAGGTATCGGGGCCGACGCCTTTGTCACGCAGGTAGTGGGCCAGGCGGTTGGCTTGCTGGTGCAGCC  >  NC_002947/4827479‑4827771
                                                                                                                                                   |                                                                                                                                                 
cAGGTGCAACTGGTCGAGGGCAATGGCGCTGACACCGTCGACCTGCGGCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCAc                                                                                                                                                 >  2:970566/1‑150 (MQ=255)
   gTGCAACTGGTCGAGGGCAATGGCGCTGACACCGTCGACCTGCGGCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTa                                                                                                                                              >  1:471639/1‑150 (MQ=255)
     gCAACTGGTCGAGGGCAATGGCGCTGACACCGTCGACCTGCGGCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAGCAGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGATCGGCAGGGTAATCGACATCCAGCGGCACGTAgg                                                                                                                                            >  1:222736/1‑150 (MQ=255)
       aaCTGGTCGAGGGCAATGGCGCTGACACCGTCGACCTGCGGCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCa                                                                                                                                          >  1:338876/1‑150 (MQ=255)
                aGGGCAATGGCGCTGACACCGTCGACCTGCGGCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCGGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGcc                                                                                                                                 >  2:899884/1‑150 (MQ=255)
                         gcgcTGACACCGTCGACCTGCGGCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATg                                                                                                                        <  1:98957/150‑1 (MQ=255)
                          cgcTGACACCGTCGACCTGCGGCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATgg                                                                                                                       <  2:696748/150‑1 (MQ=255)
                           gttGACACCGTCGACCTGCGGCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGc                                                                                                                      >  2:650820/3‑150 (MQ=255)
                                cacCGTCGACCTGCGGCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGACTTGAGAATGGCTAGCa                                                                                                                 <  1:869816/150‑1 (MQ=255)
                                         ccTGCGGCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAg                                                                                                                                             <  1:186631/113‑1 (MQ=255)
                                         ccTGCGGCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAg                                                                                                                                             >  2:186631/1‑113 (MQ=255)
                                         ccTGCGGCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCa                                                                                                        <  2:437711/150‑1 (MQ=255)
                                          cTGCGGCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGTACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAg                                                                                                       >  2:670238/1‑150 (MQ=255)
                                               gCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGTACGTAGGCACCGCCGGCCTTga                                                                                                                             >  2:1226776/1‑123 (MQ=255)
                                               gCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGTACGTAGGCACCGCCGGCCTTga                                                                                                                             <  1:1226776/123‑1 (MQ=255)
                                               gCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAAtt                                                                                                  <  1:650820/150‑1 (MQ=255)
                                                      gCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGa                                                                                           <  1:859348/150‑1 (MQ=255)
                                                         cagcagGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAgcg                                                                                        >  1:501318/1‑150 (MQ=255)
                                                         cagcagGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCAAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAgcg                                                                                        >  2:130653/1‑150 (MQ=255)
                                                            cagGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTc                                                                                     >  1:168235/1‑150 (MQ=255)
                                                              ggCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGa                                                                                   >  1:20167/1‑150 (MQ=255)
                                                               gctgctGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGAt                                                                                  <  2:20167/150‑1 (MQ=255)
                                                               gctgctGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGAt                                                                                  <  2:15497/150‑1 (MQ=255)
                                                                    tgtgGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGAACGGACGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGGGCTCGATGGCaa                                                                             >  2:979914/1‑150 (MQ=255)
                                                                          tCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCca                                                                       >  2:817544/1‑150 (MQ=255)
                                                                            aGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCcaca                                                                     <  1:164626/150‑1 (MQ=255)
                                                                             gTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCcacac                                                                    <  2:168235/150‑1 (MQ=255)
                                                                                gCAAACGGGCGTCGCAGTAGGCCAGCAGGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAGTGGCTACCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAgg                                                                 <  2:222736/150‑1 (MQ=255)
                                                                                       ggCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAGGTATCggg                                                          <  2:1216323/150‑1 (MQ=255)
                                                                                       ggCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAGGTATCggg                                                          <  2:501318/150‑1 (MQ=255)
                                                                                               aGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAGGTATCGGGGCCGACGc                                                  >  1:1269945/1‑150 (MQ=255)
                                                                                                            gTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGAt                                                                                  >  1:857968/1‑105 (MQ=255)
                                                                                                            gTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGAt                                                                                  <  2:857968/105‑1 (MQ=255)
                                                                                                             tAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAg                                                                  <  2:764327/120‑1 (MQ=255)
                                                                                                             tAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAg                                                                  >  1:764327/1‑120 (MQ=255)
                                                                                                                      cGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGAAATGTCGTGTAACGCTCGATGGCAATGGCCACCCAGGTATCGGGGCCGACGCCGTTGTCACGCAGGTAGTGGGcc                           >  1:1186048/1‑150 (MQ=255)
                                                                                                                       gTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAGGTATCGGGGCCGACGCCTTTGTCACGCAGGTAGTGGGCCa                          >  2:957062/1‑150 (MQ=255)
                                                                                                                        ttCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATg                                                                                 <  1:240381/94‑1 (MQ=255)
                                                                                                                        ttCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATg                                                                                 >  2:240381/1‑94 (MQ=255)
                                                                                                                          cGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAGGTATCGGGGCCGACGCCTTTGTCACGCAGGTAGTGGGCCAGGc                       <  1:970566/150‑1 (MQ=255)
                                                                                                                             cggggTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAGGTATCGGGGCCGACGCCTTTGTCACGCAGGTAGTGGGCCAGGCGGt                    <  1:899884/148‑1 (MQ=255)
                                                                                                                                    aTCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAGGTATCGGGGCCGACGCCTTTGTCACGCAGGTAGTGGGCCAGGCGGtt                   >  1:4809/1‑144 (MQ=255)
                                                                                                                                    aTCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAGGTATCGGGGCCGACGCCTTTGTCACGCAGGTAGTGGGCCAGGCGGtt                   <  2:4809/144‑1 (MQ=255)
                                                                                                                                    aTCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAGGTATCGGGGCCGACGCCTTTGTCACGCAGGTAGTGGGCCAGGCGGTTGGCTTg             >  1:401934/1‑150 (MQ=255)
                                                                                                                                      cGACATCCAGCGGTACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAGGTAAAGGGGACGCCCCCTTTGTCCCGCAGGTACTGGGCCAGGCGGTTGGCTTGCt           <  1:670238/150‑1 (MQ=255)
                                                                                                                                             cAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAGGTATCGGGGCCGACGCCTTTGAGTCACAGGTAGTGGGCCAGGCGGTTGGCTTGCTGGTGCAg    >  1:708079/1‑150 (MQ=255)
                                                                                                                                               gCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAGGTATCGGGGCCGACGCCTTTGTCACGCAGGTAGTGGGCCAGGCGGTTGGCTTGCTGGTGCAGcc  >  1:562983/1‑150 (MQ=255)
                                                                                                                                               gCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAGGTATCGGGGCCGACGCCTTTGTCACGCAGGTAGTGGGCCAGGCGGTTGGCTTGCTGGTGCAGcc  <  1:817544/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CAGGTGCAACTGGTCGAGGGCAATGGCGCTGACACCGTCGACCTGCGGCAGCTTGCCCAGCAGGCTGCTGTGGCTCAGTAGCAAACGGGCGTCGCAGTCGGCCAGCATGTAGGCCAGGCGTTCGGCAGGGTAATCGACATCCAGCGGCACGTAGGCACCGCCGGCCTTGAGAATGGCTAGCAGGCCGACCAGCAATTGCGGTGAGCGCTCGATGGCAATGGCCACACAGGTATCGGGGCCGACGCCTTTGTCACGCAGGTAGTGGGCCAGGCGGTTGGCTTGCTGGTGCAGCC  >  NC_002947/4827479‑4827771

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: