Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,125,401 A→C 9.1% D65A (GAT→GCT)  PP_4512 → DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,125,4010AC 9.1% 117.8 / 3.0 42D65A (GAT→GCT) PP_4512DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base A (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

CTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTT‑CGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCTTTACGGTGATCCTCGAAAGCCTGC  >  NC_002947/5125255‑5125544
                                                                                                                                                   |                                                                                                                                                
cttcGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTGGGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGAt                                                                                                                                                 >  1:278412/1‑149 (MQ=255)
cttcGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTT‑CGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGCTg                                                                                                                                                >  8:269627/1‑149 (MQ=255)
          aCCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTT‑CGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGcgccg                                                                                                                                         >  3:102482/1‑146 (MQ=255)
          aCCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTT‑CGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGcgccg                                                                                                                                         <  4:102482/146‑1 (MQ=255)
           ccGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTT‑CGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCct                                                                                                                                      <  1:113510/148‑1 (MQ=255)
                gTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTT‑CGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAG‑CAGGCGGCGCGCATCCTGgctgct                                                                                                                               <  6:121758/149‑1 (MQ=255)
                     cGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTT‑CGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAG‑CAGGCGGCGCGCATCCTGGCTGCTcag                                                                                                                            <  4:262404/147‑3 (MQ=255)
                      ggCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTT‑CGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCggtgg                                                                                                                         >  8:66865/1‑149 (MQ=255)
                      ggCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTT‑CGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCggtg                                                                                                                          <  7:7948/148‑1 (MQ=255)
                      ggCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTT‑CGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCggtg                                                                                                                          >  8:7948/1‑148 (MQ=255)
                      ggCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTT‑CGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCggtg                                                                                                                          <  1:79553/148‑1 (MQ=255)
                       gCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTT‑CGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCggtggt                                                                                                                        <  8:119731/149‑1 (MQ=255)
                                ccaccaGCGCCAGCCCGGCGATCATCGTCTT‑CGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGctgctg                                                                                                               <  7:66865/149‑1 (MQ=255)
                                                ggCGATCATCGTCTT‑CGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCa                                                                                               <  5:282519/149‑1 (MQ=255)
                                                     tcatcGTCTT‑CGGGCGCGCCACCTTCGCGGTACTCGCGCTGGCGCTGTTCGCCAGCCTGGCCTGTGCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGTCCGGCCACTGGGTCAGcttc                                                                                          <  6:159398/149‑1 (MQ=255)
                                                       atcGTCTT‑CGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGcttct                                                                                         <  2:278412/148‑1 (MQ=255)
                                                                         cACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTgg                                                                       >  6:111220/1‑149 (MQ=255)
                                                                                 ccGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGt                                                                >  8:297472/1‑148 (MQ=255)
                                                                                         gcgcTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGc                                                       >  4:328353/1‑149 (MQ=255)
                                                                                             tGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCAcc                                                   >  2:157349/1‑149 (MQ=255)
                                                                                                cccTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTc                                                <  5:111220/149‑1 (MQ=255)
                                                                                                 ccTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGcc                                                      >  5:354155/1‑142 (MQ=255)
                                                                                                 ccTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGcc                                                      <  6:354155/142‑1 (MQ=255)
                                                                                                        gCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGc                                        >  8:75659/1‑149 (MQ=255)
                                                                                                           aGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTggg                                                                      >  8:288798/1‑116 (MQ=255)
                                                                                                           aGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTggg                                                                      <  7:288798/116‑1 (MQ=255)
                                                                                                            gCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGc                                    >  1:213855/1‑149 (MQ=255)
                                                                                                                      tCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAg                                     >  4:41971/1‑138 (MQ=255)
                                                                                                                      tCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAg                                     <  3:41971/138‑1 (MQ=255)
                                                                                                                        cGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCtt                         >  8:136254/1‑148 (MQ=255)
                                                                                                                                 aaCGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGCTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCTTTACGGTGATc               >  5:28393/1‑149 (MQ=255)
                                                                                                                                  aCGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCTTTACGGTGATcc              <  3:328353/149‑1 (MQ=255)
                                                                                                                                    gCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTNGGATTCGCCNGCTTCCCGGCCTTTACGGTGATCCTc            >  8:87100/1‑149 (MQ=255)
                                                                                                                                    gCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCTTTACGGTGATCCTc            <  7:136254/149‑1 (MQ=255)
                                                                                                                                    gCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCTTTACGGTGATCCTc            >  8:82384/1‑149 (MQ=255)
                                                                                                                                       cgAGCAG‑CAGGCGGCGCGCATCCTGGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCTTTACGGTGATCCTCGa          >  5:172977/2‑148 (MQ=255)
                                                                                                                                         agcagcCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCTTTACGGCGATCCTCGaaa        <  7:269627/148‑1 (MQ=255)
                                                                                                                                          gcagcCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCTTTACGGTGATCCTCGAAAGc      <  7:87100/149‑1 (MQ=255)
                                                                                                                                           cagcCAGGAGGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCTTTACGGTGATCCTCGAAAGc      <  7:82384/148‑1 (MQ=255)
                                                                                                                                            agcCAGGATGCGAGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCTTTACGGTGATCCTCGAAAGCct    >  4:90514/1‑149 (MQ=255)
                                                                                                                                             gcCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCTTTACGGTGATCCTCGAAAGCctg   <  7:75659/149‑1 (MQ=255)
                                                                                                                                              ccAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCTTTACGGTGATCCTCGAAAGCctgc  >  3:297205/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTT‑CGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCTTTACGGTGATCCTCGAAAGCCTGC  >  NC_002947/5125255‑5125544

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑