Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,707,318 A→G 12.0% H355R (CAT→CGT)  srmB → ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,707,3180AG12.0% 75.2 / 2.0 34H355R (CAT→CGT) srmBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base A (10/19);  new base G (2/2);  total (12/22)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

TGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTATTGATGAGTTACGCCCGAAAACGCGTGCGCCAAGCGAAAAGCAGACCGGCAAGCCATCGAAGAAAGTACTGGCTAAA  >  CP009273/2707172‑2707457
                                                                                                                                                  |                                                                                                                                           
tGTTGCCGGGCGCGGCATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATg                                                                                                                                           <  1:405355/149‑1 (MQ=255)
   tGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGAcc                                                                                                                                        <  1:457304/149‑1 (MQ=255)
           gcgGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCGTGACCATctgctg                                                                                                                                >  1:448108/1‑149 (MQ=255)
           gcgGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATctgctg                                                                                                                                >  1:278757/1‑149 (MQ=255)
           gcgGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAAGTTCGCTGGTGGAAGCCCGTGACCACctgctg                                                                                                                                >  1:448112/1‑149 (MQ=255)
               tATCGACATTCCGGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTgg                                                                                                                            <  2:281676/149‑1 (MQ=255)
                             ccGTCAGCCACGTCTTTAACTTCGATCTGCCGCGCAGTGGCGATACTTCTTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGTCCATCTGCTGCTGGGTAAAGTAGGCCGc                                                                                                              <  1:257933/148‑1 (MQ=255)
                                tCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCtat                                                                                                           >  2:9843/1‑149 (MQ=255)
                                  acccACGTGTTGACCTTCGATATGCCGCGCAGTGGCGATCCTTATTCGCTCCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTatat                                                                                                         <  2:281365/147‑1 (MQ=255)
                                       cGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTAAAGAACTATCCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAag                                                                                                    >  1:322240/1‑149 (MQ=255)
                                              cactcggattggCCGCGCAGTGGCGATACGTATTTGCACCGTATCGGACGTACCACGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCCGGACCAGCTTCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAAt                                                                                             <  1:67677/138‑1 (MQ=255)
                                                        tGCCGCGCCTTGGCGCTACTTCTTTGCACCGTTTCGGACGGCCCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTTCAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTcgc                                                                                   <  2:351368/149‑1 (MQ=255)
                                                          ccGCGCAGTGGCGATACTTATTTGCCCCGTATCGGGCGAAACCCACGCGCCTGGCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGt                                                                                 >  1:297743/1‑149 (MQ=255)
                                                            gcgcAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTa                                                                               >  1:93370/1‑149 (MQ=255)
                                                                   ggCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTATtgatga                                                                        <  1:414682/149‑1 (MQ=255)
                                                                      gATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATGTTGCAGAGCCAATTAAAGCTCGCGTTATTGGTGAGtt                                                                     >  1:10939/1‑149 (MQ=255)
                                                                       aTACTTATTTGCACCGTATCGGACGTGGCGCGCGGGCCGGTCGGAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTATTGATGAGTTa                                                                    <  2:294084/149‑1 (MQ=255)
                                                                          cTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTATTGATGAGTTACGc                                                                 <  1:3329/149‑1 (MQ=255)
                                                                          cTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATAGTGAAGAGCCAATTAAAGCTCGCGTTATTGATGAGTTACGc                                                                 <  2:117337/149‑1 (MQ=255)
                                                                                          cggactgcccGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCACGACCCTCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTATTGATGAGTTACGCCCGAAAACGCGTgcgc                                                 <  2:449405/141‑1 (MQ=255)
                                                                                          cGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTATTGATGAGTTACGCCCGAAAACGCGTgcgc                                                 <  1:210086/149‑1 (MQ=255)
                                                                                            gACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTATTGATGAGTTACGCCCGAAAACGCGTGCGCCa                                               <  1:411993/149‑1 (MQ=255)
                                                                                                  ccGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGTGAAAGTAGGTCGCTATATTGATGTGTCAATTCTCGTTCGTGTTATTGATGAGTTACGCCCGAAAACGCGTGCGCCAAGCGaa                                         <  1:67463/149‑1 (MQ=255)
                                                                                                      gcgcgcCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTATTGATGAGTTACGCCCGAAAACGCGTGCGCCAAGCGAAAAGc                                     >  1:408352/1‑149 (MQ=255)
                                                                                                       cgcgcCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTATTGATGAGTTACGCCCGAAAACGCGTGCGCCAAGCGAAAAGcg                                    >  1:193942/1‑148 (MQ=255)
                                                                                                            cGGGCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTATTGATGAGTTACGCCCGAAAACGCGTGCGCCAAGCGAAAAGCAGACCg                               <  1:39205/149‑1 (MQ=255)
                                                                                                               tCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTACCGTAGGCCGCGATATTGGAGAGCCGATTAAAGCTCGCGTTATGGATGAGTAACGCCCGAAAACGCGGGCCCAAAGCGAAAAACAGACCGGCa                            >  2:94376/1‑149 (MQ=255)
                                                                                                                                  tCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTATTGATGAGTTACGCCCGAAAACGCGTGCGCCAAGCGAAAAGCAGACCGGCAAGCCATCGAAGAAAGTACt         <  2:108561/149‑1 (MQ=255)
                                                                                                                                    gcggtgagAAGCCCATGACCATCTGTTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTATTGATGAGTTACGCCCGAAAACGCGTGCGCCAAGCGAAAAGCAGACCGGCAAGCCATCGAAGAAAGTACTgg       <  2:341254/142‑1 (MQ=255)
                                                                                                                                     cTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTATTGATGAGTTACGCCCGAAAACGCGTGCGCCAAGCGAAAAGCAGACCGGCAAGCCATCGACGCGAGTACTAGc      >  1:359053/1‑149 (MQ=255)
                                                                                                                                      tggtggCAGCCCGTGACCATCGGCTGCTGGGTAAGGTAGGCCGCTATATTGAAGCGCCGATTAAAGCTCGCGTTATTGATGAGTTACGCCCGAAAACGCGTGCGCCTAGCGAAAAGCAGACCGGCAAGCCATCGAAGAAAGTACTGGCt     <  2:448112/149‑1 (MQ=255)
                                                                                                                                      tggtggAAGCCCGTGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTATTGATGAGTTACGCCCGAAAACGCGTGCGCCAAGCGAAAAGCAGACCGGCAAGCCATCGAAGAAAGTACTGGCt     <  2:448108/149‑1 (MQ=255)
                                                                                                                                      tgggggAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTATTGATGAGTTACGCCCGAAAACGCGTGCGCCAAGCGAAAAGCAGACCGGCAAGCCATCGAAGAAAGTACTGGCt     <  2:278757/149‑1 (MQ=255)
                                                                                                                                         tggAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGCAGAGCCCATTAAAGCTCGCGTTATTGATGAGTTACGCCCGAAAACGCGTGCGCCAAGCGAAAAGCAGACCGGCAAGCCATCGAAGAAAGTACTGGCTaaa  <  1:295201/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
TGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGCTATATTGAAGAGCCAATTAAAGCTCGCGTTATTGATGAGTTACGCCCGAAAACGCGTGCGCCAAGCGAAAAGCAGACCGGCAAGCCATCGAAGAAAGTACTGGCTAAA  >  CP009273/2707172‑2707457

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