Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,013,589 G→A 10.6% noncoding (346/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,5890GA10.6% 123.3 / 4.9 47noncoding (346/385 nt)PP_mr44CrcZ ncRNA
Reads supporting (aligned to +/- strand):  ref base G (19/23);  new base A (2/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 = 2.31e-01

GCAGGTCAGCGCTCAAAAAAATACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGATCAGCGCCGTGGCTTCGCGCCCCATGCCAGC  >  NC_002947/4013446‑4013737
                                                                                                                                                |                                                                                                                                                    
gCAGGTCAGCGCTCAAAAAAATACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTcc                                                                                                                                                 <  1:183126/149‑1 (MQ=255)
gCAGGTCAGCGCTCAAAAAAATACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTcc                                                                                                                                                 <  2:105284/149‑1 (MQ=255)
      caggccTCAAAAAAATACCTTTGCTCCTGACCCCGCCGGGGGGGCGCCCGCAACACCATGCCCACCCCCCAAAAACAACAAAAGGCCGGGCCCAAAAATAACAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCCGAAGCTCCCCTCa                                                                                                                                            <  1:76879/144‑1 (MQ=255)
        gcgcTCAAAAAAATACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTg                                                                                                                                         <  6:36562/149‑1 (MQ=255)
              aaaaaaaTACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCt                                                                                                                                   >  4:149674/1‑149 (MQ=255)
              aaaaaaaTACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCt                                                                                                                                   >  2:673/1‑149 (MQ=255)
                         tttGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAg                                                                                                                                           >  3:92295/1‑130 (MQ=255)
                         tttGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAg                                                                                                                                           <  4:92295/130‑1 (MQ=255)
                             ctctTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCAGCCCTCTTcctcct                                                                                                                     >  3:22970/1‑148 (MQ=255)
                              tctTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTc                                                                                                                   <  4:322089/149‑1 (MQ=255)
                               ctTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTc                                                                                                                   <  4:22970/148‑1 (MQ=255)
                               ctTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTGCGGCTCCCCTCAGGGCTTCCTGCCCTCTTCCGCCTTCt                                                                                                                  >  7:281027/1‑149 (MQ=255)
                                ttGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAACGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTc                                                                                                                   <  7:132068/147‑1 (MQ=255)
                                  gACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACACCGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTg                                                                                                                 <  5:73322/147‑1 (MQ=255)
                                    ccccGCCTGGGGGTCGCCCGCAACACCAGGACAAGTCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTc                                                                                                                                     >  6:293133/1‑125 (MQ=255)
                                    ccccGCCTGGGGGTCGCCCGCAACACCAGGACAAGTCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTc                                                                                                                                     <  5:293133/125‑1 (MQ=255)
                                     cccGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCa                                                                                                                                            >  6:93321/1‑117 (MQ=255)
                                     cccGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCa                                                                                                                                            <  5:93321/117‑1 (MQ=255)
                                      ccGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTtctc                                                                                                           <  5:65984/149‑1 (MQ=255)
                                               ggTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTg                                                                                                  <  1:207744/149‑1 (MQ=255)
                                                 tCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATc                                                                                                     >  8:31742/1‑144 (MQ=255)
                                                 tCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATc                                                                                                     <  7:31742/144‑1 (MQ=255)
                                                     ccGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTc                                                                                            >  6:340560/1‑149 (MQ=255)
                                                                ggACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTcca                                                                                 >  1:60972/1‑149 (MQ=255)
                                                                   cAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTccacca                                                                              <  5:340560/149‑1 (MQ=255)
                                                                          aCAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTGATCCTGCCCGTCGCCCGTCTGCCGCAGCCACgg                                                                       >  8:253418/1‑149 (MQ=255)
                                                                              aaaCAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATc                                                                   <  7:253418/149‑1 (MQ=255)
                                                                               aacaacAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCCGAAGCTCCCCTc                                                                                                                                             <  4:77531/74‑1 (MQ=38)
                                                                               aacaacAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCCGAAGCTCCCCTc                                                                                                                                             >  3:77531/1‑74 (MQ=38)
                                                                                acaacaATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTtc                                                                                     >  6:151363/1‑129 (MQ=255)
                                                                                acaacaATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTtc                                                                                     <  5:151363/129‑1 (MQ=255)
                                                                                  aacaaTAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGATTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTATCCACCAGCAACGGCATCGTGc                                                               >  7:284351/1‑149 (MQ=255)
                                                                                          gCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGcc                                                              <  8:271291/142‑1 (MQ=255)
                                                                                          gCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGcc                                                              >  7:271291/1‑142 (MQ=255)
                                                                                                 ccAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTc                                                                                            <  1:167681/105‑1 (MQ=255)
                                                                                                 ccAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTc                                                                                            >  2:167681/1‑105 (MQ=255)
                                                                                                  cAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGc                                               <  1:673/149‑1 (MQ=255)
                                                                                                   aataatAAAAAAAGAGCAGGCAACAACGCTTTGAGGGCAGCTTCGGCTCCCCTCAGTGCTTCCTGACCTCTTCCTCCTTCTGCTTCTACTCATCCTGACCTTCGACATTCACCACCAGCAACGGCATCCTGCCCAATACACGCAGGGc                                               >  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)
                                                                                                                                                ggCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGATCAGCGCCGTGGCTTCGCGCCCCATGCCAgc  >  5:266989/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCAGGTCAGCGCTCAAAAAAATACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGATCAGCGCCGTGGCTTCGCGCCCCATGCCAGC  >  NC_002947/4013446‑4013737

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: