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.95.0% 71.5 / ‑1.3 20coding (554/816 nt)BW25113_RS01680hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (10/9);  total (11/9)
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.29e-01

GTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATT  >  NZ_CP009273/340352‑340518
                                                                                        |                                                                              
gTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGaa                                                                               >  1:230372/1‑90 (MQ=255)
                   gTATTCTTCTTGGATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATgtc                                                           >  2:171674/1‑90 (MQ=255)
                   gTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATgt                                                            >  2:207180/1‑90 (MQ=255)
                     aTTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATgtcgt                                                         >  2:320818/1‑90 (MQ=255)
                                  aCCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagac                                            >  2:373366/1‑90 (MQ=255)
                                     atatGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacaga                                         >  1:25475/1‑90 (MQ=255)
                                         gCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgcg                                     <  1:91437/90‑1 (MQ=255)
                                             ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTa                                                      >  1:334127/1‑69 (MQ=255)
                                             ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTa                                                      <  2:334127/69‑1 (MQ=255)
                                             ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                 >  1:338122/1‑90 (MQ=255)
                                             ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                 >  1:161970/1‑90 (MQ=255)
                                                     aaTTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAg                         <  2:317893/90‑1 (MQ=255)
                                                          tCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCt                    <  2:25475/90‑1 (MQ=255)
                                                                aaGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCa              <  2:250070/90‑1 (MQ=255)
                                                                 aGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt             <  2:305356/90‑1 (MQ=255)
                                                                 aGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt             <  2:318330/90‑1 (MQ=255)
                                                                      aaGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt        <  2:100432/90‑1 (MQ=255)
                                                                      aaGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt        <  1:292054/90‑1 (MQ=255)
                                                                            gtATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattattatt  >  1:373842/1‑90 (MQ=255)
                                                                            gtATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattattatt  >  2:201074/1‑90 (MQ=255)
                                                                            gtATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattattatt  >  1:159580/1‑90 (MQ=255)
                                                                                        |                                                                              
GTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATT  >  NZ_CP009273/340352‑340518

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