Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,036,818 C→T 9.2% R175C (CGC→TGC)  PP_0895 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,036,8180CT 9.2% 119.4 / 2.4 44R175C (CGC→TGC) PP_0895hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/19);  new base T (2/2);  total (23/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 = 6.21e-01

GCCGAGTTCGATTACTGCGCGTTCTGCGCGGTCGAGCGACCCGGTGCGCACCCGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCGAGCTGCACGAAACCTGGGACGCCTACGCCGACCACCTGCAAGGCCTGTG  >  NC_002947/1036674‑1036965
                                                                                                                                                 |                                                                                                                                                   
gCCGAGTTCGATTACTGCGCGTTCTGCGCGGTCGAGCGACCCGGTGCGCACCCGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGa                                                                                                                                                 >  4:199345/1‑149 (MQ=255)
gCCGAGTTCGATTACTGCGCGTTCTGCGCGGTCGAGCGACCCGGTGCGCACCCGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGa                                                                                                                                                 >  2:2076/1‑149 (MQ=255)
  cGAGTTCGATTACTGCGCGTTCTGCGCGGTCGAGCGACCCGGTGCGCACCCGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGAcc                                                                                                                                               >  1:314178/1‑149 (MQ=255)
     gTTCGATTACTGCGCGTTCTGCGCGGTCGAGCGACGCGGTGCGCACCCGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTgg                                                                                                                                            <  7:291978/149‑1 (MQ=255)
     gTTCGATTACTGCGCGTTCTGCGCGGTCGAGCGACCCGGTGCGCACCCGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTgg                                                                                                                                            <  6:361524/149‑1 (MQ=255)
                 cgcgTTCTGCGCGGTCGAGCGACCCGGTGCGCACCCGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGaaa                                                                                                                                <  3:319688/149‑1 (MQ=255)
                 cgcgTTCTGCGCGGTCGAGCGACCCGGTGCGCACCCGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGa                                                                                                                                  <  1:49884/147‑1 (MQ=255)
                          cgcgGTCGAGCGACCCGGTGCGCACCCGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGa                                                                                                                         <  2:3108/147‑1 (MQ=255)
                                   gCGACCCGGTGCGCACCCGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGt                                                                                                              <  1:320370/149‑1 (MQ=255)
                                       cccGGTGCGCACCCGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCg                                                                                                            >  8:231288/1‑147 (MQ=255)
                                       cccGGTGCGCACCCGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCGTCCTGCACGACCCATCTCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCg                                                                                                            <  7:231288/147‑1 (MQ=255)
                                          ggTGCGCACCCGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAg                                                                                                                     <  1:62126/135‑1 (MQ=255)
                                          ggTGCGCACCCGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAg                                                                                                                     >  2:62126/1‑135 (MQ=255)
                                                 aCCCGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAg                                                                                                >  6:74231/1‑149 (MQ=255)
                                                  cccGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAgg                                                                                               <  5:74231/149‑1 (MQ=255)
                                                            ccAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAg                                                                                                                     >  5:216005/1‑117 (MQ=255)
                                                            ccAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAg                                                                                                                     <  6:216005/117‑1 (MQ=255)
                                                            ccAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGt                                                                                     <  5:202136/149‑1 (MQ=255)
                                                              aGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTa                                                                                    >  1:84659/1‑149 (MQ=255)
                                                                       cAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGcc                                                                           <  2:84659/148‑1 (MQ=255)
                                                                       cAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCt                                                                          <  3:199345/149‑1 (MQ=255)
                                                                        aaaCCCCTGCACGGCTTCTGGCGCACCCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTg                                                                         >  8:296902/1‑149 (MQ=255)
                                                                        aaaCCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGTGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGc                                                                            >  1:286806/1‑146 (MQ=255)
                                                                        aaaCCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGTGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGc                                                                            <  2:286806/146‑1 (MQ=255)
                                                                        aaaCCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTAt                                                                                   >  3:259195/1‑139 (MQ=255)
                                                                        aaaCCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTAt                                                                                   <  4:259195/139‑1 (MQ=255)
                                                                         aaCCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTgg                                                                        <  2:314178/149‑1 (MQ=255)
                                                                                       ttCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTAc                                                          >  4:342407/1‑149 (MQ=255)
                                                                                          tGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGc                                                        >  6:226246/1‑148 (MQ=255)
                                                                                              gcAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGTGCGACCAGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATc                                                    >  3:147401/1‑148 (MQ=255)
                                                                                                      ggCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCGAGCTGCAc                                           <  3:342407/149‑1 (MQ=255)
                                                                                                       gCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCGAGCTGCAc                                           >  6:177570/1‑148 (MQ=255)
                                                                                                          tCCTGCACGAACCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCAGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCGAGCTGCACGaaa                                       >  5:125192/1‑149 (MQ=255)
                                                                                                             tGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCGAGCTGCACGAAACCt                                    >  2:311013/1‑149 (MQ=255)
                                                                                                                cACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCGAGCTGCAc                                           <  8:39971/140‑1 (MQ=255)
                                                                                                                cACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCGAGCTGCAc                                           >  7:39971/1‑140 (MQ=255)
                                                                                                                   gACCCATCGCTGCCTACCACCTACGCCTGGAGCGACCTGGACGAGCAGGAAAGCTCCCACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCAAGCTGCACGAAACCTGGGAc                               >  6:182164/1‑149 (MQ=255)
                                                                                                                        aTCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGATCCGACAAGCTGATGTCGTTCTGGCTCACGGCAGTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCGAGCTGACCGAAACCTGGGACTCCTa                          >  2:332094/1‑149 (MQ=255)
                                                                                                                           gCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCAAGCTGCACGAAACCTGGGACGCCTACg                        <  5:182164/148‑1 (MQ=255)
                                                                                                                               cGTACCACCTACGCCTGGTGCGACCAGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCGAGCTGCACGAAACCTGGGACGCCTACGCCGac                   <  4:147401/149‑1 (MQ=255)
                                                                                                                                        tACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCGAGCTGCACGAAACCTGGGACGCCTACGCCGACCACCTGc            >  6:180121/1‑147 (MQ=255)
                                                                                                                                           gCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCGAGCTGCACGAAACCTGGGACGCCTACGCCGACCACCTGCAAGGc       >  8:114269/1‑149 (MQ=255)
                                                                                                                                             cTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCGAGCTGCACGAAACCTGGGACGCCTACGCCGACCACCTGCAAGGCCt     <  5:177570/149‑1 (MQ=255)
                                                                                                                                             cTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCGAGCTGCACGAAACCTGGGACGCCTACGCCGACCACCTGCAAGGCCt     >  6:199660/1‑149 (MQ=255)
                                                                                                                                                gcgcGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCGAGCTGCACGAAACCTGGGACGCCTACGCCGACCACCTGCAAGGCCtgtg  <  5:199660/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCCGAGTTCGATTACTGCGCGTTCTGCGCGGTCGAGCGACCCGGTGCGCACCCGCGACGGCCAGCCGAGTACAAACCCCTGCACGGCTTCTGGCGCAACCGCGGCTTCCTGCACGACCCATCGCTGCGTACCACCTACGCCTGGCGCGACCTGGACGAGCAGGAAAGCTCCGACAAGCTGATGTCGTTCTGGCTCAAGGCACTGCCGGTATGATCCGCCTGGCAGCCTGCCAGTACGCCATCGAGCTGCACGAAACCTGGGACGCCTACGCCGACCACCTGCAAGGCCTGTG  >  NC_002947/1036674‑1036965

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: