Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,557,839 T→G intergenic (+108/‑43) BW25113_RS12770 → / → BW25113_RS12775 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,557,8390TG82.3% 85.0 / 14.2 34intergenic (+108/‑43)BW25113_RS12770/BW25113_RS12775hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base G (10/18);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGA  >  NZ_CP009273/2557705‑2557973
                                                                                                                                      |                                                                                                                                      
ggtagAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGtt                                                                                                                                     <  2:284494/138‑1 (MQ=255)
ggtagAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGtt                                                                                                                                     <  2:451756/138‑1 (MQ=255)
ggtagAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGtt                                                                                                                                     <  1:603666/138‑1 (MQ=255)
  tagAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAAc                                                                                                                                  <  2:71587/139‑1 (MQ=255)
          ttGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAAt                                                                                                                          >  1:587682/1‑139 (MQ=255)
            gAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATg                                                                                                                         >  1:243177/1‑138 (MQ=255)
                          attgCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATgg                                                                                                          <  1:651300/139‑1 (MQ=255)
                                 aTAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATaa                                                                                                   <  2:133668/139‑1 (MQ=255)
                                              aacaacCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTa                                                                                      >  1:357118/1‑139 (MQ=255)
                                                            ccGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAAttt                                                                        >  2:651490/1‑139 (MQ=255)
                                                            ccGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAAttt                                                                        <  2:587682/139‑1 (MQ=255)
                                                              ggCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCaa                                                                                                                           >  2:79719/1‑86 (MQ=255)
                                                              ggCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCaa                                                                                                                           <  1:79719/86‑1 (MQ=255)
                                                              ggCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGc                                                                      >  1:203055/1‑139 (MQ=255)
                                                              ggCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGc                                                                      >  2:413655/1‑139 (MQ=255)
                                                              ggCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGc                                                                      >  2:339951/1‑139 (MQ=255)
                                                                 cccTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGAc                                                                   <  2:393232/139‑1 (MQ=255)
                                                                   cTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACAt                                                                 >  2:721442/1‑139 (MQ=255)
                                                                      gCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCt                                                              >  2:534685/1‑139 (MQ=255)
                                                                                   tgATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGc                                                  <  1:326150/138‑1 (MQ=255)
                                                                                     aTCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGa                                                <  1:413655/138‑1 (MQ=255)
                                                                                      tCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATg                                              >  2:158997/1‑139 (MQ=255)
                                                                                          cGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATgg                                             <  2:357118/136‑1 (MQ=255)
                                                                                            tCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGa                                        <  1:721442/139‑1 (MQ=255)
                                                                                            tCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGa                                        <  1:339951/139‑1 (MQ=255)
                                                                                             cAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGATATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGa                                        <  1:328205/138‑1 (MQ=255)
                                                                                               atatTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGAc                                                                   >  1:750520/1‑109 (MQ=255)
                                                                                               atatTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGAc                                                                   <  2:750520/109‑1 (MQ=255)
                                                                                                      cTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAgg                              <  2:203055/139‑1 (MQ=255)
                                                                                                       ttCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGGTTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGa                                        >  2:568485/1‑128 (MQ=255)
                                                                                                       ttCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGa                                        <  1:568485/128‑1 (MQ=255)
                                                                                                                 cTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGt                   <  2:438176/139‑1 (MQ=255)
                                                                                                                          caccaATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCaaa          <  1:315894/139‑1 (MQ=255)
                                                                                                                                  cgcgGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagga  <  1:241268/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
GGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGA  >  NZ_CP009273/2557705‑2557973

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