Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,391,329 A→G 8.5% E82E (GAA→GAG)  PP_5707 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,391,3290AG 8.5% 128.1 / 2.9 47E82E (GAA→GAG) PP_5707hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/19);  new base G (2/2);  total (26/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 = 9.39e-01

CAAAACGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATCGCCATCA  >  NC_002947/5391186‑5391475
                                                                                                                                               |                                                                                                                                                  
cAAAACGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCgg                                                                                                                                               <  4:108362/149‑1 (MQ=255)
 aaaaCTCCCGCGAACAGGCCTCGGAACTCATGGGCTATCCGCGGGACCAGACCAACGACGGACTGATGGAAGACAAGCCCACACTGATCTAGCCGGCGCTCTACATCTGTGCGTTGGGCAAAGCACTACTCGATGATGCCGAACTCgt                                                                                                                                               >  8:257801/1‑147 (MQ=255)
    aCGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACATGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCat                                                                                                                                            >  6:273393/1‑148 (MQ=255)
        ccGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGc                                                                                                                                       >  1:145435/1‑149 (MQ=255)
              aCAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGaa                                                                                                                                 <  7:69824/149‑1 (MQ=255)
                    cTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACg                                                                                                                            >  7:124625/1‑148 (MQ=255)
                     tCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCa                                                                                                                          >  6:280830/1‑149 (MQ=255)
                      cGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCAc                                                                                                                         <  2:301758/149‑1 (MQ=255)
                        gAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACcg                                                                                                                       <  1:212845/149‑1 (MQ=255)
                             gATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAGTTCGGcatcatc                                                                                                                                        >  4:132121/1‑126 (MQ=255)
                             gATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAGTTCGGcatcatc                                                                                                                                        <  3:132121/127‑2 (MQ=255)
                               tGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGaa                                                                                                                                                    >  4:72106/1‑113 (MQ=255)
                               tGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGaa                                                                                                                                                    <  3:72106/113‑1 (MQ=255)
                                 ggCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCtt                                                                                                               >  3:292260/1‑148 (MQ=255)
                                      tGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAggg                                                                                                         <  5:300385/149‑1 (MQ=255)
                                                     cATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGc                                                                                          >  4:142306/1‑149 (MQ=255)
                                                                 gATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAgggg                                                                                                        >  2:280414/1‑123 (MQ=255)
                                                                 gATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAgggg                                                                                                        <  1:280414/123‑1 (MQ=255)
                                                                  aTGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATTCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCtttaa                                                                             >  8:230649/1‑147 (MQ=255)
                                                                   tGGAAGACAAGCCCAAACTGATGTGGGCGGGGCTCTACCTCTGCGCTTTGGGCAAAGCGCTACTGTATGATGCCGAACTCGGCATGATGCGTTGAATAACGCACCGCGATGCTTAAAGGGGGCATCGTCAGCCGCTTGGCCCCTTTTAc                                                                            <  7:257801/149‑1 (MQ=255)
                                                                   tGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTa                                                                             <  4:292260/148‑1 (MQ=255)
                                                                    ggAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACg                                                                           >  2:257665/1‑149 (MQ=255)
                                                                    ggAAGACAAACCCAAACTGATATGCGCGGCGCACTACCTCTGTGCTTTCGGCAAATCGCTACTGGATGATGACGAACTCGGCATGATGCGTAGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTTGCCCCTTATACg                                                                           >  2:164193/1‑149 (MQ=255)
                                                                               ccAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATc                                                                 >  4:170590/1‑148 (MQ=255)
                                                                                cAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGGCGCTTTTACGCGCCACGAtaga                                                               >  8:125194/1‑146 (MQ=255)
                                                                                       aTCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAg                                                        >  2:5616/1‑149 (MQ=255)
                                                                                            ggcggcGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAg                                                        <  1:5616/144‑1 (MQ=255)
                                                                                              cggcgCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCg                                                 <  5:273393/149‑1 (MQ=255)
                                                                                                    aCTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCGGTCttt                                           <  7:230649/149‑1 (MQ=255)
                                                                                                        ccTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCGGTCTTTCTCa                                       <  3:170590/149‑1 (MQ=255)
                                                                                                           ctGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCg                                                                                           >  7:302035/1‑94 (MQ=255)
                                                                                                           ctGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCg                                                                                           <  8:302035/94‑1 (MQ=255)
                                                                                                            tgtgCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCtt                                                                                >  8:238397/1‑104 (MQ=255)
                                                                                                            tgtgCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCtt                                                                                <  7:238397/104‑1 (MQ=255)
                                                                                                                 tttGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTc                              >  7:98366/1‑149 (MQ=255)
                                                                                                                     gtcAAAGCGCTACTGGATGATGCCGAACTCGGCATTATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCCTCGTCAGCCTCTTGGCCCATTTTACGCCACACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCa                           <  1:164193/146‑1 (MQ=255)
                                                                                                                      gCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCAtt                         >  3:269206/1‑149 (MQ=255)
                                                                                                                        aaaGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTa                       >  6:115023/1‑149 (MQ=255)
                                                                                                                           gcgcTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATgc                    <  4:269206/149‑1 (MQ=255)
                                                                                                                                 cTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGcgccg                >  8:229118/1‑147 (MQ=255)
                                                                                                                                  cgcatcatgCCGAGTTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCgcag             >  5:235656/7‑149 (MQ=255)
                                                                                                                                     atgatgCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCAt          >  5:8115/1‑149 (MQ=255)
                                                                                                                                       gatgCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATCg        >  1:189266/1‑149 (MQ=255)
                                                                                                                                        atgCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATCCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATCGc       <  5:280830/149‑1 (MQ=255)
                                                                                                                                         tgCCGAGTTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATCGcc      <  6:235656/149‑1 (MQ=255)
                                                                                                                                         tgCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATCGc       >  7:38655/1‑148 (MQ=255)
                                                                                                                                            cGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCg                                                 >  4:20985/1‑103 (MQ=255)
                                                                                                                                            cGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCg                                                 <  3:20985/103‑1 (MQ=255)
                                                                                                                                             gAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGGCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATCGCCATCa  <  8:124625/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CAAAACGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATCGCCATCA  >  NC_002947/5391186‑5391475

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