Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 3,208,093 T→C F563S (TTC→TCC)  rpoD → RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,208,0930TC98.2% 180.3 / NA 54F563S (TTC→TCC) rpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/30);  minor base G (0/1);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTT  >  CP009273/3208028‑3208151
                                                                 |                                                          
tgcgGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCACAAGTTCTGCGTATGCGTTCCg                                                          <  1:824348/68‑1 (MQ=255)
 gcgGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCgg                                                         <  1:2961336/68‑1 (MQ=255)
  cgGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGt                                                        <  1:3810836/68‑1 (MQ=255)
       aCGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGa                                                   <  1:3098139/68‑1 (MQ=255)
       aCGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGa                                                   <  1:2996067/68‑1 (MQ=255)
       aCGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGa                                                   <  1:2868888/68‑1 (MQ=255)
       aCGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGa                                                   <  1:2001990/68‑1 (MQ=255)
       aCGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGa                                                   <  1:1746134/68‑1 (MQ=255)
         gCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGata                                                 <  1:2918516/68‑1 (MQ=255)
         gCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGata                                                 <  1:182849/68‑1 (MQ=255)
         gCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGATCTGCGTATGCGTTCCGGTATCGata                                                 <  1:2617627/68‑1 (MQ=255)
           acgacgTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATg                                               >  1:53660/1‑68 (MQ=255)
           acgacgTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATg                                               >  1:2736203/1‑68 (MQ=255)
              acgTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAac                                            >  1:3368829/1‑68 (MQ=255)
                 tGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACAcc                                         >  1:3850375/1‑68 (MQ=255)
                 tGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACAcc                                         <  1:2228762/68‑1 (MQ=255)
                       ctggcCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTac                                   >  1:429487/1‑68 (MQ=255)
                         ggcCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTacac                                 <  1:775215/68‑1 (MQ=255)
                          gcCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACg                                <  1:4229383/68‑1 (MQ=255)
                          gcCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACg                                <  1:3590073/68‑1 (MQ=255)
                          gcCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGGTCCGGTATCGATATGAACACCGACTACACg                                <  1:3032927/68‑1 (MQ=255)
                            cTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCt                              <  1:1059272/68‑1 (MQ=255)
                               aCCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGga                           <  1:2312428/68‑1 (MQ=255)
                                ccGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGATGgaa                          >  1:2971431/1‑68 (MQ=255)
                                ccGCGCGTGAAGCAAAAGTTCTGCGTATGCGGTCCGGTATCGATATGAACACCGACTACACGCGGgaa                          >  1:2678618/1‑68 (MQ=255)
                                 cgcgcgTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGgaag                         >  1:3167784/1‑68 (MQ=255)
                                 cgcgcgTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGgaag                         >  1:3930703/1‑68 (MQ=255)
                                  gcgcgTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGgaaga                        <  1:3076959/68‑1 (MQ=255)
                                     cgTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGt                     >  1:4065627/1‑68 (MQ=255)
                                      gagAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCTATATGAACACCGACTACACGCTGGAAGAAGTg                    <  1:619414/66‑1 (MQ=255)
                                      gTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCTATATGAACACCGACTACACGCTGGAAGAAGTg                    <  1:2014783/68‑1 (MQ=255)
                                      gTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTg                    <  1:1478878/68‑1 (MQ=255)
                                      gTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTg                    <  1:2268706/68‑1 (MQ=255)
                                      gTGAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTg                    <  1:391155/68‑1 (MQ=255)
                                        gAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTggg                  <  1:2906502/68‑1 (MQ=255)
                                        gAAGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTcgg                  <  1:1621991/68‑1 (MQ=255)
                                          aGCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTa                >  1:491269/1‑68 (MQ=255)
                                           gCAAAAGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTaa               >  1:507336/1‑68 (MQ=255)
                                              aaaGTTCTTCGTATGCGTTGCGGTATCGATATGATCAGCGACTAGACGCTGGAAGAAGTGTGTAAACa            <  1:2799622/68‑1 (MQ=255)
                                              aaaGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACa            <  1:2714311/68‑1 (MQ=255)
                                               aaGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAg           <  1:13958/68‑1 (MQ=255)
                                                aGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGt          <  1:3786335/68‑1 (MQ=255)
                                                aGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGt          >  1:3557505/1‑68 (MQ=255)
                                                aGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGt          >  1:3329576/1‑68 (MQ=255)
                                                aGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGt          <  1:2654892/68‑1 (MQ=255)
                                                aGTTCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGt          >  1:1561188/1‑68 (MQ=255)
                                                  ttCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTACACAGTTc        >  1:1933657/1‑68 (MQ=255)
                                                  ttCTGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTc        >  1:1939557/1‑68 (MQ=255)
                                                     tGCGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGAc     >  1:1259595/1‑68 (MQ=255)
                                                       cGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGt   >  1:1591307/1‑68 (MQ=255)
                                                       cGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGt   >  1:2155170/1‑68 (MQ=255)
                                                       cGTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGt   >  1:1810477/1‑68 (MQ=255)
                                                        gTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGtt  >  1:1340606/1‑68 (MQ=255)
                                                        gTATGCGTTCCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGtt  >  1:3591887/1‑68 (MQ=255)
                                                                 |                                                          
TGCGGCAACGCACGACGTGCTGGCTGGCCTGACCGCGCGTGAAGCAAAAGTTCTGCGTATGCGTTTCGGTATCGATATGAACACCGACTACACGCTGGAAGAAGTGGGTAAACAGTTCGACGTT  >  CP009273/3208028‑3208151

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: