Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,172,502 G→T 8.9% D444Y (GAT→TAT)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,172,5020GT 8.9% 1872.9 / 127.3 557D444Y (GAT→TAT) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (305/201);  new base T (27/23);  total (333/224)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 608 total aligned reads displayed.

AAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGT  >  CP009273/4172436‑4172569
                                                                  |                                                                   
aaaGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCgatga                                                                    >  1:36537337/1‑68 (MQ=255)
aaaGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCgatga                                                                    <  1:28640345/68‑1 (MQ=255)
 aaGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCgatgat                                                                   <  1:20239910/68‑1 (MQ=255)
 aaGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCgatgat                                                                   <  1:13837297/68‑1 (MQ=255)
 aaGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCgatgat                                                                   <  1:30036150/68‑1 (MQ=255)
   gacgacATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGata                                                                  >  1:25634949/1‑67 (MQ=255)
     cgacATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGatat                                                                 <  1:27894287/66‑1 (MQ=255)
     cgacATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCg                                                               >  1:23086383/1‑68 (MQ=255)
      gacATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGa                                                              <  1:33173503/68‑1 (MQ=255)
      gacATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGa                                                              <  1:15743934/68‑1 (MQ=255)
      gacATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGa                                                              <  1:5636267/68‑1 (MQ=255)
       acATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGa                                                              >  1:13631206/1‑67 (MQ=255)
          tcatTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGacca                                                           >  1:27466939/1‑67 (MQ=255)
          tcatTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGacca                                                           >  1:42720441/1‑67 (MQ=255)
          tcatTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGacca                                                           >  1:24922318/1‑67 (MQ=255)
          tcatTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGacac                                                           >  1:30585056/1‑65 (MQ=255)
          tcatTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGacaa                                                           >  1:30274877/1‑65 (MQ=255)
           catTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGaccacc                                                         <  1:9119203/68‑1 (MQ=255)
           catTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGaccacc                                                         <  1:13128929/68‑1 (MQ=255)
            atTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGaccacc                                                         >  1:21517635/1‑67 (MQ=255)
            atTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCt                                                        <  1:26399510/68‑1 (MQ=255)
            atTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCt                                                        <  1:38582953/68‑1 (MQ=255)
            atTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCt                                                        <  1:8187725/68‑1 (MQ=255)
            atTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCt                                                        <  1:8440942/68‑1 (MQ=255)
            atTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCt                                                        <  1:376511/68‑1 (MQ=255)
            atTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCt                                                        <  1:3693369/68‑1 (MQ=255)
            atTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCt                                                        <  1:34407994/68‑1 (MQ=255)
              tGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCg                                                      <  1:3832753/68‑1 (MQ=255)
              tGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCg                                                      >  1:10432334/1‑68 (MQ=255)
              tGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCg                                                      >  1:12489187/1‑68 (MQ=255)
              tGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCg                                                      <  1:40323239/68‑1 (MQ=255)
                 tGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCa                                                   >  1:26593125/1‑68 (MQ=255)
                 tGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCa                                                   >  1:7156485/1‑68 (MQ=255)
                 tGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCa                                                   >  1:41102748/1‑68 (MQ=255)
                 tGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCa                                                   >  1:12649394/1‑68 (MQ=255)
                 tGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCa                                                   >  1:15273927/1‑68 (MQ=255)
                 tGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCa                                                   >  1:24020852/1‑68 (MQ=255)
                  gTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCa                                                   >  1:33384107/1‑67 (MQ=255)
                  gTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCa                                                   >  1:17225269/1‑67 (MQ=255)
                  gTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCa                                                   >  1:3142727/1‑67 (MQ=255)
                    tATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAAcc                                                >  1:29011444/1‑68 (MQ=255)
                     aTGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcg                                               <  1:1414025/68‑1 (MQ=255)
                     aTGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcg                                               <  1:27180850/68‑1 (MQ=255)
                     aTGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcg                                               <  1:16496733/68‑1 (MQ=255)
                     aTGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcg                                               <  1:14793497/68‑1 (MQ=255)
                     aTGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcg                                               >  1:5129374/1‑68 (MQ=255)
                      tGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACcgt                                              <  1:17596347/68‑1 (MQ=255)
                      tGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcgt                                              <  1:4186389/68‑1 (MQ=255)
                      tGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcgt                                              <  1:11193549/68‑1 (MQ=255)
                      tGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcgt                                              <  1:4711821/68‑1 (MQ=255)
                      tGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcgt                                              <  1:22313052/68‑1 (MQ=255)
                      tGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcgt                                              <  1:21204259/68‑1 (MQ=255)
                      tGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcgt                                              <  1:39390908/68‑1 (MQ=255)
                      tGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcgt                                              <  1:7464250/68‑1 (MQ=255)
                      tGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcgt                                              <  1:10080504/68‑1 (MQ=255)
                      tGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcgt                                              <  1:37422212/68‑1 (MQ=255)
                        aaaaaGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcgtcg                                            >  1:35772412/1‑68 (MQ=255)
                         aaaaGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcgtcg                                            <  1:32178307/67‑1 (MQ=255)
                          aaaGCTCATCGATATCCGTAACTGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTa                                          <  1:43525889/68‑1 (MQ=255)
                          aaaGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTa                                          <  1:16059784/68‑1 (MQ=255)
                          aaaGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTa                                          <  1:18041021/68‑1 (MQ=255)
                            aGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTAt                                         >  1:12341760/1‑67 (MQ=255)
                            aGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTAt                                         >  1:28278631/1‑67 (MQ=255)
                            aGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTAt                                         >  1:41347504/1‑67 (MQ=255)
                            aGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTAt                                         >  1:35426776/1‑67 (MQ=255)
                                  tCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTAtccgttcc                                  <  1:31792908/68‑1 (MQ=255)
                                    gATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTAtccgttccgt                                <  1:16856683/68‑1 (MQ=255)
                                    gATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTAtccgttccgt                                >  1:19842235/1‑68 (MQ=255)
                                       atCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTg                              >  1:22751738/1‑67 (MQ=255)
                                       atCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTg                              >  1:15113428/1‑67 (MQ=255)
                                       atCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGACGTATCCGTTCCGTTg                              >  1:11984789/1‑67 (MQ=255)
                                        tCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGc                            >  1:34734574/1‑68 (MQ=255)
                                         ccGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCg                           >  1:7792815/1‑68 (MQ=255)
                                         ccGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCg                           >  1:29879342/1‑68 (MQ=255)
                                            tAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGaaa                        >  1:40887985/1‑68 (MQ=255)
                                              aCGGTAAAGGCGAAGTCGATTATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATg                      >  1:2074462/1‑68 (MQ=255)
                                              aCGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATg                      >  1:18279159/1‑68 (MQ=255)
                                              aCGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATg                      >  1:6682396/1‑68 (MQ=255)
                                              aCGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATg                      >  1:6246844/1‑68 (MQ=255)
                                                ggTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGc                    >  1:600363/1‑68 (MQ=255)
                                                  tAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCg                   <  1:34992394/67‑1 (MQ=255)
                                                      ggCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaa               >  1:2163174/1‑67 (MQ=255)
                                                      ggCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaa               >  1:18946323/1‑67 (MQ=255)
                                                      ggCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaa               >  1:36078080/1‑67 (MQ=255)
                                                      ggCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaa               >  1:36303830/1‑67 (MQ=255)
                                                      ggCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaa               >  1:42425181/1‑67 (MQ=255)
                                                      ggCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaa               >  1:40087152/1‑67 (MQ=255)
                                                       gCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaaa              >  1:25196213/1‑67 (MQ=255)
                                                       gCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaaa              >  1:235749/1‑67 (MQ=255)
                                                        cGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAAcc            >  1:32689671/1‑68 (MQ=255)
                                                        cGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAAcc            >  1:1161027/1‑68 (MQ=255)
                                                        cGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAAcc            >  1:34166844/1‑68 (MQ=255)
                                                        cGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTACCCGTTCCGTTGGCGAAATGGCGGaaaa              <  1:18532114/66‑1 (MQ=255)
                                                           aGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGt         <  1:33720155/68‑1 (MQ=255)
                                                             tCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGtt        >  1:23243330/1‑67 (MQ=255)
                                                             tCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGtt        >  1:26134197/1‑67 (MQ=255)
                                                                 tgatATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCcgcg   >  1:1100076/1‑68 (MQ=255)
                                                                 tgatATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCcgcg   >  1:43090893/1‑68 (MQ=255)
                                                                 tgatATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCcgc    <  1:1935651/67‑1 (MQ=255)
                                                                 tgatATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCcgc    <  1:24539520/67‑1 (MQ=255)
                                                                  |                                                                   
AAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGT  >  CP009273/4172436‑4172569

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