Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,539,664 C→T 8.3% G441D (GGC→GAC)  PP_4872 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,539,6640CT 8.3% 133.1 / 2.5 48G441D (GGC→GAC) PP_4872hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/22);  new base T (2/2);  total (24/24)
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

GTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCA  >  NC_002947/5539525‑5539811
                                                                                                                                             |                                                                                                                                                   
gTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCt                                                                                                                                            >  2:167444/1‑149 (MQ=255)
        tGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACg                                                                                                                                     >  4:13373/1‑148 (MQ=255)
          aGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCAt                                                                                                                                  <  1:206959/149‑1 (MQ=255)
            gTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACGTTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCAtc                                                                                                                                 >  1:60153/1‑149 (MQ=255)
            gTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCAtct                                                                                                                                >  6:79318/1‑149 (MQ=255)
             ttttCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCAtctc                                                                                                                               >  5:177073/1‑149 (MQ=255)
                    aGATACCCGGGAATACGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGAtt                                                                                                                        <  5:79318/149‑1 (MQ=255)
                     gATACCCGGGAATACGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGAttt                                                                                                                       >  4:78038/1‑149 (MQ=255)
                      aTACCCGGGAATACGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTc                                                                                                                      <  6:249774/149‑1 (MQ=255)
                        gcccGGGAATACGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTTCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCAAAGGTCCGAGACAACGCC‑GGTGAAGTCCGGCAACTTTTCACGCATCTCTTCGATTTCGt                                                                                                                    >  4:136132/2‑149 (MQ=255)
                        gcccGGGAATACGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGCCCACGCC‑GGTGAAGTCCGGCCCCTTGTCACGCATCTCTTCGATTTCGt                                                                                                                    >  8:212300/2‑149 (MQ=255)
                              gAATACGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGAtt                                                                                                                        >  6:16061/1‑139 (MQ=255)
                              gAATACGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGAtt                                                                                                                        <  5:16061/139‑1 (MQ=255)
                              gAATACGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGa                                                                                                              >  8:65016/1‑149 (MQ=255)
                                aTACGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATc                                                                                                            <  4:141240/149‑1 (MQ=255)
                                   cGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGTCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGAc                                                                                                         <  1:167444/149‑1 (MQ=255)
                                     cACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTc                                                                                                       <  3:78038/149‑1 (MQ=255)
                                         acgGCTTGCGGCAGTAC‑TTCTCGATTTCGTGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGc                                                                                                   <  8:190354/149‑1 (MQ=255)
                                                   gCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGTCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCt                                                                                         <  3:136132/149‑1 (MQ=255)
                                                   gCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGTCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCt                                                                                         <  7:212300/149‑1 (MQ=255)
                                                   gCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCt                                                                                         >  2:156491/1‑149 (MQ=255)
                                                           ttCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGc                                                                                                                  <  6:32541/117‑1 (MQ=255)
                                                           ttCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGc                                                                                                                  >  5:32541/1‑117 (MQ=255)
                                                           ttCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGgatgat                                                                                  >  7:3429/1‑149 (MQ=255)
                                                               cGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCgg                                                                              >  2:119054/1‑149 (MQ=255)
                                                                         gCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCt                                                                    >  3:256135/1‑149 (MQ=255)
                                                                                    cGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAg                                                         <  1:156491/149‑1 (MQ=255)
                                                                                          ggTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGATCACCGAGCacgaa                                                   <  7:65016/149‑5 (MQ=255)
                                                                                          ggTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGATCACCGAGCacgaa                                                   <  8:3429/149‑5 (MQ=255)
                                                                                           gTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCATGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGATCACCGAGCa                                                       <  4:36734/144‑1 (MQ=255)
                                                                                           gTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGATCACCGAGCa                                                       >  3:36734/1‑144 (MQ=255)
                                                                                           gTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACCCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCTTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAgg                                                  >  4:53729/1‑149 (MQ=255)
                                                                                           gTACATGTTGGCAGTCGGCCAGCCAAGGTCCGAGACCACGCCGGGTGAAGCCAGGCACCTTGTCACGAATCTCTTCGATTTCGTGCAGGATCGAAACGTGAGAGCGGTTATGGATGATGCGGCCTTCGTGCACGGTGATCGAGCAGAAg                                                   >  1:56377/1‑149 (MQ=255)
                                                                                               aTGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACa                                              <  2:60153/149‑1 (MQ=255)
                                                                                               aTGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACa                                              >  7:221455/1‑149 (MQ=255)
                                                                                                  ttGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGc                                            >  3:73289/1‑148 (MQ=255)
                                                                                                              ccGCCAAGGGCCGAGACCACGCC‑GGTGAAGCGCGGCACCTTGTCACGGATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCTGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGcc                               <  4:73289/149‑1 (MQ=255)
                                                                                                                 ccAAGGGCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCTTCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCTGTGATCGAGCAGAAGGt                                                 <  3:185245/128‑1 (MQ=255)
                                                                                                                     gttCCGAGACCACGCC‑GTTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCATGATCGACTCTTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCTGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCatt                        <  2:56377/147‑2 (MQ=255)
                                                                                                                     ggTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCatg                        <  4:256135/149‑1 (MQ=255)
                                                                                                                         cGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCTGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCatgatg                     <  6:177073/148‑1 (MQ=255)
                                                                                                                               cACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGc                                                        >  6:241711/1‑107 (MQ=255)
                                                                                                                               cACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGc                                                        <  5:241711/107‑1 (MQ=255)
                                                                                                                               cACGCC‑GGTGAAGCCCGGCACCTTGGCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCGCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCg              <  3:13373/149‑1 (MQ=255)
                                                                                                                                          aaGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGat    >  1:201639/1‑149 (MQ=255)
                                                                                                                                          aaGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGa     >  2:157333/1‑148 (MQ=255)
                                                                                                                                           aGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGatc   <  1:157333/149‑1 (MQ=255)
                                                                                                                                            gCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGatca  >  1:86318/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTAC‑TTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCCGCCAAGGTCCGAGACCACGCC‑GGTGAAGCCCGGCACCTTGTCACGCATCTCTTCGATTTCGTGCAGGATCGACTCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCA  >  NC_002947/5539525‑5539811

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