Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,183,620 A→C 8.6% P466P (CCA→CCC)  PP_3680 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,183,6200AC 8.6% 138.0 / 3.6 46P466P (CCA→CCC) PP_3680hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/17);  new base C (2/2);  total (27/19)
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.96e-01

GTTCGAGAAGCCGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTATTCGAC  >  NC_002947/4183474‑4183767
                                                                                                                                                  |                                                                                                                                                   
gttcgAGAAGCCGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCAc                                                                                                                                                    >  1:260230/1‑148 (MQ=255)
    gagaAGCCGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAggg                                                                                                                                               >  2:163321/1‑149 (MQ=255)
        aGCCGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGgc                                                                                                                                           <  2:260230/149‑1 (MQ=255)
          ccGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGgcgc                                                                                                                                         >  3:111075/1‑149 (MQ=255)
          ccGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGgcgc                                                                                                                                         >  3:106494/1‑149 (MQ=255)
            gTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCCCCCCTAGTGGAgcgatc                                                                                                                                       >  1:313678/1‑146 (MQ=255)
            gTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCCCCCCTAGTGGAgcgatc                                                                                                                                       <  2:313678/149‑4 (MQ=255)
                  gCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCaa                                                                                                                                 >  4:250736/1‑149 (MQ=255)
                          cGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCt                                                                                                                         <  3:250736/149‑1 (MQ=255)
                              ctctACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTggc                                                                                                                      <  8:62271/148‑1 (MQ=255)
                                       gCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGaa                                                                                                            >  3:88572/1‑149 (MQ=255)
                                                 tggtACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGt                                                                                                  <  1:163321/149‑1 (MQ=255)
                                                    tACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTc                                                                                                                                   <  1:145440/113‑1 (MQ=255)
                                                    tACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTc                                                                                                                                   >  2:145440/1‑113 (MQ=255)
                                                               ggggATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGgctg                                                                                                                 >  4:82905/1‑120 (MQ=255)
                                                               ggggATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGgctg                                                                                                                 <  3:82905/120‑1 (MQ=255)
                                                                 ggATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGAccc                                                                                  >  1:89639/1‑149 (MQ=255)
                                                                       tGAACAAAACCTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATAGCTACAATAGGCGGGCGCCATGTCAATCACATCTGTCGGCTGCTGAACGCGCTGGGTAAACCTCATATCACCCTGCTCg                                                                            >  4:264386/1‑149 (MQ=255)
                                                                        gaagaaATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTGat                                                                                         >  8:37657/1‑135 (MQ=255)
                                                                         aagaaATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGAc                                                                          >  5:85467/1‑149 (MQ=255)
                                                                          agaaATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGAcc                                                                         <  5:275836/149‑1 (MQ=255)
                                                                                         tCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCAc                                                                                                                                                    >  4:76432/1‑59 (MQ=255)
                                                                                         tCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCAc                                                                                                                                                    <  3:76432/59‑1 (MQ=255)
                                                                                              tCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCaa                                                     <  4:106494/149‑1 (MQ=255)
                                                                                              tCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCaa                                                     <  4:111075/149‑1 (MQ=255)
                                                                                                 ttCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCaa                                                     <  3:119649/146‑1 (MQ=255)
                                                                                                 ttCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCaa                                                     >  4:119649/1‑146 (MQ=255)
                                                                                                 ttCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGt                                                  >  1:139210/1‑149 (MQ=255)
                                                                                                      gCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTg                                             >  2:268010/1‑149 (MQ=255)
                                                                                                            ggAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCg                                       >  7:149298/1‑149 (MQ=255)
                                                                                                                 tCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACCCCTCATATGACCCTGTTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGtgtg                                   >  8:189026/1‑148 (MQ=255)
                                                                                                                 tCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGCCCCTGCTCGACCTAGACGTCGGGCGCCATCAAGGTGGCTGGGGGCGtgtg                                   >  4:10908/1‑148 (MQ=255)
                                                                                                                  cccGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGGCAATCACTTCTGGCGGCTGCGGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGGGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGaa                                 <  2:188066/149‑1 (MQ=255)
                                                                                                                   ccGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAAc                                <  4:88572/149‑1 (MQ=255)
                                                                                                                         aCGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGtt                          <  3:264386/149‑1 (MQ=255)
                                                                                                                          cGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTc                         <  1:268010/149‑1 (MQ=255)
                                                                                                                                 gCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACt                    >  2:315628/1‑147 (MQ=255)
                                                                                                                                    gTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCt               <  2:139210/149‑1 (MQ=255)
                                                                                                                                       catggcGCCCCCCTAGTGGAGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTc                         >  5:36847/6‑136 (MQ=255)
                                                                                                                                       catggcGCCCCCCTAGTGGAGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTc                         <  6:36847/131‑1 (MQ=255)
                                                                                                                                         gATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATg          >  7:99220/1‑149 (MQ=255)
                                                                                                                                             gCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTAtt      <  6:85467/149‑1 (MQ=255)
                                                                                                                                             gCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTAtt      >  3:155777/1‑149 (MQ=255)
                                                                                                                                              cTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACATCATATGACCCTGCTCGACCTAGACGTGGGGCGACATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTATTc     >  4:282598/1‑149 (MQ=255)
                                                                                                                                               tCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTATTc     >  6:153661/1‑148 (MQ=255)
                                                                                                                                                 cACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTATTCGAc  >  5:232156/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTTCGAGAAGCCGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTATTCGAC  >  NC_002947/4183474‑4183767

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