Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,643,753 G→A G177G (GGC→GGT rlmN ← dual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,643,7530GA98.5% 216.4 / ‑5.0 66G177G (GGC→GGTrlmNdual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (31/33);  total (31/35)
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.88e-01

GCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAATTTCCGAC  >  NC_000913/2643656‑2643849
                                                                                                 |                                                                                                
gCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGTTTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCgcccatacc                                                                                                <  2:189807/100‑1 (MQ=255)
gCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCgcccatacc                                                                                                <  2:920878/100‑1 (MQ=255)
 cGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCgcccataccc                                                                                               >  2:285601/1‑100 (MQ=255)
    gACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCcatc                                                                                            <  1:1278568/100‑1 (MQ=255)
      cgcgACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCcatcat                                                                                          >  1:805139/1‑100 (MQ=255)
       gcgACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCatcatc                                                                                         <  2:992246/100‑1 (MQ=255)
         gACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCAc                                                                                       <  2:981196/100‑1 (MQ=255)
              ttAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGt                                                                                  >  2:642428/1‑100 (MQ=255)
               tAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGATGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGtt                                                                                 >  2:146012/1‑100 (MQ=255)
                 gACAGGTCAAAATCTAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAg                                                                               <  2:1576992/100‑1 (MQ=255)
                 gACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAg                                                                               >  2:1667186/1‑100 (MQ=255)
                 gACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAg                                                                               >  1:815770/1‑100 (MQ=255)
                    aGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGa                                                                            >  2:181709/1‑100 (MQ=255)
                     ggCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGAt                                                                           >  2:651993/1‑100 (MQ=255)
                       ccAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCg                                                                         >  1:349091/1‑100 (MQ=255)
                       ccAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCg                                                                         <  2:273143/100‑1 (MQ=255)
                        cAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCgg                                                                        <  2:1218195/100‑1 (MQ=255)
                          aaaCCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCCTCATCACTACGTTAGTGATCGGAc                                                                      >  1:468096/1‑100 (MQ=255)
                          aaaCCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGAc                                                                      <  1:285599/100‑1 (MQ=255)
                          aaaCCGAAATCATCAAGCATGATTTCCATCGCAGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGAc                                                                      >  1:745144/1‑100 (MQ=255)
                            aCCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGc                                                                    <  1:598865/100‑1 (MQ=255)
                              cGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTg                                                                  <  2:295288/100‑1 (MQ=255)
                              cGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTg                                                                  >  1:610922/1‑100 (MQ=255)
                                 aaTCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGAcc                                                               <  2:874702/100‑1 (MQ=255)
                                 aaTCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGAcc                                                               <  2:1487800/100‑1 (MQ=255)
                                  atcatcAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCg                                                              <  1:330531/100‑1 (MQ=255)
                                    catcaAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGt                                                            >  1:1410220/1‑100 (MQ=255)
                                    catcaAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGt                                                            >  1:1802450/1‑100 (MQ=255)
                                     atcaAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTg                                                           >  1:6988/1‑100 (MQ=255)
                                     ataaatcATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCAGACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTg                                                           <  2:326354/94‑1 (MQ=255)
                                         agcaccaTTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATTCCCACCCTCACTACGTTAGTGCTCGGAcgctgaccgctgac                                                         <  1:989677/92‑1 (MQ=255)
                                         aGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACtt                                                       <  1:394874/100‑1 (MQ=255)
                                            aTGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAATGATCGGACGCTGACCGGTGACtttt                                                     >  2:1153643/1‑99 (MQ=255)
                                             tGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTtgc                                                   <  1:20425/100‑1 (MQ=255)
                                               aTTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTtgctg                                                 <  2:951700/100‑1 (MQ=255)
                                               aTTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTtgctg                                                 <  1:1735953/100‑1 (MQ=255)
                                               aTTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTtgctg                                                 >  2:132008/1‑100 (MQ=255)
                                                  tCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTgcgc                                              >  1:1531964/1‑100 (MQ=255)
                                                  tCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTgcgc                                              >  2:1506608/1‑100 (MQ=255)
                                                    cATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCcg                                            >  2:480608/1‑100 (MQ=255)
                                                     aTCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCcga                                           >  1:917746/1‑100 (MQ=255)
                                                      tCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGTCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCcgac                                          <  2:456351/100‑1 (MQ=255)
                                                      tCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCcgac                                          <  2:1079243/100‑1 (MQ=255)
                                                            gAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGAtttt                                    >  1:1598939/1‑100 (MQ=255)
                                                            gAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGAtttt                                    >  2:546440/1‑100 (MQ=255)
                                                              acgacgTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTcg                                  <  1:314699/100‑1 (MQ=255)
                                                              acgacgTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTcg                                  <  2:823610/100‑1 (MQ=255)
                                                               cgacgTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTcgc                                 >  1:1679730/1‑100 (MQ=255)
                                                               cgacgTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTcgc                                 >  2:404812/1‑100 (MQ=255)
                                                                gacgTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTcgcc                                >  1:151924/1‑100 (MQ=255)
                                                                 acgTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTcgccg                               <  1:651989/100‑1 (MQ=255)
                                                                  cgTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTcgccgc                              <  1:213455/100‑1 (MQ=255)
                                                                  cgTTGTTCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTcgccgc                              >  2:665644/1‑100 (MQ=255)
                                                                    ttgttCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCAc                            <  1:483240/100‑1 (MQ=255)
                                                                     tgttCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACg                           <  1:592109/100‑1 (MQ=255)
                                                                       ttCAGGTTGAGCAGCTGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGcc                         <  1:380802/100‑1 (MQ=255)
                                                                       ttCAGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGcc                         <  2:5087/100‑1 (MQ=255)
                                                                          aGGTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCcaca                      <  2:1817084/100‑1 (MQ=255)
                                                                           ggTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCacac                     >  1:1171660/1‑100 (MQ=255)
                                                                           ggTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCacac                     >  1:868990/1‑100 (MQ=255)
                                                                            gTTGAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACAcc                    <  1:352586/100‑1 (MQ=255)
                                                                               gAGCAGCGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTgg                 <  1:26765/100‑1 (MQ=255)
                                                                                agcagcGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGc                >  2:88530/1‑100 (MQ=255)
                                                                                agcagcGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCAACGATTTTCGCCGCACGCCACACCTg                  >  2:262859/1‑98 (MQ=255)
                                                                                 gcagcGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGcc               <  2:1794234/100‑1 (MQ=255)
                                                                                 gcagcGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGcc               <  2:917752/100‑1 (MQ=255)
                                                                                   agcGGCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGa             <  1:1012429/100‑1 (MQ=255)
                                                                                       gCTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAAt         >  2:1740243/1‑100 (MQ=255)
                                                                                        cTCGCCCATACCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAAtt        <  2:851994/100‑1 (MQ=255)
                                                                                           gctcatacccatCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAATTTcc     <  2:1753065/97‑1 (MQ=255)
                                                                                           gcccatacccatCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAATTTcc     >  1:477027/1‑100 (MQ=255)
                                                                                            cccatacccatCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAATTTCCg    >  2:734098/1‑100 (MQ=255)
                                                                                            cccatacccatCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAATTTCCg    >  2:734361/1‑100 (MQ=255)
                                                                                            cccatacccatCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAATTTCCg    >  1:986462/1‑100 (MQ=255)
                                                                                              catacccatCATCACTACGTTAGTGATCGTACGCTGACCGGTTACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAATTTCCGac  <  2:488400/100‑1 (MQ=255)
                                                                                                 |                                                                                                
GCGTGACGCGACGTTTAGACAGGCCAAAACCGAAATCATCAAGCATGATTTCCATCGCCGGAACGACGTTGTTCAGGTTGAGCAGCGGCTCGCCCATGCCCATCATCACTACGTTAGTGATCGGACGCTGACCGGTGACTTTTGCTGCGCCGACGATTTTCGCCGCACGCCACACCTGGCCGATAATTTCCGAC  >  NC_000913/2643656‑2643849

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