Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,410,649 A→G I124V (ATC→GTC)  dusB → tRNA‑dihydrouridine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,410,6490AG96.0% 257.6 / 0.5 76I124V (ATC→GTC) dusBtRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base G (42/31);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 8.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

CCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCC  >  NC_000913/3410552‑3410745
                                                                                                 |                                                                                                
ccAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTc                                                                                                <  1:2268217/100‑1 (MQ=255)
ccAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTc                                                                                                <  1:1690403/100‑1 (MQ=255)
   gATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCtt                                                                                             >  1:1162742/1‑100 (MQ=255)
     ttattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTAc                                                                                           <  2:579593/100‑1 (MQ=255)
     ttattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTAc                                                                                           <  2:579371/100‑1 (MQ=255)
       attGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCg                                                                                         >  1:1888786/1‑100 (MQ=255)
          gATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAgg                                                                                      >  1:12877/1‑100 (MQ=255)
          gATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAgg                                                                                      >  1:17748/1‑100 (MQ=255)
            tatCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGgtc                                                                                    >  1:1205626/1‑100 (MQ=255)
             atCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGgtcg                                                                                   <  2:1527786/100‑1 (MQ=255)
              tCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGgtcgt                                                                                  >  2:916079/1‑100 (MQ=255)
               naaTATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGgtcgtc                                                                                 >  1:42922/2‑100 (MQ=255)
                   atGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATg                                                                             >  2:1614083/1‑100 (MQ=255)
                    tGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGc                                                                            <  2:1264092/100‑1 (MQ=255)
                    tGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGc                                                                            <  2:74677/100‑1 (MQ=255)
                    tGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGc                                                                            <  1:574588/100‑1 (MQ=255)
                      ggTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAg                                                                          <  2:379523/100‑1 (MQ=255)
                      ggTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAg                                                                          >  1:403882/1‑100 (MQ=255)
                         tGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTgg                                                                       >  1:1279657/1‑100 (MQ=255)
                         tGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTgg                                                                       <  2:1162740/100‑1 (MQ=255)
                         tGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTgg                                                                       <  1:784476/100‑1 (MQ=255)
                           cccGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGAc                                                                     <  1:1658449/100‑1 (MQ=255)
                            ccGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGTAGTGGACg                                                                    >  2:2149652/1‑100 (MQ=255)
                             cGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGt                                                                   >  2:567777/1‑100 (MQ=255)
                                cTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTcc                                                                <  1:509370/100‑1 (MQ=255)
                                  aaaaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTg                                                              >  1:1060520/1‑100 (MQ=255)
                                  aaaaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTg                                                              >  2:152652/1‑100 (MQ=255)
                                    aaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGtt                                                            >  1:2125432/1‑100 (MQ=255)
                                     aaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTa                                                           <  1:1831826/100‑1 (MQ=255)
                                       aGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTAcc                                                         <  1:683455/100‑1 (MQ=255)
                                        gTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTAccc                                                        >  2:1517129/1‑100 (MQ=255)
                                        gTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTAccc                                                        >  2:500397/1‑100 (MQ=255)
                                        gTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTAccc                                                        >  2:1013801/1‑100 (MQ=255)
                                          gAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTg                                                      <  2:844580/100‑1 (MQ=255)
                                            aTCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGaa                                                    >  2:1031673/1‑100 (MQ=255)
                                            aTCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGaa                                                    >  2:1031320/1‑100 (MQ=255)
                                             tCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAg                                                   >  2:1110198/1‑100 (MQ=255)
                                             tCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAg                                                   <  2:1205623/100‑1 (MQ=255)
                                              cGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGa                                                  >  2:2115203/1‑100 (MQ=255)
                                                 aaGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTc                                               <  1:1013803/100‑1 (MQ=255)
                                                     tCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCAc                                           <  1:582708/100‑1 (MQ=255)
                                                      cGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCAcc                                          <  1:916081/100‑1 (MQ=255)
                                                      cGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCAcc                                          >  1:2325518/1‑100 (MQ=255)
                                                       gCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCg                                         >  2:1413145/1‑100 (MQ=255)
                                                          ggTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCt                                      >  2:1399827/1‑100 (MQ=255)
                                                           gTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTg                                     >  1:118601/1‑100 (MQ=255)
                                                              cAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGc                                  >  1:856209/1‑100 (MQ=255)
                                                              cAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGc                                  >  1:1233865/1‑100 (MQ=255)
                                                                gCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCAc                                >  2:76069/1‑100 (MQ=255)
                                                                gCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCAc                                >  2:1415566/1‑100 (MQ=255)
                                                                  ccTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCg                              <  1:500679/100‑1 (MQ=255)
                                                                  ccTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCg                              >  2:2121716/1‑100 (MQ=255)
                                                                   ctctTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCgg                             <  2:492334/100‑1 (MQ=255)
                                                                    tctTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGa                            <  1:1110200/100‑1 (MQ=255)
                                                                       tGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAaca                         <  1:1882039/100‑1 (MQ=255)
                                                                       tGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAaca                         >  2:1582049/1‑100 (MQ=255)
                                                                         cAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACAcc                       <  2:1888779/100‑1 (MQ=255)
                                                                         cAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACAcc                       >  1:1062639/1‑100 (MQ=255)
                                                                          aGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCg                      <  2:17747/100‑1 (MQ=255)
                                                                          aGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCg                      <  2:12877/100‑1 (MQ=255)
                                                                            tACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTTAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTa                    <  2:1279654/100‑1 (MQ=255)
                                                                              cccGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAAc                  <  2:1436545/100‑1 (MQ=255)
                                                                                 ggATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGc               >  2:1458482/1‑100 (MQ=255)
                                                                                  gATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCg              <  1:704274/100‑1 (MQ=255)
                                                                                   aTGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGa             <  1:76071/100‑1 (MQ=255)
                                                                                     tctgTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAag           <  2:556682/97‑1 (MQ=255)
                                                                                     gTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAag           <  1:1399830/100‑1 (MQ=255)
                                                                                       ngTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAaga          >  2:433123/2‑99 (MQ=255)
                                                                                       cGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAagag         >  2:2123465/1‑100 (MQ=255)
                                                                                       cGTTAAATCGGTCCTTACCGAGGTCGTCAAGGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAagag         >  2:1537851/1‑100 (MQ=255)
                                                                                        gTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAgaga        >  1:556683/1‑100 (MQ=255)
                                                                                          tAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGAtt      >  1:462806/1‑100 (MQ=255)
                                                                                            aaTCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGc    >  2:877654/1‑100 (MQ=255)
                                                                                             aTCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGcn   <  2:42920/100‑2 (MQ=255)
                                                                                              tCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGccc  >  1:1068080/1‑100 (MQ=255)
                                                                                              tCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGccc  >  1:1671338/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCC  >  NC_000913/3410552‑3410745

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