Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,773,272 C→A intergenic (‑90/‑46) yfjV ← / → yfjW CP4‑57 prophage; putative arsenite transporter, N‑terminal fragment/CP4‑57 prophage; uncharacterized protein YfjW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,773,2720CA93.3% 42.8 / ‑3.1 15intergenic (‑90/‑46)yfjV/yfjWCP4‑57 prophage; putative arsenite transporter, N‑terminal fragment/CP4‑57 prophage; uncharacterized protein YfjW
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (6/8);  total (7/8)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TCGTTGATAGAGTCTTCAGCCATCAACACTGTATCGAGTAATGAAATTGGTTGCTAGAAAGTCGATGATTTGTTGCGTATGATTAAAAATCAACCTTAAAGGCTTTATACATGGATTTGCTATGGCTTTTGTATCAGAAAGAGAAATTGTA  >  NC_000913/2773197‑2773347
                                                                           |                                                                           
tCGTTGATAGAGTCTTCAGCCATCAACACTGTATCGAGTAATGAAATTGGTTGCTAGAAAGTCGATGATTTGTTGAGTAt                                                                         >  2:499540/1‑80 (MQ=255)
  gTTGATAGAGTCTTCAGCCATCAACACTGTATCGAGTAATGAAATTGGTTGCTAGAAAGTCGATGATTTGTTGAGTATGa                                                                       <  2:499522/80‑1 (MQ=255)
     gATAGAGTCTTCAGCCATCAACACTGTATCGAGTAATGAAATTGGTTGCTAGAAAGTCGATGATTTGTTGAGTATGATTa                                                                    <  1:499523/80‑1 (MQ=255)
                   ccATCAACACTGTATCGAGTAATGAAATTGGTTGCTAGAAAGTCGATGATTTGTTGAGTATGATTAAAAATCAACCTTaaa                                                     <  1:499524/81‑1 (MQ=255)
                        aaCACTGTATCGAGTAATGAAATTGGTTGCTAGAAAGTCGATGATTTGTTGAGTATGATTAAAAATCAACCTTAAAGGCtt                                                >  2:499537/1‑81 (MQ=255)
                                   gAGTAATGAAATTGGTTGCTAGAAAGTCGATGATTTGTTGAGTATGATTAAAAATCAACCTTAAAGGCTTTATACATGGAt                                     >  1:499538/1‑81 (MQ=255)
                                    aGTAATGAAATTGGTTGCTAGAAAGTCGATGGTTTGTTGAGTATGATTAAAAATCAACCTTAAAGGCTTTATACATGGAtt                                    <  2:499525/81‑1 (MQ=255)
                                       aaTGAAATTGGTTGCTAGAAAGTCGATGATTTGTTGCGTATGATTAAAAATCAACCTTAAAGGCTTTATACATGGATTTGc                                 >  2:499543/1‑81 (MQ=255)
                                                ggTTGCTAGAAAGTCGATGATTTGTTGAGTATGATTAAAAATCAACCTTAAAGGCTTTATACATGGATTTGCTATGGCttt                        >  1:499542/1‑81 (MQ=255)
                                                 gTTGCTAGAAAGTCGATGATTTGTTGAGTATGATTAAAAATCAACCTTAAAGGCTTTATACATGGATTTGCTATGGCtttt                       <  1:499526/81‑1 (MQ=255)
                                                   tGCTAGAAAGTCGATGATTTGTTGAGTATGATTAAAAATCAACCTTAAAGGCTTTATACATGGATTTGCTATGGCTTTTg                      >  2:499545/1‑80 (MQ=255)
                                                         aaaGTCGATGATTTGTTGAGTATGATTAAAAATCAACCTTAAAGGCTTTATACATGGATTTGCTATGGCTTTTGTATCAg                >  2:499544/1‑80 (MQ=255)
                                                            gTCGATGATTTGTTGAGTATGATTAAAAATCAACCTTAAAGGCTTTATACATGGATTTGCTATGGCTTTTGTATCAGAAag            <  2:499527/81‑1 (MQ=255)
                                                                   atTTGTTGAGTATGATTAAAAATCAACCTTAAAGGCTTTATACATGGATTTGCTATGGCTTTTGTATCAGAAAGAGAAAtt     <  1:499528/81‑1 (MQ=255)
                                                                       gttgAGTATGATTAAAAATCAACCTTAAAGGCTTTATACATGGATTTGCTATGGCTTTTGTATCAGAAAGAGAAATTGTa  <  2:499529/80‑1 (MQ=255)
                                                                           |                                                                           
TCGTTGATAGAGTCTTCAGCCATCAACACTGTATCGAGTAATGAAATTGGTTGCTAGAAAGTCGATGATTTGTTGCGTATGATTAAAAATCAACCTTAAAGGCTTTATACATGGATTTGCTATGGCTTTTGTATCAGAAAGAGAAATTGTA  >  NC_000913/2773197‑2773347

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