Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,005,787 G→T 10.2% V713V (GTG→GTT)  PP_4410 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,005,7870GT10.2% 113.3 / 3.2 39V713V (GTG→GTT) PP_4410hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (17/18);  new base T (2/2);  total (19/20)
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.49e-01

CTGAACTTCCAAATCTGGTTTAAACAGAATACCCAGGAGGCAAATCCCATTGTAGGCTTCCTTCCTCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAACTCCCCACGCTGCC‑‑CGAGCCTCCTTCGTTACAGACATGTCGCTTAATCTCCCACCAGAAGCGGTATCATTGC  >  NC_002947/5005644‑5005929
                                                                                                                                               |                                                                                                                                                
cTGAACTTCCAAATCTGGTTTAAACAGAATACCCAGGAGGCAAATCCCATTGTAGGCTTCCTTCCTCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGAtt                                                                                                                                               >  3:297841/1‑147 (MQ=255)
cTGAACTTCCAAATCTGGTTTAAACAGAATACCCAGGAGGCAAATCCCATTGTAGGCTTCCTTCCTCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTaa                                                                                                                                             <  1:129265/149‑1 (MQ=255)
 tGAACTTCCAAATCTGGTTTAAACAGAATACCCAGGAGGCAAATCCCATTGTAGGCTTCCTTCCTCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAg                                                                                                                                            >  6:60971/1‑149 (MQ=255)
  gAACTTCCAAATCTGGTTTAAACAGAATACCCAGGAGGCAAATCCCATTGTAGGCTTCCTTCCTCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTTATTAAgg                                                                                                                                           >  6:54692/1‑149 (MQ=255)
  gAACTTCCAAATCTGGTTTAAACAGAATACCCAGGAGGCAAATCCCATTGTAGGCTTCCTTCCTCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTTATTAAgg                                                                                                                                           >  6:60001/1‑149 (MQ=255)
   aaCTTCCAAATCTGGTTTAAACAGAATACCCAGGAGGCAAATCCCATTGTAGGCTTCCTTCCTCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAg                                                                                                                                            <  5:60971/147‑1 (MQ=255)
                    tAAACAGAATACCCAGGAGGCAAATCCCATTGTAGGCTTCCTTCCTCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACaataat                                                                                                                          <  8:24770/148‑1 (MQ=255)
                             tACCCAGGAGGCAAATCCCATTGTAGGCTTCCTTCCTCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAgg                                                                                                                <  8:103832/149‑1 (MQ=255)
                               cccAGGAGGCAAATCCCATTGTAGGCTTCCTTCCTCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGga                                                                                                              <  5:115652/149‑1 (MQ=255)
                                ccAGGAGGCAAATCCCATTGTAGGCTTCCTTCCTCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAgg                                                                                                                >  1:316076/1‑146 (MQ=255)
                                                 ttGTAGGCTTCCTTCCTCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATctca                                                                                            <  2:292956/149‑1 (MQ=255)
                                                        cttccttcCTCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTTATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTNAggg                                                                                     <  5:60001/149‑1 (MQ=255)
                                                        cttccttcCTCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTTATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAggg                                                                                     <  5:54692/149‑1 (MQ=255)
                                                         ttccttccTCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCAATCCTGACAAGCTCAAGGTGATTAAGGGCTCCACTAACAAGAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGc                                                                                    >  2:98338/1‑149 (MQ=255)
                                                            cttcctCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTa                                                                                 <  4:172912/149‑1 (MQ=255)
                                                               cctCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTAt                                                                                <  2:321407/147‑1 (MQ=255)
                                                               cctCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTAt                                                                                <  1:357622/147‑1 (MQ=255)
                                                                 tctTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCaa                                                                            >  8:340474/1‑149 (MQ=255)
                                                                   ttAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCaa                                                                            <  2:316076/147‑1 (MQ=255)
                                                                             cccTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTaa                                                                                                                                 <  3:229623/84‑1 (MQ=255)
                                                                             cccTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTaa                                                                                                                                 >  4:229623/1‑84 (MQ=255)
                                                                              ccTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAACTCCCCACGCTGc                                                               >  7:111037/1‑149 (MQ=255)
                                                                                  cGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAAc                                                                           >  5:209610/1‑133 (MQ=255)
                                                                                  cGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAAc                                                                           <  6:209610/133‑1 (MQ=255)
                                                                                         cTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAggg                                                                                                                                          <  7:84007/63‑1 (MQ=255)
                                                                                         cTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAggg                                                                                                                                          >  8:84007/1‑63 (MQ=255)
                                                                                                           caggATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAACTCCCCACGCTGCC‑‑CGAGCCTCCTTCGTTACAGACATGTCGc                                <  8:158143/147‑1 (MQ=255)
                                                                                                              gATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAACTCCCCACGCTGCC‑‑CGAGCCTCCTTCGTTACAGACATGTCGCtt                              >  6:200095/1‑148 (MQ=255)
                                                                                                              gATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCAAATCTCACTCAGGGCTTATGCAACTCGCCACGCTGCC‑‑GGAGCCGCCTTCGTTACAGACATGTCGCTTa                             >  8:283817/1‑149 (MQ=255)
                                                                                                                  aacaaAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAACTCACCACGCTGCC‑‑CGAGCCTCCTTCGTTACAGACATGTCGCTTAAtct                         >  2:144623/1‑149 (MQ=255)
                                                                                                                    caaAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAACTCCCCACGCTGCC‑‑CGAGCCTCCTTCGTTACAGACATGTCGCTTAATCTcc                       >  4:20145/1‑149 (MQ=255)
                                                                                                                       aaTCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAACTCCCCACGCTGCC‑‑CGAGCCTCCTTCGTTACAGACATGTCGc                                <  2:8994/137‑1 (MQ=255)
                                                                                                                       aaTCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAACTCCCCACGCTGCC‑‑CGAGCCTCCTTCGTTACAGACATGTCGc                                >  1:8994/1‑137 (MQ=255)
                                                                                                                         tccatccTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAACTCCCCACGCTGCC‑‑CGAGCCTCCTTCGTTACAGACATGTCGCTTAATCTCccacca                  <  7:340474/149‑1 (MQ=255)
                                                                                                                            atccTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAACTCCCCACGCTGCC‑‑CGAGCCTCCTTCGTTACAGACATGTCGCTTAATCTCCCACCAg                 >  6:862/1‑147 (MQ=255)
                                                                                                                             tccTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAACCCCCCACGCTGCC‑‑CGAGCCTCCTTCGTTACAGACATGTCGCTTAATCTCCCACCAGAAg              >  3:161904/1‑149 (MQ=255)
                                                                                                                              ccTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAACTCCCCACGCTGCC‑‑CGAGCCTCCTTCGTTACAGACATGTCGCTTAATCTCCCACCAGAAGc             >  6:3437/1‑149 (MQ=255)
                                                                                                                                     aggtcaaggtGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAACTCCCCACGCTGCC‑‑CGAGCCTCCTTCGTTACAGACATGTCGCTTAATCTCCCACCAGAAGCGGTATCa      <  4:297841/149‑1 (MQ=255)
                                                                                                                                           aggtGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAACTCCCCACGCTGCCCTCGAGCCTCCTTCGTTACAGACATGTCGCTTAATCTCCCACCAGAAGCGGTATCATTGc  <  1:98338/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CTGAACTTCCAAATCTGGTTTAAACAGAATACCCAGGAGGCAAATCCCATTGTAGGCTTCCTTCCTCTTAAAGAAGGCCCTACGAAATCCTTTTTCACCTGGGAAGATTGGATCAACAAAATCCATCCTGACAAGGTCAAGGTGATTAAGGGCTCCACTAACAATAATATATATCAGGGAGACTACTGCCCAATCTCACTCAGGGCTTATGCAACTCCCCACGCTGCC‑‑CGAGCCTCCTTCGTTACAGACATGTCGCTTAATCTCCCACCAGAAGCGGTATCATTGC  >  NC_002947/5005644‑5005929

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑