Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,765,794 (C)7→6 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,7880C.98.0% 198.1 / ‑2.9 51intergenic (‑12/‑130)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 C (0/0);  major base . (16/33);  minor base T (0/1);  total (16/35)
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

GGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATG  >  NC_000913/2765650‑2765920
                                                                                                                                          |                                                                                                                                    
ggCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcc                                                                                                                                     <  2:138596/140‑1 (MQ=255)
   tctctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccc                                                                                                                                  <  2:15956/140‑1 (MQ=255)
    ctctGCTCCTTCTGTTTCGACGCTGAATTTTATTTCGTCAGCTGGCGAGTAATTACTAAGTATCGAGTTGCCGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcctccc                                                                                                                                 <  2:392217/140‑1 (MQ=255)
     tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccccg                                                                                                                                >  1:786483/1‑139 (MQ=255)
     tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccccg                                                                                                                                >  1:62705/1‑139 (MQ=255)
     tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccccg                                                                                                                                >  1:452995/1‑139 (MQ=255)
     tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcc                                                                                                                                     >  1:21166/1‑135 (MQ=255)
     tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcc                                                                                                                                     <  2:21166/135‑1 (MQ=255)
               cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATc                                                                                                                     <  1:534743/140‑1 (MQ=255)
               cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATc                                                                                                                     <  2:1926/140‑1 (MQ=255)
               cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATc                                                                                                                     <  2:452995/140‑1 (MQ=255)
               cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATc                                                                                                                     >  2:752283/1‑140 (MQ=255)
               cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATc                                                                                                                     <  2:778616/140‑1 (MQ=255)
               cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATc                                                                                                                     >  2:798354/1‑140 (MQ=255)
                  tttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTg                                                                                                                  <  1:752283/140‑1 (MQ=255)
                     cGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGccc                                                                                                               <  2:513877/140‑1 (MQ=255)
                       aCGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCta                                                                                                             >  2:35738/1‑140 (MQ=255)
                       aCGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCta                                                                                                             <  2:32737/140‑1 (MQ=255)
                            gAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGtt                                                                                                        >  2:42575/1‑140 (MQ=255)
                                 ttttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGt                                                                                                   <  1:606913/140‑1 (MQ=255)
                                 ttttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGt                                                                                                   <  2:62705/140‑1 (MQ=255)
                                    tCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                                                                <  1:817103/140‑1 (MQ=255)
                                    tCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACAATCCTGCCCTATAGTTGGAGTGcc                                                                                                <  1:76304/140‑1 (MQ=255)
                                     cTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccc                                                                                               <  1:679795/140‑1 (MQ=255)
                                       tGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccccc                                                                                             >  2:833194/1‑140 (MQ=255)
                                       tGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccccc                                                                                             >  2:173938/1‑140 (MQ=255)
                                                    ataatTACTAAGTATAGAGTTGCTGCTTGCTGAGATCTGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCt                                                                                <  1:256463/140‑1 (MQ=255)
                                                          aCTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAAttt                                                                          >  1:599277/1‑140 (MQ=255)
                                                            tAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                                                        <  1:724813/140‑1 (MQ=255)
                                                            tAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                                                        <  2:599277/140‑1 (MQ=255)
                                                            tAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGCGGCTAAATTTAg                                                                        <  2:37308/140‑1 (MQ=255)
                                                               gTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCaa                                                                     <  1:173938/140‑1 (MQ=255)
                                                                 aTCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                   <  1:798354/140‑1 (MQ=255)
                                                                    gAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCaa                                                                <  1:833194/140‑1 (MQ=255)
                                                                     aGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                                                        <  1:685052/131‑1 (MQ=255)
                                                                     aGTTGCTGCTTGGAGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                                                        >  2:685052/1‑131 (MQ=255)
                                                                      gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                              <  1:167391/140‑1 (MQ=255)
                                                                      gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                              <  2:807110/140‑1 (MQ=255)
                                                                      gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                              >  2:72348/1‑140 (MQ=255)
                                                                      gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                              <  2:17250/140‑1 (MQ=255)
                                                                       ttGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                                                             >  1:191623/1‑140 (MQ=255)
                                                                        tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                                                            <  1:35738/140‑1 (MQ=255)
                                                                          ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGGTGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                                                          >  1:662575/1‑140 (MQ=255)
                                                                          ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTAGCTAAATTTTGCAAATCAATACACt                                                          >  1:299056/1‑140 (MQ=255)
                                                                                   gagaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                 <  2:786483/140‑1 (MQ=255)
                                                                                     gaTCAGGAAGAGGTTGCCGAAACGAGGGACATCATCCTCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                               <  1:378534/140‑1 (MQ=255)
                                                                                     gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCATGTGCCTATATTTAGCAGATCAGTACACTTCAGGGGGGta                                               >  2:259960/1‑140 (MQ=255)
                                                                                            aaGAGGTTACCGAAACGGTGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGt                                        <  2:662575/140‑1 (MQ=255)
                                                                                                 gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                                   >  2:481240/1‑140 (MQ=255)
                                                                                                 gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                                   <  1:626853/140‑1 (MQ=255)
                                                                                                 gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                                   <  1:72348/140‑1 (MQ=255)
                                                                                                 gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGATTATTGTAGAGt                                   >  2:626853/1‑140 (MQ=255)
                                                                                                 gTTACCGAAACGATGGACATCATCCCCCATTCTGTAGACGGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                                   <  2:299056/140‑1 (MQ=255)
                                                                                                  ttACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAGATCAATACACTTAAGGGGTGTATTATTGTAGAGtt                                  >  2:395501/1‑140 (MQ=255)
                                                                                                      cGAAACGATGGACATCAACCCCTATTTTGTGAAGTTCCCACCGGTAGATGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAAGGGGGTATTATTGTAGAGTTTccc                              <  2:393925/140‑1 (MQ=255)
                                                                                                          aCGATGGACATCATCACCCATTATGAGGACTGCCCCCTGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCata                          <  1:259960/140‑1 (MQ=255)
                                                                                                             aTGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGt                       <  1:481240/140‑1 (MQ=255)
                                                                                                                    tcatcCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg                <  1:42575/140‑1 (MQ=255)
                                                                                                                                  gtgGACTGCCCCCCGCTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATg  <  1:395501/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                    
GGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATG  >  NC_000913/2765650‑2765920

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