Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,184,613 A→T I207N (ATC→AAC)  phoQ ← two‑component system sensor histidine kinase PhoQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,184,6130AT95.2% 67.0 / ‑3.4 21I207N (ATC→AAC) phoQtwo‑component system sensor histidine kinase PhoQ
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (7/13);  total (7/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

ATGTTCTTCCAGTTCGCGGACTTCTTTTGCCAGGGCTTCGATGGGGCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGATCACTAACAGCAGATTGGCTGAGAGCACATAGATAAACCAGCTCCAGACCATATAGGAACTTTTTAGCTCCACCGGAATGGTATCC  >  NZ_CP009273/1184531‑1184699
                                                                                  |                                                                                      
aTGTTCTTCCAGTTCGCGGACTTCTTTTGCCAGGGCTTCGATGGGGCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGTTCACTaa                                                                                 >  1:431146/1‑90 (MQ=255)
  gTTCTTCCAGTTCGCGGACTTCTTTTGCCAGGGCTTCGATGGGGCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGTTCACTAAca                                                                               <  2:480146/90‑1 (MQ=255)
           gTTCGCGGACTTCTTTTGCCAGGGCTTCGATGGGGCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGTTCACTAACAGCAGATTgg                                                                      >  1:83851/1‑90 (MQ=255)
                 ggACTTCTTTTGCCAGGGCTTCGATGGGGCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGTTCACTAACAGCAGATTGGCTgaga                                                                <  1:689258/90‑1 (MQ=255)
                 ggACTTCTTTTGCCAGGGCTTCGATGGGGCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGTTCACTAACAGCAGATTGGCTgaga                                                                <  2:361994/90‑1 (MQ=255)
                  gACTTCTTTTGCCAGGGCTTCGATGGGGCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGTTCACTAACAGCAGATTGGCTGagag                                                               <  1:279230/90‑1 (MQ=255)
                      tcttTTGCCAGGGCTTCGATGGGGCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGTTCACTAACAGCAGATTGGCTGAGAGcaca                                                           <  1:267717/90‑1 (MQ=255)
                      tcttTTGCCAGGGCTTCGATGGGGCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGATCACTAACAGCAGATTGGCTGAGAGcaca                                                           <  1:702/90‑1 (MQ=255)
                              cAGGGCTTCGATGGGGCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGTTCACTAACAGCAGATTGGCTGAGAGCACATAGATaaa                                                   >  1:34560/1‑90 (MQ=255)
                                gggCTTCGATGGGGCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGTTCACTAACAGCAGATTGGCTGAGAGCACATAGATAAAcc                                                 <  2:34560/90‑1 (MQ=255)
                                      cGATGGGGCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGTTCACTAACAGCAGATTGGCTGAGAGcac                                                            <  2:209567/73‑1 (MQ=255)
                                      cGATGGGGCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGTTCACTAACAGCAGATTGGCTGAGAGcac                                                            >  1:209567/1‑73 (MQ=255)
                                        aTGGGGCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGTTCACTAACAGCAGATTGGCTGAGAGCACATAGATAAACCAGCTCCAg                                         <  2:575599/90‑1 (MQ=255)
                                           gggCGTAAACTCCACCAGGCGGCGACCCACAGCATCGGGTTCACTAACAGCAGATTGGCTGAGAGCACATAGATAAACCAGCTCCAGAcc                                      <  1:556249/90‑1 (MQ=255)
                                           gggCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGTTCACTAACAGCAGATTGGCTGAGAGCACATAGATAAACCAGCTCCAGAcc                                      >  1:434723/1‑90 (MQ=255)
                                           gggCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGTTCACTAACAGCAGATTGGCTGAGAGCACATAGATAAACCAGCTCCAGAcc                                      >  1:133563/1‑90 (MQ=255)
                                                           aGGCGGCGACCCACAGCAGCGGGTTCACTAACAGCAGATTGGCTGAGAGCACATAGATAAACCAGCTCCAGACCATATAGGAACTTTTTa                      >  2:698453/1‑90 (MQ=255)
                                                                gcgACCCACAGCAGCGGGTTCACTAACAGCAGATTGGCTGAGAGCACATAGATAAACCAGCTCCAGACCATATAGGAACTTTTTAGCTcc                 <  1:698453/90‑1 (MQ=255)
                                                                gcgACCCACAGCAGCGGGTTCACTAACAGCAGATTGGCTGAGAGCACATAGATAAACCAGCTCCAGACCATATAGGAACTTTTTAGCTcc                 <  1:278288/90‑1 (MQ=255)
                                                                        cagcagCGGGTTCACTAACAGCAGATTGGCTGAGAGCACATAGATAAACCAGCTCCAGACCATATAGGAACTTTTTAGCTCCACCGGAAt         <  2:133563/90‑1 (MQ=255)
                                                                               gggTTCACTAACAGCAGATTGGCTGAGAGCACATAGATAAACCAGCTCCAGACCATATAGGAACTTTTTAGCTCCACCGGAATGGTATcc  <  2:83851/90‑1 (MQ=255)
                                                                                  |                                                                                      
ATGTTCTTCCAGTTCGCGGACTTCTTTTGCCAGGGCTTCGATGGGGCGTAAACTCCACCAGGCGGCGACCCACAGCAGCGGGATCACTAACAGCAGATTGGCTGAGAGCACATAGATAAACCAGCTCCAGACCATATAGGAACTTTTTAGCTCCACCGGAATGGTATCC  >  NZ_CP009273/1184531‑1184699

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