Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,593,877 T→C A89A (GCA→GCG rpoH ← RNA polymerase sigma factor RpoH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,593,8770TC96.3% 177.0 / ‑0.7 54A89A (GCA→GCGrpoHRNA polymerase sigma factor RpoH
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (24/28);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

TTTGCGCTGCGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACAAACCGCAGGTGAGACAGGATCAGCGTTTTAGCTGCTTCCAGATCGC  >  NZ_CP009273/3593742‑3594010
                                                                                                                                       |                                                                                                                                     
tttGCGCTGCGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCt                                                                                                                                    >  2:215944/1‑139 (MQ=255)
        gcgctTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAg                                                                                                                            >  1:799090/1‑139 (MQ=255)
        gcgctTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAg                                                                                                                            >  2:288234/1‑139 (MQ=255)
          gctTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGc                                                                                                                          >  1:717502/1‑139 (MQ=255)
          gctTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGc                                                                                                                          <  1:516155/139‑1 (MQ=255)
              tggtggTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTtcat                                                                                                                               >  2:461609/1‑130 (MQ=255)
              tggtggTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTtcat                                                                                                                               <  1:461609/130‑1 (MQ=255)
                   gtCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTAc                                                                                                                 >  1:553018/1‑139 (MQ=255)
                   gtCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTAc                                                                                                                 >  1:279422/1‑139 (MQ=255)
                                 gATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCaaa                                                                                                   <  1:24341/139‑1 (MQ=255)
                                 gATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCaaa                                                                                                   <  1:414849/139‑1 (MQ=255)
                                           ttACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCtgtg                                                                                         >  1:248335/1‑139 (MQ=255)
                                              cGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGc                                                                                       <  1:215944/138‑1 (MQ=255)
                                              cGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGc                                                                                       <  1:449035/138‑1 (MQ=255)
                                              cGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGc                                                                                       <  1:830981/138‑1 (MQ=255)
                                                     cGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCa                                                                                 >  2:810580/1‑137 (MQ=255)
                                                     cGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCa                                                                                 >  2:764498/1‑137 (MQ=255)
                                                       tATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGc                                                                             <  2:266444/139‑1 (MQ=255)
                                                       tATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGc                                                                             <  2:464825/139‑1 (MQ=255)
                                                       tATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGc                                                                             <  1:810580/139‑1 (MQ=255)
                                                       tATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGc                                                                             <  1:194685/139‑1 (MQ=255)
                                                        aTTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGcc                                                                            >  1:608720/1‑139 (MQ=255)
                                                            gTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCa                                                                        <  2:717502/139‑1 (MQ=255)
                                                               gATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATaa                                                                     >  2:430057/1‑139 (MQ=255)
                                                                 tctcTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAgg                                                                                                                           <  1:143727/83‑1 (MQ=255)
                                                                 tctcTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAgg                                                                                                                           >  2:143727/1‑83 (MQ=255)
                                                                       cTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGa                                                               <  2:109382/137‑1 (MQ=255)
                                                                       cTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGa                                                               <  2:64960/137‑1 (MQ=255)
                                                                            aTCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAatat                                                        >  1:512589/1‑139 (MQ=255)
                                                                            aTCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGc                                                                                       >  2:866539/1‑108 (MQ=255)
                                                                            aTCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGc                                                                                       <  1:866539/108‑1 (MQ=255)
                                                                               cAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGaa                                                     <  1:724998/139‑1 (MQ=255)
                                                                               cAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGaa                                                     <  1:319556/139‑1 (MQ=255)
                                                                               cAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGaa                                                     >  2:307166/1‑139 (MQ=255)
                                                                               cAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGaa                                                     >  2:495481/1‑139 (MQ=255)
                                                                                 gTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGaac                                                    >  1:3902/1‑138 (MQ=255)
                                                                                    aaCGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGaacaac                                                 <  1:612303/138‑1 (MQ=255)
                                                                                      cGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACaaa                                              <  2:512589/139‑1 (MQ=255)
                                                                                      cGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAat                                                          <  1:60002/127‑1 (MQ=255)
                                                                                      cGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAat                                                          >  2:60002/1‑127 (MQ=255)
                                                                                        gCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACAAAcc                                            <  2:308710/139‑1 (MQ=255)
                                                                                         cGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACAAACCg                                           <  2:77426/139‑1 (MQ=255)
                                                                                              gagaCCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACAAACCGCAGGt                                      >  1:590033/1‑139 (MQ=255)
                                                                                              gagaCCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACAAACCGCAGGt                                      <  2:799090/139‑1 (MQ=255)
                                                                                                   cAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGc                                                                             <  2:432826/95‑1 (MQ=255)
                                                                                                   cAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGc                                                                             >  1:432826/1‑95 (MQ=255)
                                                                                                   cAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACAAACCGCAGGTgaga                                  >  1:2624/1‑138 (MQ=255)
                                                                                                      gcgcACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTAc                                                                 <  1:770963/104‑1 (MQ=255)
                                                                                                      gcgcACACCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTAc                                                                 >  2:770963/1‑104 (MQ=255)
                                                                                                          acacCCACTTCCGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACAAACCGCAGGTGAGACAGGAt                            >  1:255722/1‑137 (MQ=255)
                                                                                                                    ccGGGTTGAAACGGCGCACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACAAACCGCAGGTGAGACAGGATCAGCGTTTTAg                 <  1:288234/138‑1 (MQ=255)
                                                                                                                                  cgcACCGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACAAACCGCAGGTGAGACAGGATCAGCGTTTTAGCTGCTTCCAGATCGc  <  1:685907/139‑1 (MQ=255)
                                                                                                                                    cACCGCTTTCACCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACAAACCGCAGGTGAGACAGGATCAGCGTTTTAgct               >  2:616663/1‑124 (MQ=255)
                                                                                                                                    cACCGCTTTCACCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACAAACCGCAGGTGAGACAGGATCAGCGTTTTAgct               <  1:616663/124‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
TTTGCGCTGCGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACAAACCGCAGGTGAGACAGGATCAGCGTTTTAGCTGCTTCCAGATCGC  >  NZ_CP009273/3593742‑3594010

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