Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 340,444 (A)5→4 coding (558/816 nt) BW25113_RS01680 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273340,4400A.96.5% 110.7 / ‑1.8 29coding (554/816 nt)BW25113_RS01680hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (18/10);  total (19/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

GTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTA  >  NZ_CP009273/340352‑340525
                                                                                        |                                                                                     
gTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGaa                                                                                      >  1:302373/1‑90 (MQ=255)
          aTTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTcg                                                                           >  1:238550/1‑90 (MQ=255)
          aTTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTcg                                                                           >  1:405686/1‑90 (MQ=255)
          aTTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAGAGTGTATGCATGTGGAAAATGGCTTcgc                                                                          >  2:107704/1‑90 (MQ=255)
                      ttcttcTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATgtcgtc                                                               <  1:226576/90‑1 (MQ=255)
                          tcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTAtt                                                           >  1:191748/1‑90 (MQ=255)
                          tcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTAtt                                                           >  1:2420/1‑90 (MQ=255)
                                     atatGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacaga                                                >  2:93097/1‑90 (MQ=255)
                                        tGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgc                                             >  2:324885/1‑90 (MQ=255)
                                         cccTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgcg                                            <  2:9086/89‑1 (MQ=255)
                                             ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                        >  1:387880/1‑90 (MQ=255)
                                             ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                        >  1:332846/1‑90 (MQ=255)
                                             ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                        >  2:360756/1‑90 (MQ=255)
                                             ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                        >  1:504416/1‑90 (MQ=255)
                                             ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                        >  1:53083/1‑90 (MQ=255)
                                                                aaGGGAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                    <  1:107704/90‑1 (MQ=255)
                                                                 aGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                    <  2:238550/90‑1 (MQ=255)
                                                                 aGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                    <  2:423303/90‑1 (MQ=255)
                                                                      aaGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt               <  1:84364/90‑1 (MQ=255)
                                                                      aaGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt               <  2:332846/90‑1 (MQ=255)
                                                                      aaGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt               <  2:504416/90‑1 (MQ=255)
                                                                          gtgtATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCtt                          <  1:487933/75‑1 (MQ=255)
                                                                          gtgtATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCtt                          >  2:487933/1‑75 (MQ=255)
                                                                          gtgtATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattatta           >  2:289012/1‑90 (MQ=255)
                                                                            gtATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgc                                             <  1:70858/54‑1 (MQ=255)
                                                                            gtATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgc                                             >  2:70858/1‑54 (MQ=255)
                                                                                  aTGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACtt   >  1:252244/1‑90 (MQ=255)
                                                                                  aTGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACtt   >  1:467710/1‑90 (MQ=255)
                                                                                  aTGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACtt   >  1:42211/1‑90 (MQ=255)
                                                                                    gtgGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTa  >  1:295792/1‑90 (MQ=255)
                                                                                        |                                                                                     
GTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTA  >  NZ_CP009273/340352‑340525

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