Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,013,587 T→G 12.7% noncoding (344/385 nt) PP_mr44 → CrcZ ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,013,5870TG12.7% 110.5 / 9.5 47noncoding (344/385 nt)PP_mr44CrcZ ncRNA
Reads supporting (aligned to +/- strand):  ref base T (18/23);  new base G (3/3);  total (21/26)
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.60e-01

GCAGGTCAGCGCTCAAAAAAATACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGATCAGCGCCGTGGCTTCGCGCCCCATGC  >  NC_002947/4013446‑4013733
                                                                                                                                              |                                                                                                                                                  
gCAGGTCAGCGCTCAAAAAAATACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTcc                                                                                                                                             <  1:183126/149‑1 (MQ=255)
gCAGGTCAGCGCTCAAAAAAATACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTcc                                                                                                                                             <  2:105284/149‑1 (MQ=255)
      caggccTCAAAAAAATACCTTTGCTCCTGACCCCGCCGGGGGGGCGCCCGCAACACCATGCCCACCCCCCAAAAACAACAAAAGGCCGGGCCCAAAAATAACAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCCGAAGCTCCCCTCa                                                                                                                                        <  1:76879/144‑1 (MQ=255)
        gcgcTCAAAAAAATACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTg                                                                                                                                     <  6:36562/149‑1 (MQ=255)
              aaaaaaaTACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCt                                                                                                                               >  2:673/1‑149 (MQ=255)
              aaaaaaaTACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCt                                                                                                                               >  4:149674/1‑149 (MQ=255)
                         tttGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAg                                                                                                                                       <  4:92295/130‑1 (MQ=255)
                         tttGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAg                                                                                                                                       >  3:92295/1‑130 (MQ=255)
                             ctctTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCAGCCCTCTTcctcct                                                                                                                 >  3:22970/1‑148 (MQ=255)
                              tctTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTc                                                                                                               <  4:322089/149‑1 (MQ=255)
                               ctTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTc                                                                                                               <  4:22970/148‑1 (MQ=255)
                               ctTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTGCGGCTCCCCTCAGGGCTTCCTGCCCTCTTCCGCCTTCt                                                                                                              >  7:281027/1‑149 (MQ=255)
                                ttGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAACGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTc                                                                                                               <  7:132068/147‑1 (MQ=255)
                                  gACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACACCGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTg                                                                                                             <  5:73322/147‑1 (MQ=255)
                                    ccccGCCTGGGGGTCGCCCGCAACACCAGGACAAGTCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTc                                                                                                                                 >  6:293133/1‑125 (MQ=255)
                                    ccccGCCTGGGGGTCGCCCGCAACACCAGGACAAGTCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTc                                                                                                                                 <  5:293133/125‑1 (MQ=255)
                                     cccGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCa                                                                                                                                        <  5:93321/117‑1 (MQ=255)
                                     cccGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCa                                                                                                                                        >  6:93321/1‑117 (MQ=255)
                                      ccGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTtctc                                                                                                       <  5:65984/149‑1 (MQ=255)
                                               ggTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTg                                                                                              <  1:207744/149‑1 (MQ=255)
                                                 tCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATc                                                                                                 >  8:31742/1‑144 (MQ=255)
                                                 tCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATc                                                                                                 <  7:31742/144‑1 (MQ=255)
                                                     ccGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTc                                                                                        >  6:340560/1‑149 (MQ=255)
                                                                ggACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTcca                                                                             >  1:60972/1‑149 (MQ=255)
                                                                   cAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTccacca                                                                          <  5:340560/149‑1 (MQ=255)
                                                                          aCAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTGATCCTGCCCGTCGCCCGTCTGCCGCAGCCACgg                                                                   >  8:253418/1‑149 (MQ=255)
                                                                              aaaCAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATc                                                               <  7:253418/149‑1 (MQ=255)
                                                                               aacaacAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCCGAAGCTCCCCTc                                                                                                                                         <  4:77531/74‑1 (MQ=38)
                                                                               aacaacAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCCGAAGCTCCCCTc                                                                                                                                         >  3:77531/1‑74 (MQ=38)
                                                                                acaacaATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTtc                                                                                 <  5:151363/129‑1 (MQ=255)
                                                                                acaacaATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTtc                                                                                 >  6:151363/1‑129 (MQ=255)
                                                                                  aacaaTAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGATTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTATCCACCAGCAACGGCATCGTGc                                                           >  7:284351/1‑149 (MQ=255)
                                                                                          gCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGcc                                                          <  8:271291/142‑1 (MQ=255)
                                                                                          gCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGcc                                                          >  7:271291/1‑142 (MQ=255)
                                                                                                 ccAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTc                                                                                        <  1:167681/105‑1 (MQ=255)
                                                                                                 ccAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTc                                                                                        >  2:167681/1‑105 (MQ=255)
                                                                                                  cAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGc                                           <  1:673/149‑1 (MQ=255)
                                                                                                   aataatAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGCAGCTTCGGCTCCCCTCAGTGCTTCCTGACCTCTTCCTCCTTCTGCTTCTACTCATCCTGACCTTCGACATTCACCACCAGCAACGGCATCCTGCCCAATACACGCAGGGc                                           >  4:194100/1‑148 (MQ=255)
                                                                                                                      ggCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGATCAgc                        >  7:189293/1‑149 (MQ=255)
                                                                                                                                cTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGATCAGCGCCGTGGCtt              <  7:121107/149‑1 (MQ=255)
                                                                                                                                 tttGAGGGGAGCTTCGTCTCCCCTCAGTGCTTCCAGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTCCTGATCAGCGCCGTGGCTTc             >  1:129020/1‑149 (MQ=255)
                                                                                                                                 tttGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGATCAGCGCCGTGGCTTc             >  3:329466/1‑149 (MQ=255)
                                                                                                                                   tGAGGGGAGCCGAAGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGa                             <  3:104758/131‑1 (MQ=255)
                                                                                                                                   tGAGGGGAGCCGAAGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGa                             >  4:104758/1‑131 (MQ=255)
                                                                                                                                    gAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGATCAGCGCCGTGGCTTcgc           <  3:149674/148‑1 (MQ=255)
                                                                                                                                       gtgAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGGTGATCAGCGCCGTGGCTTCGCGcc        <  2:129020/146‑1 (MQ=255)
                                                                                                                                            cTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGATCAGCGCCGTGGCTTCGCGCCCCATGc  <  8:284351/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GCAGGTCAGCGCTCAAAAAAATACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAA‑AAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGATCAGCGCCGTGGCTTCGCGCCCCATGC  >  NC_002947/4013446‑4013733

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