Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 894,431 Δ1 bp 13.8% coding (458/1128 nt) rlmC → 23S rRNA m(5)U747 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273894,4300A.13.8% 73.3 / 7.6 29coding (457/1128 nt)rlmC23S rRNA m(5)U747 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (10/15);  new base . (2/2);  total (12/17)
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.57e-01

ATCGCCCGAGCGGGGTTAACGCCCTACAACGTGGCGCGTAAGCGTGGCGAACTGAAATACATTCTGCTGACTGAAAGCCAGAGTGATGGAGGCATGATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTACACGAACAACTACCGCAGCTGAAAGTTATTACCGTCAATATTCAGCCGGTACATATGGCGATTATGGAAGGGGAGACGGAGATCTACCTGACCGAACAACAGGCACTGGCGGAGCGTTTTAA  >  CP009273/894286‑894569
                                                                                                                                                |                                                                                                                                           
aTCGCCCGAGCGGGGTTAACGCCCTACAACGTGGCGCGTAAGCGTGGCGAACTGAAATACATTCTGCTGACTGAAAGCCAGAGTGATGGAGGCATGATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGgc                                                                                                                                         >  1:217566/1‑149 (MQ=255)
    cccGAGCGGGGTTAACGCCCTACAACGTGGCGCGTAAGCGTGGCGAACTGAAATACATTCTGCTGACTGAAAGCCAGAGTGATGGAGGCATGATGCTGCGCTTTGCACTGCGTTCTGATACCAAGCTGGCGCAACGGCGTAAg                                                                                                                                           >  2:298457/1‑143 (MQ=255)
     ccGAGCGGGGTTAACGCCCTACAACGTGGCGCGTAAGCGTGGCGAACTGAAATACATTCTGCTGACTGAAAGCCAGAGTGATGGAGGCATGATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGc                                                                                                                                    <  1:67999/149‑1 (MQ=255)
         gCGGGGTTAACGCCCTACAACGGGGCGCGTAAGCGTGGCGAACTGAAATACATTCTGCTGACTGAAAGCCAGAGTGATGGAGGCATCATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCt                                                                                                                                      <  2:191424/143‑1 (MQ=255)
                  aCGCCCTACAACGTGGCGCGTAAGCGTGGCGAACTGAAATACATTCTGCTGACTGAAAGCCAGAGTGATGGAGGCATGATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAGGCGCTGCCGTGGTTACACGaa                                                                                                                      >  1:105180/1‑149 (MQ=255)
                  aCGCCCTACAACGTGGCGCGTAAGCGTGGCGAACTGAAATACATTCTGCTGACTGAAAGCCAGAGTGATGGAGGCATGATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAGGCGCTGCCGTGGTTACACGaa                                                                                                                      >  1:105019/1‑149 (MQ=255)
                    gCCCTACAACGTGGCGCGTAAGCGTGGCGAACTGAAATACATTCTGCTGACTGAAAGCCAGAGTGATGGAGGCATGATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCGGCCGTGGttgc                                                                                                                           >  2:272550/1‑141 (MQ=255)
                    gCCCGACAACGTGGCGCGTAAGCGTGGCGAACTGAAATACATTCTGCTGACTGAAAGCCAGAGTGATGGAGGCATGATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTACACGaac                                                                                                                     <  1:287714/149‑1 (MQ=255)
                     cccTACAACGTGGCGCGTAAGCGTGGCGAACTGAAATACATTCTGCTGACTGAAAGCCAGAGTGATGGAGGCATGATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTaca                                                                                                                          <  2:24591/143‑1 (MQ=255)
                         aCAACGTGGCGCGTAAGCGTGGCGAACTGCAATACATTCTGCTGACTGAAAGCCAGAGGGATGGAGGCATGATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTGCACgga                                                                                                                      >  2:210471/1‑141 (MQ=255)
                             cGTGGCGCGTAAGCGTGGCGAACTGAAATACATTCTGCTGACTGAAAGCCAGAGTGATGGAGGCATGATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTACACGAACAACTACCGc                                                                                                            <  1:200624/149‑1 (MQ=255)
                               tGGCGCGTAAGCGTGGCGAACTGAAATACATTCTGCTGACTGAAAGCCAGAGTGATGGAGGCATGATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTACACGAACAACTACCGCAg                                                                                                          <  1:62442/149‑1 (MQ=255)
                                                                    gactgcACGCCAGAGTGACGGAGGCATTATGCTGCGCTGTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTACACGAACAACTACCGCAGCTGAAAGTTATTACCGTCAATATTCAGCCgg                                                                           <  2:281182/143‑1 (MQ=255)
                                                                    gactgaAAGCCAGAGTGATGGAGGCATGATGCGGCGCTTTTTACGGCGTTCTGATACCAAGCTGGCGCATCTGCGTAAGGCGCTGCCGTGGTTACACGAACAACTACCGCAGCTGAAAGTTATTACCGTCAATATTCAGCCgg                                                                           <  2:129410/143‑1 (MQ=255)
                                                                     actgaAAGCCAGAGTGATGGAGGCATGATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTACACGAACAACTACCGCAGCTGAAAGTTATTACCGTCAATATTCAGCCGGTACatat                                                                    >  1:191972/1‑149 (MQ=255)
                                                                       tgaAAGCCAGAGTGATGGAGGCATGATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTACACGAACAACTACCGCAGCTGAAAGTTATTACCGTCAATATTCAGCCGGTACATATgg                                                                  >  1:322672/1‑149 (MQ=255)
                                                                             ccAGAGGTATGGAGGCATGATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTACACGAACAACTACCGCAGCTGAAAGTTATTACCGTCAATATTCAGCCGGTACATATgg                                                                  <  2:198389/143‑1 (MQ=255)
                                                                              cAGAGTGATGGAGGCATGATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTACACGAACAACTACCGCAGCTGAGAGTTATTACCGTCAATATTCAGCCGGTACATATGGCGATTAt                                                           >  1:256013/1‑149 (MQ=255)
                                                                              cAGAGTGAGGGAGGCATGATGCTGTAGTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGATACACGAACAACTACCGCAGCTGAAAGTTATTACCGTCAATATTCAGCCGGTACATATGGCGATTAt                                                           <  1:201261/149‑1 (MQ=255)
                                                                                          gggatgatgCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTACACGAACAACTACCGCAGCTGAAAGTTATTACCGTCAATATTCCGCCGGTACATATGGCGATTATGGAAgg                                                     <  2:172357/140‑1 (MQ=255)
                                                                                           gcaggctgctgcGCTTTGGTCGGCGTTCTAGTACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTACACGAACAACTACCGCAGCGGAAAGTTATGACCGTCAATATTGAGCCGGTACATATGGCGATTATGGAAggg                                                    <  2:125845/137‑1 (MQ=255)
                                                                                                        ttttGTACTCGCTTCGGATACCAAGCTGGCGTAACTGCGTAAGTCGCTGCCGTGCTTACACGAACAACTCTCGCACCTGAAAGTTATTACCGTCAATATTCAGCCGGTACATATGGCCATTATGGAAGGGTAGACGGAGATCt                                       <  2:23068/142‑1 (MQ=255)
                                                                                                        cTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTACACGAACAACTACCGCAGCTGAAAGTTATTACCGTCAATATTCAGCCGGTACATATGGCGATTATGGAAGGGGAGACGGAGATCt                                       <  2:172527/143‑1 (MQ=255)
                                                                                                                    ttCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTACACGAACAACTACCGCAGCTGAAAGTTATTACCGTCAATATTCAGCCGGTACATATGGCGATTATGGTAGGGGCGACGGGGATCTCCCTGACCGaac                           >  2:265670/1‑143 (MQ=255)
                                                                                                                        gATACCAAGCTGGCGCAACTGCGTAGGCGCTGCCGTGGTTACACGAACAACTACCGCAGCTGAAAGTTATTACCGTCAATATTCAGCCGGTACATATGGCGATTATGGAAGGGGAGACGGAGATCTACCTGACCGAACAACAg                      <  2:92718/143‑1 (MQ=255)
                                                                                                                        gATACCAAGCTGGCGCAACTGCGTAGGCGCTGCCGTGGTTACACGAACAACTACCGCAGCTGAAAGTTATTACCGTCAATATTCAGCCGGTACATATGGCGATTATGGAAGGGGAGACGGAGAGCTACCTGACCGAACAACAg                      <  2:92858/143‑1 (MQ=255)
                                                                                                                         gtACCAAGCTGGCGCAACTGCGTAAGGCGCGGCCGTGGTTACGCGAACAACTACCGCAGCTGAAAGTTATTACGGTCAAGATTCAGCCGGTACATATGGCGATTATGGAAGGGGAGACGGAGACCTACCTGACCGAACAACAg                      >  2:263228/2‑143 (MQ=255)
                                                                                                                          tACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTACACGAACAACTACCGCAGCTGAAAGTTATTACCGTCAATATTCAGCCGGTACATATGGCGATTATGGAAGGGGAGACGGAGATCTACCTGACCGAACAACAgg                     <  2:290884/143‑1 (MQ=255)
                                                                                                                                             cGTAAGGCGCTGCCGTGGTTACACGAACAACTACCGCAGCTGAAAGTTATTACCGTCAATATTCAGCCGGTACATATGGCGATTATGGAAGGGGAGACGGAGATCTACCTGACGGAACAACAGGCACTGGGGGGGGGTTTTaa  >  2:52740/1‑143 (MQ=255)
                                                                                                                                                |                                                                                                                                           
ATCGCCCGAGCGGGGTTAACGCCCTACAACGTGGCGCGTAAGCGTGGCGAACTGAAATACATTCTGCTGACTGAAAGCCAGAGTGATGGAGGCATGATGCTGCGCTTTGTACTGCGTTCTGATACCAAGCTGGCGCAACTGCGTAAGGCGCTGCCGTGGTTACACGAACAACTACCGCAGCTGAAAGTTATTACCGTCAATATTCAGCCGGTACATATGGCGATTATGGAAGGGGAGACGGAGATCTACCTGACCGAACAACAGGCACTGGCGGAGCGTTTTAA  >  CP009273/894286‑894569

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: