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.76.5% 25.6 / 15.5 17coding (554/816 nt)BW25113_RS01680hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base . (6/7);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATT  >  NZ_CP009273/340362‑340518
                                                                              |                                                                              
aTTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTcg                                                                    >  1:370702/1‑90 (MQ=255)
             tcttctTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCt                                                       >  1:273428/1‑90 (MQ=255)
                   tGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTcc                                                  >  1:77670/1‑90 (MQ=255)
                              tGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgc                                      <  1:511521/90‑1 (MQ=255)
                              tGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgc                                      <  1:90858/90‑1 (MQ=255)
                               gCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgcg                                     <  2:10156/90‑1 (MQ=255)
                                 cTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACgcgc                                    <  2:399515/90‑1 (MQ=255)
                                   ttGCTGAGAATTTTCGCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                 >  2:497871/1‑90 (MQ=255)
                                   ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                 >  1:246859/1‑90 (MQ=255)
                                   ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                 >  2:470760/1‑90 (MQ=255)
                                           aaTTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAg                         <  1:315981/90‑1 (MQ=255)
                                                       aGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt             <  2:389656/90‑1 (MQ=255)
                                                       aGGGAAAGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt             <  2:455298/90‑1 (MQ=255)
                                                        gggAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt             <  1:167891/90‑1 (MQ=255)
                                                            aaGAGTGTATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt        <  2:251796/90‑1 (MQ=255)
                                                             agagTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt        <  2:301688/90‑1 (MQ=255)
                                                                  gtATGCATGTGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattattatt  >  2:558884/1‑90 (MQ=255)
                                                                              |                                                                              
ATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATT  >  NZ_CP009273/340362‑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: