Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,229,964 G→A 8.8% R82W (CGG→TGG)  PP_3707 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,229,9640GA 8.8% 124.2 / 2.9 46R82W (CGG→TGG) PP_3707hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (16/25);  new base A (2/2);  total (18/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

TCATGACCAGATATCTTCTCAATCAATTCGCAAATGCCGACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCAGACTGACTCACCACCATAAATGATCTGA‑GCTCAGACCAATCGTCGCCATCAATAAAATAAGAATTTCCGCAAA  >  NC_002947/4229825‑4230110
                                                                                                                                           |                                                                                                                                                    
tCATGACCAGATATCTTCTCAATCAATTCGCAAATGCCGACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTgg                                                                                                                                             <  2:334853/149‑1 (MQ=255)
         gATATCTTCTCAATCAATTCGCAAATGCCGACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACctc                                                                                                                                    <  4:280250/149‑1 (MQ=255)
          atatCTTATCAATCAATTCGCAAATGCCGACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACctct                                                                                                                                   >  4:182977/1‑149 (MQ=255)
               ttctCAATCAATTCGCAAATGCCGACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑Cgg                                                                                                                             >  3:172684/1‑149 (MQ=255)
                  tcaatcaaTTCGCAAATGCCGACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAg                                                                                                                          >  2:27276/1‑149 (MQ=255)
                    aatcaatTCGCAAATGCCGACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTa                                                                                                                           <  4:160996/146‑1 (MQ=255)
                    aatcaatTCGCAAATGCCGACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTa                                                                                                                           <  4:125216/146‑1 (MQ=255)
                     atcaatTCGCAAATGCCGACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTc                                                                                                                       <  2:274844/149‑1 (MQ=255)
                      tcaatTCGCAAATGCCGACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCtctct‑c                                                                                                                               <  6:308666/140‑1 (MQ=255)
                      tcaatTCGCAAATGCCGACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCtctct‑c                                                                                                                               >  5:308666/1‑140 (MQ=255)
                         atTCGCTAATGCAGACTTGTTAGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAATTTTTTTGCTCGATCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGa                                                                                                                   <  5:273578/149‑1 (MQ=255)
                         atTCGCAAATGCCGACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGa                                                                                                                   <  4:172684/149‑1 (MQ=255)
                              cAAATGCCGACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGtt                                                                                                              <  5:208477/149‑1 (MQ=255)
                              cAAATGCCGACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGtt                                                                                                              <  2:355343/149‑1 (MQ=255)
                               aaaTGCCGACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTc                                                                                                             <  3:182977/149‑1 (MQ=255)
                                 aTGCCGTCTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTccc                                                                                                           <  1:234017/149‑1 (MQ=255)
                                      gACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGa                                                                                                      >  7:373587/1‑149 (MQ=255)
                                       aCTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCACAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCTAt                                                                                                     >  6:317058/1‑149 (MQ=255)
                                       aCTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTAGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCACAAGCTTGGCAGAACCTCTCA‑CGGTAGTTCCTGAGCGTTCACATCCAt                                                                                                     >  6:317055/1‑149 (MQ=255)
                                           gTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTc                                                                                                  <  2:185277/148‑1 (MQ=255)
                                            ttCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCa                                                                                                 >  3:365292/1‑148 (MQ=255)
                                             tCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATAGTACGAGATGGGATGCAGCCTCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCaa                                                                                                <  5:1039/148‑1 (MQ=255)
                                              cGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCaaa                                                                                               <  3:308025/148‑1 (MQ=255)
                                                     gCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTta                                                                                       <  5:28540/149‑1 (MQ=255)
                                                       gAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTtata                                                                                     >  3:227209/1‑149 (MQ=255)
                                                       gAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTtata                                                                                     <  1:27276/149‑1 (MQ=255)
                                                            aGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATc                                                                                >  2:168453/1‑149 (MQ=255)
                                                                aGCTCGCGACATCCTGGACGTACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCa                                                                            <  8:331686/149‑1 (MQ=255)
                                                                 gCTCGCGACATCCTTGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCACAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCa                                                                            <  5:317055/148‑1 (MQ=255)
                                                                     gcgACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCACAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCTATGTCAAATACGCTTATATAATCTTct                                                                            <  5:317058/144‑2 (MQ=255)
                                                                            ccTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCAGACTGACTcac                                                                 >  2:4660/1‑148 (MQ=255)
                                                                             cTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCAGACTGACTcacca                                                               <  8:55913/149‑1 (MQ=255)
                                                                                               tttAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCAGACTGACTcac                                                                 >  6:44991/1‑129 (MQ=255)
                                                                                               tttAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCAGACTGACTCACCACCATAAATGATCTGA‑GCt                                            >  7:114196/1‑149 (MQ=255)
                                                                                                         ttGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCAGACTGACTCACCACCATAAATGATCTGA‑GCTCAGACCAAt                                   <  4:227209/148‑1 (MQ=255)
                                                                                                          tGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCAGACTGACTCACCACCATAAATGATCTGA‑GCTCAGACCAAtc                                  <  6:251272/148‑1 (MQ=255)
                                                                                                             tCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCAGACTGACTCACCACCATAAATGATCTGA‑GCTCAGACCAAtcgtcg                              >  7:232041/1‑149 (MQ=255)
                                                                                                                aCCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCA‑CGGTAGTTCCGGATCGTTCCCATCGATATCAAAAACGCTTATATAAACTACAGCCAGACTCACAACCATAAATGATCATA‑GCTCAGAACAATCGTCGCCa                           >  3:312708/1‑149 (MQ=255)
                                                                                                                     tCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATAAAATCTTCAGACTGACTCACCACCATAAATGATCTGA‑GCTCAGACCAATCGTCGCCATCaat                      <  8:373587/149‑1 (MQ=255)
                                                                                                                       gTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCAGACTGACTCACCACCATAAATGATCTGACGCTCAGACCAATCGTCGCCATCaata                     <  1:4660/149‑1 (MQ=255)
                                                                                                                             gaTGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCTCCGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCAGACTGACTCACCACCATAAATGATCTGA‑GCTCAGACCAATCGTCGCCATCaataaaat                 >  1:64103/1‑147 (MQ=255)
                                                                                                                                gggATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCAGACTGACTCACCACCATAAATGATCTGA‑GCTCAGAACAATCGTCGCCATCAATAAAATAAGAAt           >  3:363717/1‑149 (MQ=255)
                                                                                                                                   aTGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCAGACTGACTCACCACCATAAATGATCTGA‑GCTCAGACCAATCGTCGCCATCAATAAAATAAGAATTTc        <  2:64103/149‑1 (MQ=255)
                                                                                                                                       aGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCAGACTGACTCACCACCATAAATGATCTGA‑GCTCAGACCAATCGTCGCCATCAATAAAATAAGAATTTCCGc     >  1:180654/1‑148 (MQ=255)
                                                                                                                                         ccGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCAGACTGACTCACCACCATAAATGATCTGA‑GCTCAGACCAATCGTCGCCATCAATAAAATAAGAATTTCCGCaaa  <  8:114196/149‑1 (MQ=255)
                                                                                                                                          cGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCAGACTGACTCACCACCATAAATGATCTGA‑GCTCAGACCAATCGTCGcc                            <  4:85871/122‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
TCATGACCAGATATCTTCTCAATCAATTCGCAAATGCCGACTTGTTCGTACTGGCGAATCAGTGAGCTCGCGACATCCTGGACGAACTTATATGCTTTAAATTTTTTGCTCGACCTATCGTACGAGATGGGATGCAGCCGCAAGCTTGGCAGAACCTCTCT‑CGGTAGTTCCTGAGCGTTCCCATCGATGTCAAATACGCTTATATAATCTTCAGACTGACTCACCACCATAAATGATCTGA‑GCTCAGACCAATCGTCGCCATCAATAAAATAAGAATTTCCGCAAA  >  NC_002947/4229825‑4230110

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