Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,182,442 G→C V400L (GTT→CTT)  rpoB → RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,182,4420GC98.1% 173.0 / ‑2.7 53V400L (GTT→CTT) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (32/20);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

AAGCAGCTGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCT  >  NC_000913/4182380‑4182503
                                                              |                                                             
aaGCAGCTGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCtt                                                             >  1:3056331/1‑65 (MQ=255)
 agcagcTGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTg                                                            >  1:3525461/1‑65 (MQ=255)
 agcagcTGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTg                                                            >  1:4195065/1‑65 (MQ=255)
 agcagcTGAAAGCCTGTTCGAGAACCTGTTCTCATCCGAAGACCGTTATGACTTGTCTGCGCTTg                                                            >  1:1324581/1‑65 (MQ=255)
   cagcTGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGggtggt                                                          >  1:1991603/1‑65 (MQ=255)
   cagcTGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGt                                                          <  1:2912251/65‑1 (MQ=255)
   cagcTGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGt                                                          >  1:1585858/1‑65 (MQ=255)
   cagcTGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGt                                                          <  1:4335929/65‑1 (MQ=255)
   cagcTGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGt                                                          <  1:1709577/65‑1 (MQ=255)
   cagcTGAAAGCCTGTTCGAGAACCTGTTCTTCGCCGAAGACCGTTATGACTTGTCTGCGCTTGGt                                                          <  1:1638668/65‑1 (MQ=255)
      cTGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGt                                                       <  1:425952/65‑1 (MQ=255)
           aGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGaa                                                  <  1:217310/65‑1 (MQ=255)
            gCCTGTTCGAGAACTTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGCATGacg                                                 >  1:2487364/1‑63 (MQ=255)
            gCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTTGTCGTATGact                                                 >  1:3653347/1‑63 (MQ=255)
            gCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTTGTCGTATGAAg                                                 >  1:2087195/1‑65 (MQ=255)
            gCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATgcct                                                 >  1:974107/1‑62 (MQ=255)
            gCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGacg                                                 >  1:4140282/1‑63 (MQ=255)
            gCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGCAg                                                 >  1:3447189/1‑65 (MQ=255)
            gCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGCAg                                                 >  1:4338482/1‑65 (MQ=255)
            gCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGCAg                                                 >  1:4345337/1‑65 (MQ=255)
            gCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAg                                                 >  1:1084836/1‑65 (MQ=255)
            gCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAg                                                 >  1:3945027/1‑65 (MQ=255)
            gCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAg                                                 >  1:2145652/1‑65 (MQ=255)
            gCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAg                                                 >  1:3865610/1‑65 (MQ=255)
            gCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAg                                                 >  1:2298444/1‑65 (MQ=255)
            gCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAg                                                 >  1:3540052/1‑65 (MQ=255)
            gCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAg                                                 >  1:4222762/1‑65 (MQ=255)
            gCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAg                                                 >  1:4231282/1‑65 (MQ=255)
            gCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGAATTGTCTGCGCTTGGTCGTATgccg                                                 >  1:966770/1‑62 (MQ=255)
             ccTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAGt                                                >  1:204462/1‑65 (MQ=255)
                     agaACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGGTTGGTCGTATGAAGTTCAACCGt                                        <  1:3850247/65‑1 (MQ=255)
                     agaACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGt                                        <  1:3483566/65‑1 (MQ=255)
                     agaACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGt                                        <  1:2548438/65‑1 (MQ=255)
                     agaACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGt                                        <  1:3408838/65‑1 (MQ=255)
                     agaACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGt                                        <  1:1564238/65‑1 (MQ=255)
                     agaACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGt                                        >  1:2478561/1‑65 (MQ=255)
                     agaACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGt                                        <  1:513895/65‑1 (MQ=255)
                       aaCCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGACGTTCTACCGTtc                                      >  1:878120/1‑65 (MQ=255)
                       aaCCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTtc                                      <  1:3485124/65‑1 (MQ=255)
                         ccTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGACGTTTAACGGTtctc                                    >  1:4059372/1‑65 (MQ=255)
                              tcttctCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTctgctg                               >  1:3887172/1‑65 (MQ=255)
                                ttctCCGAAGACCGTTATGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTgcg                             >  1:1703404/1‑65 (MQ=255)
                                          aCCGTTATGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGGTGCGCGAAGAAATc                   >  1:4117664/1‑65 (MQ=255)
                                          aCCGTTATGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATc                   >  1:4245485/1‑65 (MQ=255)
                                           ccGTTATGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCg                  >  1:3586255/1‑65 (MQ=255)
                                                aTGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGt             <  1:3838009/65‑1 (MQ=255)
                                                aTGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGt             <  1:2185388/65‑1 (MQ=255)
                                                aTGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGt             <  1:2727404/65‑1 (MQ=255)
                                                aTGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGt             <  1:2313899/65‑1 (MQ=255)
                                                 tGACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGtt            <  1:3289619/65‑1 (MQ=255)
                                                  gACTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTc           >  1:3342646/1‑65 (MQ=255)
                                                   cagtGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTcc          <  1:783390/62‑1 (MQ=255)
                                                   aCTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTcc          <  1:887354/65‑1 (MQ=255)
                                                   aCTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTcc          <  1:2930126/65‑1 (MQ=255)
                                                   aCTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTcc          <  1:2804600/65‑1 (MQ=255)
                                                    cTTGTCTGCGCTTGGTCGTATGACGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCg         >  1:4254836/1‑65 (MQ=255)
                                                    cTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCg         >  1:1102734/1‑65 (MQ=255)
                                                    cTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCg         >  1:3886321/1‑65 (MQ=255)
                                                    cTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCg         >  1:3803169/1‑65 (MQ=255)
                                                    cTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCg         >  1:3487322/1‑65 (MQ=255)
                                                    cTTGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCg         >  1:1376893/1‑65 (MQ=255)
                                                      tGTCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGt       >  1:80501/1‑65 (MQ=255)
                                                        tCTGCGCTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTAt     >  1:683123/1‑65 (MQ=255)
                                                           gcgcTTGGTCGTATGAAGTTTAACCGTTGTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCt  >  1:3523906/1‑65 (MQ=255)
                                                           gcgcTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATccc  >  1:4200201/1‑64 (MQ=255)
                                                              |                                                             
AAGCAGCTGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCT  >  NC_000913/4182380‑4182503

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