Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,203,159 T→C 9.5% intergenic (+382/‑91) 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,1590TC 9.5% 108.8 / 2.7 42intergenic (+382/‑91)PP_2809/PP_2810membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/19);  new base C (2/2);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGG  >  NC_002947/3203012‑3203306
                                                                                                                                                   |                                                                                                                                                   
aTGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGAt                                                                                                                                                     >  1:335013/1‑148 (MQ=255)
     tCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCC‑‑TGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGATCAACTGACCATGTcat                                                                                                                                             >  3:350677/1‑147 (MQ=255)
      ccTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCC‑‑TGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGATCAACTGACCATGTcat                                                                                                                                             <  4:350677/148‑3 (MQ=255)
               ccAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAg                                                                                                                                     >  4:246474/1‑149 (MQ=255)
               ccAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAg                                                                                                                                     >  5:13169/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                                                                                                                                            <  6:163075/103‑4 (MQ=255)
                                                      ataaatTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGATCAACTGACCATGTcatg                                                                                                                                            >  5:163075/1‑100 (MQ=255)
                                                          atTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATg                                                                                                                                 <  2:17257/110‑1 (MQ=255)
                                                          atTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATg                                                                                                                                 >  1:17257/1‑110 (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                                                                                  <  8:151764/148‑1 (MQ=255)
                                                                   gCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCaa                                                                                  <  2:335013/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                                                                                               <  6:31064/117‑1 (MQ=255)
                                                                                     aTCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGa                                                                                               >  5:31064/1‑117 (MQ=255)
                                                                                           cTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGag                                                                                              <  6:154852/112‑1 (MQ=255)
                                                                                           cTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGag                                                                                              >  5:154852/1‑112 (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)
                                                                                                                                             aGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTc                          >  6:181450/1‑130 (MQ=255)
                                                                                                                                             aGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTc                          <  5:181450/130‑1 (MQ=255)
                                                                                                                                              gTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTAtttt      >  6:200144/1‑149 (MQ=255)
                                                                                                                                               ttGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTc     <  1:178846/149‑1 (MQ=255)
                                                                                                                                                  aTCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAgg  >  2:346119/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTTCAGG  >  NC_002947/3203012‑3203306

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