Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 272,712 (G)5→6 intergenic (+58/‑135) eyeA → / → ykfC CP4‑6 prophage; small RNA EyeA/CP4‑6 prophage; PF00078 domain‑containing protein YkfC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913272,7071.G95.0% 135.6 / 0.1 40intergenic (+53/‑140)eyeA/ykfCCP4‑6 prophage; small RNA EyeA/CP4‑6 prophage; PF00078 domain‑containing protein YkfC
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (18/20);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCAGCGTCTATCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACGAGGTGAA  >  NC_000913/272571‑272841
                                                                                                                                         |                                                                                                                                        
gcagcaGCGTCTATCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAggg                                                                                                                                      <  1:661983/140‑1 (MQ=255)
  agcagcGTCTATCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAggggg                                                                                                                                   >  1:44510/1‑140 (MQ=255)
      gcGTCTATCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAggggggcat                                                                                                                               >  2:643860/1‑136 (MQ=255)
      gcGTCTATCACGCTTGCGTTGCAGACTTCTGCCAAGTTGAGGCAAGCATGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAggggggcat                                                                                                                               <  2:262869/140‑5 (MQ=255)
                     gCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGAcgcg                                                                                                                 <  2:42283/140‑1 (MQ=255)
                         tGCTGACTTCTGCCAACTTGCGACAAGCAAGCATATAGACTGCGACGGACAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCgctgg                                                                                                             <  2:241418/140‑1 (MQ=255)
                           cTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCgctggct                                                                                                           >  2:428131/1‑140 (MQ=255)
                                ttCTGCCAACTTGCGGCAAGCAAGGATAAATAGTGTGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCgg                                                                                                      <  2:192985/140‑1 (MQ=255)
                                 tCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGa                                                                                                     <  1:109503/140‑1 (MQ=255)
                                   tGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGc                                                                                                   <  1:643860/140‑1 (MQ=255)
                                   tGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGc                                                                                                   >  2:79247/1‑140 (MQ=255)
                                   tGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCCGGCTGGCGGAGc                                                                                                   >  2:347461/1‑140 (MQ=255)
                                   tGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCCGGCTGGCGGAGc                                                                                                   <  1:347461/140‑1 (MQ=255)
                                     ccAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCtt                                                                                                 >  2:591001/1‑140 (MQ=255)
                                                           aaaGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGa                                                                           <  1:591001/140‑1 (MQ=255)
                                                                  tgcgagggCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCa                                                                    <  1:453894/135‑1 (MQ=255)
                                                                  gCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCa                                                                    >  1:591175/1‑140 (MQ=255)
                                                                  gCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAACGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCa                                                                    >  1:439367/1‑140 (MQ=255)
                                                                   cGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCAt                                                                   >  1:547179/1‑140 (MQ=255)
                                                                   cGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTAGCGAAGCGTCCGTAGTCGTCCTCGATGGGGAAAGCCCAt                                                                   >  2:197741/1‑140 (MQ=255)
                                                                   cGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACCCGCTGGCTCGCGAAACATCCGTAGTCGTCCGCGATGGGGAAAGCCCAt                                                                   >  2:323510/1‑140 (MQ=255)
                                                                     aCGGGCAGCCTCCTCTGTATGCCTGGGTCTAGGCAGGGAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTa                                                                 <  1:323510/140‑1 (MQ=255)
                                                                            gCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGACTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCg                                                          <  1:398627/140‑1 (MQ=255)
                                                                             cctcctCAGTATGCCCGAGTCCCGGTAGGGAAACCGGCGAAGGTCAGCGACGGAGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGaa                                                        <  1:201751/140‑1 (MQ=255)
                                                                                    aGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAg                                                  <  2:563535/140‑1 (MQ=255)
                                                                                    aGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAg                                                  <  2:165061/140‑1 (MQ=255)
                                                                                          ccTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGttt                                            >  1:106006/1‑140 (MQ=255)
                                                                                           cTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAg                                                  <  2:661623/133‑1 (MQ=255)
                                                                                           cTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAg                                                  >  1:661623/1‑133 (MQ=255)
                                                                                           cTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGAGAAAAGCCCATTACATGGTGAAGGGAGGCAGTTTg                                           >  2:319207/1‑140 (MQ=255)
                                                                                                tCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAAtgt                                      >  1:655361/1‑140 (MQ=255)
                                                                                                    ggcaggTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGttt                                  >  1:377632/1‑140 (MQ=255)
                                                                                                          tAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACg                            >  2:151670/1‑140 (MQ=255)
                                                                                                                  ggAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTaa                    <  2:148218/140‑1 (MQ=255)
                                                                                                                  ggAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTaa                    <  1:79247/140‑1 (MQ=255)
                                                                                                                   gAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTg                                 >  1:79648/1‑126 (MQ=255)
                                                                                                                   gAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTg                                 <  2:79648/126‑1 (MQ=255)
                                                                                                                   gAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAAc                   >  1:26585/1‑140 (MQ=255)
                                                                                                                   gAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCCCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCAAAGTGAAACAGTTTGATTGTGATTGCGACGTGAATTAAc                   >  2:286527/1‑140 (MQ=255)
                                                                                                                            cGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAAc          <  2:411424/140‑1 (MQ=255)
                                                                                                                            cGACGGATGTTAAGGTGGGCATTGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAAc          <  1:260076/140‑1 (MQ=255)
                                                                                                                             gACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACg         <  2:439367/140‑1 (MQ=255)
                                                                                                                             gACGGATGTTAAGGGAGGCATGGCTGCGACGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACg         <  1:319207/140‑1 (MQ=255)
                                                                                                                                    gTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACGAGGTGaa  <  2:547179/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
GCAGCAGCGTCTATCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACGAGGTGAA  >  NC_000913/272571‑272841

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