Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,775,085 T→C 10.4% intergenic (+64/‑460) yfjW → / → ypjI CP4‑57 prophage; uncharacterized protein YfjW/putative UPF0401 protein YpjI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,775,0850TC10.4% 78.3 / 3.9 39intergenic (+64/‑460)yfjW/ypjICP4‑57 prophage; uncharacterized protein YfjW/putative UPF0401 protein YpjI
Reads supporting (aligned to +/- strand):  ref base T (19/16);  new base C (2/2);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.57e-02

TAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGC  >  NC_000913/2774945‑2775227
                                                                                                                                            |                                                                                                                                              
tAGGTCTGCTAGCAATTTCTATAAAATTCAGTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGttt                                                                                                                                              <  2:519880/143‑1 (MQ=255)
 aGGTCTGCTAGTAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGtttt                                                                                                                                             <  2:41607/143‑1 (MQ=255)
   gTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTa                                                                                                                                           >  2:298600/1‑143 (MQ=255)
        cTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTACAATATCCTAAATTAAGCCTGTTTTTAAAGCa                                                                                                                                      >  1:383829/1‑143 (MQ=255)
        cTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAACTATCCTAAGTAAGACCTGTTTTTAACGCa                                                                                                                                      >  2:421104/1‑143 (MQ=255)
         tAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGCTTTAAAAACAg                                                                                                                                     <  1:456713/143‑1 (MQ=255)
             aaTTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGCTTTAAAAACAGGCtt                                                                                                                                 <  1:410985/143‑1 (MQ=255)
                ttCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGTGCCTTTTTATGAATATTAATTAAAATATACTAAGTTAAGACTGTATTTAAAGCAggatgcat                                                                                                                              >  1:107352/1‑137 (MQ=255)
                  cTATAAAATTCATTCATTCTTGCTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCg                                                                                                                            <  1:560048/143‑1 (MQ=255)
                     taAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGa                                                                                                                         >  1:559048/1‑143 (MQ=255)
                      aaaaTTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATCTATTAAAATATCCTAAGTCAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGaa                                                                                                                        <  2:293754/143‑1 (MQ=255)
                       aaaTTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAg                                                                                                                       <  2:569777/143‑1 (MQ=255)
                            acatcattcTTAAAGTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGACTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATACTACATTAAGCATATTATTAAAGCAGGCATGAGCGTGAGGcaggc                                                                                                                  >  2:238628/4‑139 (MQ=255)
                              ttcattcTTACGTCAAGAATGAGGTTAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGGTCCTTTTTAGGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGacac                                                                                                                <  2:64608/143‑1 (MQ=255)
                                    cTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCCAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATAc                                                                                                          <  1:263923/143‑1 (MQ=255)
                                         tttAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTACGTTAAGCCTGTTTTTAAACCAAACTTGAGCGTGAAGCGGACACTCATACCTACt                                                                                                     >  1:546345/1‑143 (MQ=255)
                                         tttAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGATTGAGCGTGAAGCGGACACTCATACCTgat                                                                                                     >  2:316127/1‑141 (MQ=255)
                                            aaGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCTAATTGGAGCCTTGTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTa                                                                                                  <  2:61545/143‑1 (MQ=255)
                                              gAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTAAAGCGGACACTACTACCTGCTTTaaa                                                                                                >  2:179796/1‑143 (MQ=255)
                                                atgatgTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGt                                                                                              <  1:486531/143‑1 (MQ=255)
                                                       gaatgaaAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCACGCTTGAGCCTGAAGCGCACACGCATCCCTCCCTCAAATTTGCGttc                                                                                       >  2:405928/1‑142 (MQ=255)
                                                           gaaAATTTATTCTTGTAACGCATTTGGCTCCAATTGTAGTCTTTTTATGACTATTTATTAAAATATCCTAAGTTAGGCCTGTTTTTAAAGCAGGCTTGAGCGTGATGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGt                                                                                   <  2:227691/143‑1 (MQ=255)
                                                            aaaaTTTATTCTTGTAACGCATTTGTCTCCAATTGGAGTCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTc                                                                                  <  1:98166/143‑1 (MQ=255)
                                                                tttATTCTTTTATCGCATTTGGCTTCAATTGGAGCATTATAATGACTATGTATGAAAATATCCTCAGTTAAGCCTGTTGTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTa                                                                              <  1:182514/143‑1 (MQ=255)
                                                                       ttGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATAc                                                                       >  2:109002/1‑143 (MQ=255)
                                                                        tGTATCGCATTTGGCTCCAATTTGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAACGCAGGCTTGAGCGTGAAGCGGACACTCATCCCTGCTTTAAAGTGGCG‑TTTTGTCGTTACCCATACCt                                                                     <  2:201449/143‑1 (MQ=255)
                                                                           aaCGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTACGTTAAGCCCGTTTTTAAAGCAGGCTTGAGCGTGAAACGCACACTCATACCTACTTCAAAAAAGAGATTATGTAGTTACCCATACCtat                                                                   >  1:521394/1‑143 (MQ=255)
                                                                           aaCGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCGAAGTTAATCCTGTTTTTAAAGCTGTCTTGAGCGTGAAGCAGACACTCATACCTGATTTAAAATGGCGATTTTGTCGTTACCCATATCtat                                                                   >  1:521179/1‑143 (MQ=255)
                                                                                 tttGGCTCCAACTGGAGGCTTTTTATGACTATTTATTAAAATATCCTAGGTTAAGCCTGTTTCTAAAGCAGCAGTTAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTc                                                             <  2:35074/143‑1 (MQ=255)
                                                                                       tCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCtgttat                                                       >  1:113858/1‑141 (MQ=255)
                                                                                        cctaTTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTAGTtttgtt                                                      <  2:93765/140‑1 (MQ=255)
                                                                                          aaTTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTACTTTAAAGTGGCGTTATTGTCGTTACCCATACCTATATTTTCTGTTGTTTa                                                    >  1:529922/1‑143 (MQ=255)
                                                                                            ttGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAGGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTCCCCATACCTATATTTTCTGTCGTTTAAc                                                  >  2:17840/1‑143 (MQ=255)
                                                                                                                 ttattaAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTAAATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGaa                             >  2:433514/1‑143 (MQ=255)
                                                                                                                        aaTATGCTCAGTTATGCCTGTTTTTAAATCACGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAg                      <  2:559048/143‑1 (MQ=255)
                                                                                                                        aaTATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAg                      >  1:301241/1‑143 (MQ=255)
                                                                                                                         atatCCTAGGTTAAGCTTGTTGTTAGTGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGa                     <  2:215881/143‑1 (MQ=255)
                                                                                                                                aaGTTAAGCCTGTTTTTAAAGCATGCTTGAGCGTGAAGCGGACACTCTTACCTTCTTTAAAGTGGCGTATTTGTCTTTACGCATACCTATATTTTCTGTTGTTTAACTGCCATAACTGTTCGGTTGAAAAGCCAGAAATATGc              <  2:465709/143‑1 (MQ=255)
                                                                                                                                   tcaaGCCTGCTTTAAAAACAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTACAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTTACTGCCAGAACTGAGCGGTTGAAAAGCAAGAAATATGCCGa           >  1:298362/3‑143 (MQ=255)
                                                                                                                                   tcaaGCCTGCTTTAAAAACAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTATAACTGACAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGa           >  1:492784/3‑143 (MQ=255)
                                                                                                                                          tGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTCCCCCTCCCTATATTTTCTGTTGTTTAACTGTCAGAACTGAGCAGTTAAAAAGCCAGAAATATGCCGAGCCAGCg    >  1:192014/1‑143 (MQ=255)
                                                                                                                                            tttttAAAGAAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACTCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGc  <  1:455120/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                              
TAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGC  >  NC_000913/2774945‑2775227

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑