Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 990,177 A→T 8.5% S20C (AGC→TGC)  PP_0856 → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947990,1770AT 8.5% 136.6 / 3.5 47S20C (AGC→TGC) PP_0856lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (18/25);  new base T (2/2);  total (20/27)
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.67e-01

AAAAGCCAAGGCTGCGCTGTCCGATGACGCGGCGGTCGGTGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATC  >  NC_002947/990035‑990311
                                                                                                                                              |                                                                                                                                      
aaaaGCCAAGGCTGCGCTGTCCGATGACGCGGCGGTCGGTGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTgcagcagc                                                                                                                                   >  4:56782/1‑148 (MQ=255)
 aaaGCCAAGGCTGCGCTGTCCGATGACGCGGCGGTCGGTGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGcagcaa                                                                                                                                 <  1:332028/149‑1 (MQ=255)
 aaaGCCAAGGCTGCGCTGTCCGATGACGCGGCGGTCGGTGGCCTGCAATTGAAGCTGGATAACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGcagcaa                                                                                                                                 >  2:16214/1‑149 (MQ=255)
     ccAAGGCTGCGCTGTCCGATGACGCGGCGGTCGGTGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCTGctgcaacctg                                                                                                                             >  7:265690/1‑146 (MQ=255)
       aaGGCTGCGCTGTCCGATGACGCGGCGGTCGGTGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCTGctgcaacctgcg                                                                                                                           <  8:265690/149‑6 (MQ=255)
       aaGGCTGCGCTGTCCGATGACGCGGCGGTCGGTGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGcagcaacagc                                                                                                                             >  4:96359/1‑147 (MQ=255)
          gCTGCGCTGTCCGATGACGCGGCGGTCGGTGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCaagaa                                                                                                                        <  4:76167/149‑1 (MQ=255)
              cgcTGTCCGATGACGCGGCGGTCGGTGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGaa                                                                                                                    <  3:96359/149‑1 (MQ=255)
                cTGTCCGATGACGCGGCGGTCGGTGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACt                                                                                                                  <  7:34524/149‑1 (MQ=255)
                      gATGACGCGGCGGTCGGTGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGaa                                                                                                                    <  4:211478/141‑1 (MQ=255)
                      gATGACGCGGCGGTCGGTGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGaa                                                                                                                    >  3:159099/1‑141 (MQ=255)
                      gATGACGCGGCGGTCGGTGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGaa                                                                                                                    <  4:159099/141‑1 (MQ=255)
                      gATGACGCGGCGGTCGGTGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGaa                                                                                                                    >  3:211478/1‑141 (MQ=255)
                      gATGACGAGGCGGTCGGTGGCCTGCAATTGAAGCTGGATGACCTGGCTAAATGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCATGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGccccc                                                                                                            <  5:119807/149‑1 (MQ=255)
                       aTGACGCGGCGGTCGGTGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCg                                                                                                           >  4:100347/1‑149 (MQ=255)
                                    cggtGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGc                                                                                                         >  1:139248/1‑138 (MQ=255)
                                    cggtGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGc                                                                                                         <  2:139248/138‑1 (MQ=255)
                                     ggtGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCaa                                                                                             <  7:201883/149‑1 (MQ=255)
                                                aTTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAgg                                                                                  <  3:100347/149‑1 (MQ=255)
                                                  tGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCa                                                                                        >  2:264357/1‑141 (MQ=255)
                                                  tGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCa                                                                                        <  1:264357/141‑1 (MQ=255)
                                                    aaGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAgt                                                                              <  1:81884/149‑1 (MQ=255)
                                                          gATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAgtggtg                                                                          <  3:56782/147‑1 (MQ=255)
                                                                     tAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCa                                                               >  5:71162/1‑147 (MQ=255)
                                                                     tAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCa                                                               <  6:71162/147‑1 (MQ=255)
                                                                       aaGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTg                                                            >  8:152390/1‑148 (MQ=255)
                                                                        aGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGGGCTGAAGAAGCAGTGGa                                                          <  7:152390/149‑1 (MQ=255)
                                                                         ggggACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAg                                                         <  8:15079/149‑1 (MQ=255)
                                                                          gggACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGTAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGc                                                        >  4:246516/1‑149 (MQ=255)
                                                                           ggACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCt                                                                        >  7:80097/1‑132 (MQ=255)
                                                                           ggACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCt                                                                        <  8:80097/132‑1 (MQ=255)
                                                                                 gTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGa                                                 <  3:184022/149‑1 (MQ=255)
                                                                                             ttGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGc                                      >  3:80614/1‑148 (MQ=255)
                                                                                               ggAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAgtgg                                                                            <  5:80153/108‑1 (MQ=255)
                                                                                               ggAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAgtgg                                                                            >  6:80153/1‑108 (MQ=255)
                                                                                                gAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAgg                                  <  4:80614/149‑1 (MQ=255)
                                                                                                 aaaCATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAggg                                 >  6:130585/1‑149 (MQ=255)
                                                                                                       gcagcaGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAAc                           >  1:143822/1‑149 (MQ=255)
                                                                                                         agcagTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAAcc                          >  5:76446/1‑148 (MQ=255)
                                                                                                          gcagTGCTGACCCTGGCCGTACTGGCCGCAGGTTGTAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTa                        <  3:246516/149‑1 (MQ=255)
                                                                                                            agTGCTGACCCTGGCCGGACGGGCCGCAGGTTCCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGGCGTGCTGCCCAAGTTCACCTAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTaca                      <  7:234345/149‑1 (MQ=255)
                                                                                                                    cgcTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGATGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTg              <  2:224623/147‑1 (MQ=255)
                                                                                                                    cccTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTg              <  7:177285/149‑1 (MQ=255)
                                                                                                                                tGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTaca                      <  4:321560/129‑1 (MQ=255)
                                                                                                                                tGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTaca                      >  3:321560/1‑129 (MQ=255)
                                                                                                                                  gCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATc  >  6:335907/1‑147 (MQ=255)
                                                                                                                                       ctgTTGCTGCTGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATc                                               <  4:111576/95‑1 (MQ=255)
                                                                                                                                       ctgTTGCTGCTGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATc                                               >  3:111576/3‑97 (MQ=255)
                                                                                                                                              |                                                                                                                                      
AAAAGCCAAGGCTGCGCTGTCCGATGACGCGGCGGTCGGTGGCCTGCAATTGAAGCTGGATGACCTGGCTAAAGGGGACGCGTAAGTGATCGGTTGGAAACATGCAGCAGTGCTGACCCTGGCCGTACTGGCCGCAGGTTGCAGCAGCAACAGCAAGAAGGAACTGCCCCCGGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATC  >  NC_002947/990035‑990311

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