Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,765,794 (C)7→8 intergenic (‑18/‑124) yfjM ← / → rnlA CP4‑57 prophage; protein YfjM/CP4‑57 prophage; RNase LS, toxin of the RnlAB toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,765,7871.C96.8% 101.0 / NA 32intergenic (‑11/‑131)yfjM/rnlACP4‑57 prophage; protein YfjM/CP4‑57 prophage; RNase LS, toxin of the RnlAB toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base C (7/23);  minor base A (0/1);  total (8/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATG  >  NC_000913/2765655‑2765906
                                                                                                                                     |                                                                                                                        
tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTTTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTg                                                                                                                           >  2:394067/1‑133 (MQ=255)
tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTg                                                                                                                           >  2:88441/1‑133 (MQ=255)
    cTCGTTCTTTTTCGATGTTGAAATTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcccc                                                                                                                     <  1:230891/133‑1 (MQ=255)
    cTCGTTCTGTTTCGACGTAGAATTTTTTTTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcccc                                                                                                                     <  2:265288/133‑1 (MQ=255)
    cTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcccc                                                                                                                     <  1:318283/133‑1 (MQ=255)
       gTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccccc                                                                                                                  <  2:343534/133‑1 (MQ=255)
       gTTCTGTTTCGACGTAGAATTTGTTCTTGTTTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccccc                                                                                                                  <  1:266160/133‑1 (MQ=255)
              ttCGACGTAGCATTTCTTCTTGTCTGCTCGCGAATAATTACGAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGa                                                                                                            <  2:167340/133‑1 (MQ=255)
                   cGTAGACTTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                                                       <  2:576618/133‑1 (MQ=255)
                   cGTAGAATTTGTTCTGGTCTGCTGACGAATAATTACTAAGTCTCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                                                       <  2:261904/133‑1 (MQ=255)
                      aGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTg                                                                                                    <  1:152555/133‑1 (MQ=255)
                       gAATTTTTTCTTGGCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGc                                                                                                   <  2:544977/133‑1 (MQ=255)
                         aTTTTTTCTTTTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGccc                                                                                                 <  1:394067/133‑1 (MQ=255)
                                  tGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCACCCCGGTCGACCATCCTGCCCTATAGTTgg                                                                                        >  1:245309/1‑133 (MQ=255)
                                        ctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                                                  <  2:562877/133‑1 (MQ=255)
                                        ctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                                                  <  2:324577/133‑1 (MQ=255)
                                        ctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                                                  <  2:141296/133‑1 (MQ=255)
                                         tgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccc                                                                                 <  1:559547/133‑1 (MQ=255)
                                                 aatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCa                                                                         <  1:558050/133‑1 (MQ=255)
                                                           tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTaaa                                                               >  1:590222/1‑133 (MQ=255)
                                                           tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTaaa                                                               >  1:569505/1‑133 (MQ=255)
                                                                aGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTATACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                                          <  1:117738/133‑1 (MQ=255)
                                                                     ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                     >  1:517651/1‑133 (MQ=255)
                                                                         ttGGTGAGATCAGGAAGAGGTTACCGAATCGAGGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAAt                                                 <  2:258733/133‑1 (MQ=255)
                                                                                 aTCAGGAAGATGTTCACGAAACGATGGATATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCa                                         <  2:262699/133‑1 (MQ=255)
                                                                                  tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                        <  1:88441/133‑1 (MQ=255)
                                                                                    aGGAAGAGGTTACCGAAACGATGGTCATCATCGCCCATTTTCTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                      <  2:245309/133‑1 (MQ=255)
                                                                                     ggAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                                     <  1:516188/133‑1 (MQ=255)
                                                                                           ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCTTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtatt                               >  1:443303/1‑133 (MQ=255)
                                                                                            gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                   <  1:227617/128‑1 (MQ=255)
                                                                                            gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGATGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                   >  2:227617/1‑128 (MQ=255)
                                                                                              tACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt                            >  2:603969/1‑133 (MQ=255)
                                                                                                 cGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTa                         <  2:626407/133‑1 (MQ=255)
                                                                                                     aCGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                     <  2:603409/133‑1 (MQ=255)
                                                                                                      cGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGtt                    <  2:569505/133‑1 (MQ=255)
                                                                                                              atcatcCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCata            <  2:590222/133‑1 (MQ=255)
                                                                                                               tcatcCCCCATTTTGTGGACTGACCCCCCCTGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCatat           <  1:242893/133‑1 (MQ=255)
                                                                                                                      ccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTa    >  2:467219/1‑133 (MQ=255)
                                                                                                                       cATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTAt   <  2:131066/133‑1 (MQ=255)
                                                                                                                        aTTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg  <  2:443303/133‑1 (MQ=255)
                                                                                                                                     |                                                                                                                        
TCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATG  >  NC_000913/2765655‑2765906

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