Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,118,718 G→C 10.4% P798P (CCG→CCC)  PP_1884 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,118,7180GC10.4% 141.6 / 7.1 48P798P (CCG→CCC) PP_1884hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/22);  new base C (2/3);  total (23/25)
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.95e-01

GCAGATCGATATCGATGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTACCAATTAC  >  NC_002947/2118580‑2118864
                                                                                                                                          |                                                                                                                                                  
gCAGATCGATATCGATGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTgg                                                                                                                                          >  3:251426/1‑149 (MQ=255)
 cAGATCGATATCGATGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTggg                                                                                                                                         <  1:274961/149‑1 (MQ=255)
    atcgatatcgatGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGgca                                                                                                                                      <  3:22931/149‑1 (MQ=255)
        atatcgatGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCtgccgctgccgc                                                                                                                                             <  7:273903/138‑1 (MQ=255)
        atatcgatGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCtgccgctgccgc                                                                                                                                             >  8:273903/1‑138 (MQ=255)
        atatcgatGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTgggg                                                                                                                                        >  6:39314/1‑143 (MQ=255)
        atatcgatGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTgggg                                                                                                                                        <  5:39314/143‑1 (MQ=255)
          atcgatGAAGAGGGCTTGCCTTAATAGCATGCCGAGGCCGCTAGCCCCGCTCATTTGTTCGTTGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATAGTGCTGCCGCTGCCGCTGGGGCAGcaca                                                                                                                                 <  3:311137/148‑1 (MQ=255)
           tcgatGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCCCAGCGGCAGCGGCAGcacgat                                                                                                                               <  6:107270/149‑4 (MQ=255)
                 aaGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGc                                                                                                                          <  3:128372/148‑1 (MQ=255)
                   gACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCt                                                                                                                        >  8:98835/1‑148 (MQ=255)
                   gACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCt                                                                                                                        <  7:98835/148‑1 (MQ=255)
                   gACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCCCAGCGGCAGCGGCAGcac                                                                                                                                  <  3:333/138‑1 (MQ=255)
                   gACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCCCAGCGGCAGCGGCAGcac                                                                                                                                  >  4:333/1‑138 (MQ=255)
                    aCGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTAt                                                                                                                      >  1:174603/1‑149 (MQ=255)
                      ggCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTgggg                                                                                                                                        >  8:187623/1‑129 (MQ=255)
                      ggCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTgggg                                                                                                                                        <  7:187623/129‑1 (MQ=255)
                      ggCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTAt                                                                                                                      >  5:71799/1‑147 (MQ=255)
                        cAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAgg                                                                                                                  <  2:174603/149‑1 (MQ=255)
                          aGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGcc                                                                                                                <  2:20620/149‑1 (MQ=255)
                           gCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGccc                                                                                                               <  4:251426/149‑1 (MQ=255)
                              ttACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCtt                                                                                                             >  4:128444/1‑148 (MQ=255)
                                  tacCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAg                                                                                                          >  2:177513/1‑147 (MQ=255)
                                    ccAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTTAgaagg                                                                                                      >  3:57084/1‑149 (MQ=255)
                                    ccAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAgc                                                                                                         >  1:47472/1‑146 (MQ=255)
                                    ccAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAgc                                                                                                         <  2:47472/146‑1 (MQ=255)
                                    ccAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCCAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCAATCAGGCCCTTGAccagc                                                                                                      >  6:13540/1‑148 (MQ=255)
                                    ccAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGAGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGACAAGCCTATCAGGCCCTTGAgcagg                                                                                                      >  3:52762/1‑149 (MQ=255)
                                            aCGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGTaa                                                                                              <  3:128444/149‑1 (MQ=255)
                                                cgcgAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGTCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACggg                                                                                          <  1:70531/149‑1 (MQ=255)
                                                          gCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAg                                                                                <  6:71799/149‑1 (MQ=255)
                                                             caTTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTa                                                                             <  7:233816/149‑1 (MQ=255)
                                                                   tGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAAc                                                                       >  7:64678/1‑149 (MQ=255)
                                                                                            tctcCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCg                                                                  <  8:217613/129‑1 (MQ=255)
                                                                                            tctcCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCg                                                                  >  7:217613/1‑129 (MQ=255)
                                                                                            tctcCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCtt                                              >  4:242524/1‑149 (MQ=255)
                                                                                             ctcCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCtt                                              >  4:108521/1‑148 (MQ=255)
                                                                                                ccccGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCTACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTc                                            <  5:13540/147‑1 (MQ=255)
                                                                                                  ccGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCa                                        >  6:105095/1‑149 (MQ=255)
                                                                                                       ttACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCa                                   >  7:27044/1‑149 (MQ=255)
                                                                                                          ccccTCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGtt                                <  4:57084/149‑1 (MQ=255)
                                                                                                          cccATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGtt                                <  4:52762/149‑1 (MQ=255)
                                                                                                          cccATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGt                                 <  3:242524/148‑1 (MQ=255)
                                                                                                             aTCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAg                             >  5:149315/1‑149 (MQ=255)
                                                                                                             aTCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAg                             <  5:105095/149‑1 (MQ=255)
                                                                                                                  gCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCgggg                        <  3:108521/149‑1 (MQ=255)
                                                                                                                              cATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGcc            >  6:313082/1‑149 (MQ=255)
                                                                                                                               aTCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGcc            >  2:283348/1‑148 (MQ=255)
                                                                                                                                  gTGCTGCCCCAGCGGCAGCGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCTTACCTTTCa                                           >  7:112255/1‑114 (MQ=255)
                                                                                                                                  gTGCTGCCCCAGCGGCAGCGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCa                                           <  8:112255/114‑1 (MQ=255)
                                                                                                                                        ccgctgccgctgGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTACCAATTAc  <  1:177513/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
GCAGATCGATATCGATGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTACCAATTAC  >  NC_002947/2118580‑2118864

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