Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,203,150 G→A 9.3% intergenic (+373/‑100) PP_2809 → / → PP_2810 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,203,1500GA 9.3% 115.9 / 2.3 43intergenic (+373/‑100)PP_2809/PP_2810membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (18/21);  new base A (2/2);  total (20/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.51e-01

GGTAACACATGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGC  >  NC_002947/3203004‑3203290
                                                                                                                                                  |                                                                                                                                            
ggTAACACATGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATATATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGt                                                                                                                                             >  6:135925/1‑148 (MQ=255)
ggTAACACATGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGt                                                                                                                                             <  8:155574/148‑1 (MQ=255)
    acacATGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGt                                                                                                                                         <  2:184802/148‑1 (MQ=255)
    acacATGTTTCCTTACTTCCCAATAATGAGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGt                                                                                                                                         <  4:100137/148‑1 (MQ=255)
     cacaTGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTg                                                                                                                                       <  7:112224/149‑1 (MQ=255)
        aTGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGAt                                                                                                                                     >  1:335013/1‑148 (MQ=255)
             tCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCC‑‑TGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGATCAACTGACCATGTcat                                                                                                                             >  3:350677/1‑147 (MQ=255)
              ccTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCC‑‑TGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGATCAACTGACCATGTcat                                                                                                                             <  4:350677/148‑3 (MQ=255)
                       ccAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAg                                                                                                                     >  5:13169/1‑149 (MQ=255)
                       ccAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAg                                                                                                                     >  4:246474/1‑149 (MQ=255)
                           taatGTGCAGGCGCGATGCACTCATTAGTCGATAAATAAATTCAACTCTCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGTCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATg                                                                                                                 <  6:141153/149‑1 (MQ=255)
                           taatGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATg                                                                                                                 >  7:151764/1‑149 (MQ=255)
                                   aGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCa                                                                                                                      <  3:263868/136‑1 (MQ=255)
                                   aGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCa                                                                                                                      >  4:263868/1‑136 (MQ=255)
                                        ccAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGa                                                                                                    <  5:135925/149‑1 (MQ=255)
                                                 cATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCAc                                                                                           >  8:148273/1‑149 (MQ=255)
                                                 cATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATAATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCAc                                                                                           <  1:22721/149‑1 (MQ=255)
                                                         gATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCAcccttcgg                                                                                    <  6:13169/148‑6 (MQ=255)
                                                           taaataaaTTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTc                                                                                 <  7:148273/149‑1 (MQ=255)
                                                              ataaatTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGATCAACTGACCATGTcatg                                                                                                                            >  5:163075/1‑100 (MQ=255)
                                                              ataaatTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGATCAACTGACCATGTcatg                                                                                                                            <  6:163075/103‑4 (MQ=255)
                                                                  atTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATg                                                                                                                 >  1:17257/1‑110 (MQ=255)
                                                                  atTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATg                                                                                                                 <  2:17257/110‑1 (MQ=255)
                                                                  atTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGagcagc                                                                          <  3:246474/149‑1 (MQ=255)
                                                                         tCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAgg                                                                                    <  5:213437/132‑1 (MQ=255)
                                                                         tCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAgg                                                                                    >  6:213437/1‑132 (MQ=255)
                                                                         tCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTcc                                                                    >  3:154770/1‑148 (MQ=255)
                                                                           gCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCaa                                                                  <  2:335013/148‑1 (MQ=255)
                                                                           gCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCaa                                                                  <  8:151764/148‑1 (MQ=255)
                                                                             tGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGcc                                                               >  7:160291/1‑149 (MQ=255)
                                                                                    tctAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGccc                                                              <  5:177891/143‑1 (MQ=255)
                                                                                    tctAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGccc                                                              >  6:177891/1‑143 (MQ=255)
                                                                                             aTCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGa                                                                               >  5:31064/1‑117 (MQ=255)
                                                                                             aTCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGa                                                                               <  6:31064/117‑1 (MQ=255)
                                                                                                   cTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGag                                                                              >  5:154852/1‑112 (MQ=255)
                                                                                                   cTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGag                                                                              <  6:154852/112‑1 (MQ=255)
                                                                                                   cTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAAt                                         <  8:160291/149‑1 (MQ=255)
                                                                                                        aCGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGaaa                                     >  6:65586/1‑148 (MQ=255)
                                                                                                         cGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAAc                                    >  8:68528/1‑148 (MQ=255)
                                                                                                          gCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAAc                                    >  5:226141/1‑147 (MQ=255)
                                                                                                          gCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACaa                                  <  5:65586/149‑1 (MQ=255)
                                                                                                                                tgtTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAgtg            <  6:226141/149‑1 (MQ=255)
                                                                                                                                          atgacatgGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATccgc  >  1:228662/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
GGTAACACATGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGC  >  NC_002947/3203004‑3203290

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