Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,332,872 C→A R609S (CGC→AGC)  topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,332,8720CA98.6% 242.2 / ‑3.7 71R609S (CGC→AGC) topADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (43/27);  total (43/28)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

AAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAAC  >  NC_000913/1332776‑1332967
                                                                                                |                                                                                               
aaaGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGca                                                                                              >  1:2238355/1‑100 (MQ=255)
   gCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAg                                                                                           >  2:633861/1‑100 (MQ=255)
   gCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAg                                                                                           <  2:2298507/100‑1 (MQ=255)
     tGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTTCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCg                                                                                         <  2:1443601/100‑1 (MQ=255)
     tGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCg                                                                                         <  1:1327879/100‑1 (MQ=255)
      gAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGa                                                                                        >  1:1320448/1‑100 (MQ=255)
             aTCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCgg                                                                                 <  2:1988460/100‑1 (MQ=255)
             aTCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCgg                                                                                 >  1:1509864/1‑100 (MQ=255)
               ccGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCgggg                                                                               <  1:1015404/100‑1 (MQ=255)
               ccGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCgggg                                                                               <  2:1675486/100‑1 (MQ=255)
               ccGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCgggg                                                                               <  1:2042024/100‑1 (MQ=255)
                cGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGt                                                                              >  1:2579487/1‑100 (MQ=255)
                cGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGt                                                                              <  2:365257/100‑1 (MQ=255)
                 ggAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTa                                                                             >  2:771062/1‑100 (MQ=255)
                  gAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTAt                                                                            <  2:421679/100‑1 (MQ=255)
                   aaGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTAtt                                                                           <  1:240639/100‑1 (MQ=255)
                        ggtggtATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTg                                                                      >  2:1300575/1‑100 (MQ=255)
                         gtggtATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTgg                                                                     <  2:2644261/100‑1 (MQ=255)
                         gtggtATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTgg                                                                     <  2:2501655/100‑1 (MQ=255)
                         gtggtATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTgg                                                                     >  2:2439545/1‑100 (MQ=255)
                         gtggtATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTgg                                                                     <  1:2167173/100‑1 (MQ=255)
                          tggtATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGc                                                                    >  2:904019/1‑100 (MQ=255)
                              aTGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGtt                                                                >  1:2656895/1‑100 (MQ=255)
                                 cgcCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTg                                                             >  1:2358896/1‑100 (MQ=255)
                                 cgcCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTg                                                             >  1:1940609/1‑100 (MQ=255)
                                  gcCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTgg                                                            <  2:2404403/100‑1 (MQ=255)
                                   cccGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGc                                                           <  1:2264938/100‑1 (MQ=255)
                                    ccGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCt                                                          <  1:1157673/100‑1 (MQ=255)
                                           aGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCt                                                   >  2:1367591/1‑100 (MQ=255)
                                                gTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCACTgccgc                                              <  1:352384/100‑1 (MQ=255)
                                                 ttCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTgccgcc                                             <  1:904027/100‑1 (MQ=255)
                                                    ngACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGaa                                          >  1:1395296/2‑100 (MQ=255)
                                                       ccAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAaga                                       >  2:1267634/1‑100 (MQ=255)
                                                       ccAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAaga                                       >  1:1736632/1‑100 (MQ=255)
                                                       ccAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATGCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAaga                                       >  2:1785797/1‑100 (MQ=255)
                                                       ccAGCATTGACTGCCCGACTTGTGGTCGAAAAATGCGGCTTAGCACCGCGGGCACCGGGGGATTCCTTGGCTGTTCCGGCTATGAGATgccgc                                              >  2:2235876/1‑93 (MQ=255)
                                                        cAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAagag                                      <  1:1236978/100‑1 (MQ=255)
                                                         aGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGc                                     <  2:1513540/100‑1 (MQ=255)
                                                          gCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCg                                    >  2:2284618/1‑100 (MQ=255)
                                                          gCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCg                                    >  2:2189907/1‑100 (MQ=255)
                                                            aTTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGtt                                  <  1:1226868/100‑1 (MQ=255)
                                                              tGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGc                                <  2:1320438/100‑1 (MQ=255)
                                                                 cTGCCCGACTTGTGGTCGCAACATGGGTCTTAGCACAGCGAACGCCGGGGTATGCCTTGGCTGTTCTGGCAATACGCTGCTGCCGCAACAGCGTTGCaa                              >  2:308980/1‑99 (MQ=255)
                                                                  tGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCaaa                             >  2:2072367/1‑99 (MQ=255)
                                                                   gCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAac                           >  2:927145/1‑100 (MQ=255)
                                                                   gCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAac                           >  1:1909756/1‑100 (MQ=255)
                                                                   gCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAac                           <  1:2257191/100‑1 (MQ=255)
                                                                    cccGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAacc                          >  2:2010526/1‑100 (MQ=255)
                                                                     ccGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAacca                         >  2:1900168/1‑100 (MQ=255)
                                                                     ccGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAacca                         >  1:2688401/1‑100 (MQ=255)
                                                                     ccGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTATGGCTATGCGCTGCCGCCGAAAGAGCGATGCAAAacca                         >  2:685980/1‑100 (MQ=255)
                                                                      cGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAaccac                        <  1:2337382/100‑1 (MQ=255)
                                                                       gACATGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAaccacc                       >  1:297212/1‑100 (MQ=255)
                                                                        aCTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAAccacca                      >  2:161197/1‑100 (MQ=255)
                                                                           tgtgGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTa                   <  2:1823972/100‑1 (MQ=255)
                                                                           tgtgGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTa                   >  1:274088/1‑100 (MQ=255)
                                                                            gtgGTCGCAAAATGGGGATTAGCACAGCGAGTACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTaa                  >  2:485722/1‑100 (MQ=255)
                                                                            gtgGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTaa                  >  1:2691773/1‑100 (MQ=255)
                                                                              ggTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAAcc                >  1:1853626/1‑100 (MQ=255)
                                                                               gTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCt               >  2:682963/1‑100 (MQ=255)
                                                                                   cAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTg           >  1:391799/1‑100 (MQ=255)
                                                                                    aaaaTGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGc          <  2:2238332/100‑1 (MQ=255)
                                                                                      aaTGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCg        <  1:1300585/100‑1 (MQ=255)
                                                                                       aTGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCgg       >  1:1270317/1‑100 (MQ=255)
                                                                                       aTGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCgg       >  2:1789800/1‑100 (MQ=255)
                                                                                        tGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGa      >  1:2385579/1‑100 (MQ=255)
                                                                                         ggggATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGaa     >  1:153604/1‑100 (MQ=255)
                                                                                          gggATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGaaa    >  1:1941550/1‑100 (MQ=255)
                                                                                            gATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAAc  <  1:1367602/100‑1 (MQ=255)
                                                                                            gATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAAc  >  2:161340/1‑100 (MQ=255)
                                                                                             aTTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGTCTATGCGCTGCCGCCGAAACAGCGTTGCAGCACCACCATTTACCTGGTg           >  2:2629895/1‑90 (MQ=255)
                                                                                                |                                                                                               
AAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAAC  >  NC_000913/1332776‑1332967

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