Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 290,020 T→C E164G (GAA→GGA)  BW25113_RS01425 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273290,0200TC83.3% 67.1 / 9.4 30E164G (GAA→GGA) BW25113_RS01425hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (13/12);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

CCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTAACGAA  >  NZ_CP009273/289885‑290157
                                                                                                                                       |                                                                                                                                         
ccGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCGTTTTCCATATCGTAATTTcccc                                                                                                                                        <  2:278629/139‑1 (MQ=255)
                 gcccTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCtt                                                                                                                       <  1:236941/138‑1 (MQ=255)
                      cAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACg                                                                                                                  >  1:842488/1‑139 (MQ=255)
                             tgttCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACt                                                                                                           >  1:606450/1‑139 (MQ=255)
                                                atGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCa                                                                                        >  2:208888/1‑139 (MQ=255)
                                                 tGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCAt                                                                                       >  1:194658/1‑139 (MQ=255)
                                                      gAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAAtt                                                                                  >  2:565313/1‑139 (MQ=255)
                                                       aaTATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCttt                                                                                                                      <  2:487387/102‑1 (MQ=255)
                                                       aaTATCCTGAGACACTTTGATGCTTTACATCAAAATGCCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCttt                                                                                                                      >  1:487387/1‑102 (MQ=255)
                                                         tatCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATtcat                                                                               <  1:696289/139‑1 (MQ=255)
                                                                agaCACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGc                                                                        >  2:249283/1‑139 (MQ=255)
                                                                   cacTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGtt                                                                     <  1:823796/139‑1 (MQ=255)
                                                                      tttGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAg                                                                  >  1:254849/1‑139 (MQ=255)
                                                                                tACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGtt                                                                     <  2:83763/126‑1 (MQ=255)
                                                                                tACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGtt                                                                     >  1:83763/1‑126 (MQ=255)
                                                                                tACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTaa                                                        >  2:339355/1‑139 (MQ=255)
                                                                                     cAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGaa                                                                                    >  2:253680/1‑106 (MQ=255)
                                                                                     cAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGaa                                                                                    <  1:253680/106‑1 (MQ=255)
                                                                                         aTGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATg                                               <  1:120/139‑1 (MQ=255)
                                                                                                     gtaaatgtaaatGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTaa                                    <  1:14671/138‑1 (MQ=255)
                                                                                                     gtaaatgtaaatGGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTaa                                    <  2:606450/138‑1 (MQ=255)
                                                                                                      taaatgtaaatgGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGtt                                                                     <  2:769970/104‑1 (MQ=255)
                                                                                                      taaatgtaaatgGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGtt                                                                     >  1:769970/1‑104 (MQ=255)
                                                                                                        aatgtaaatgGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGa                                <  2:216760/139‑1 (MQ=255)
                                                                                                            taaatgGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCa                            >  2:111913/1‑139 (MQ=255)
                                                                                                               atgGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCatga                         >  1:578376/1‑139 (MQ=255)
                                                                                                                tgGCTTTTTCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCatgat                        <  1:339355/139‑1 (MQ=255)
                                                                                                                    tttttCCATATCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTc                    >  1:729443/1‑139 (MQ=255)
                                                                                                                      tttCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCg                   <  1:208888/138‑1 (MQ=255)
                                                                                                                          cATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTc              <  1:187894/139‑1 (MQ=255)
                                                                                                                            tatCGTAATTTCCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCa            <  2:254849/139‑1 (MQ=255)
                                                                                                                                       tCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTAACGaa  <  1:565313/138‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                         
CCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTAACGAA  >  NZ_CP009273/289885‑290157

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