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, uncharacterized protein/CP4‑57 prophage, RNase LS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,765,7871.C98.2% 194.0 / ‑3.4 57intergenic (‑11/‑131)yfjM/rnlACP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, RNase LS
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (20/36);  total (21/36)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.31e-01

GGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAAT  >  NC_000913/2765650‑2765919
                                                                                                                                          |                                                                                                                                     
ggCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTg                                                                                                                                        >  2:271592/1‑138 (MQ=255)
ggCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTg                                                                                                                                        <  1:271592/138‑1 (MQ=255)
   tctctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTTTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcccc                                                                                                                                  <  1:831152/139‑1 (MQ=255)
    ctctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccc                                                                                                                                 <  1:297176/139‑1 (MQ=255)
    ctctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccc                                                                                                                                 <  2:426108/139‑1 (MQ=255)
     tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcccccc                                                                                                                                <  1:370249/139‑1 (MQ=255)
              tCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGAc                                                                                                                        >  2:629862/1‑139 (MQ=255)
                  tttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                                                                    <  1:194436/139‑1 (MQ=255)
                  tttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                                                                    <  1:715794/139‑1 (MQ=255)
                  tttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                                                                    <  1:868433/139‑1 (MQ=255)
                  tttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                                                                    <  2:865616/139‑1 (MQ=255)
                  tttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                                                                    <  1:477368/139‑1 (MQ=255)
                  tttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                                                                    <  2:276932/139‑1 (MQ=255)
                        cGTAGAATTTTCTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGATTGCCCCCCCCGGTAGACGATCCTGccc                                                                                                              <  1:682767/139‑1 (MQ=255)
                          tAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCta                                                                                                            >  2:228948/1‑139 (MQ=255)
                          tAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCta                                                                                                            >  1:876866/1‑139 (MQ=255)
                                   ttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                                                   <  2:819507/139‑1 (MQ=255)
                                   ttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                                                   <  1:361520/139‑1 (MQ=255)
                                   ttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                                                   <  2:1034378/139‑1 (MQ=255)
                                   ttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                                                   <  1:169046/139‑1 (MQ=255)
                                   ttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                                                   <  2:464814/139‑1 (MQ=255)
                                                          aCTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                                                                    <  1:106414/99‑1 (MQ=255)
                                                          aCTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                                                                    >  2:106414/1‑99 (MQ=255)
                                                              aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                       <  2:646074/139‑1 (MQ=255)
                                                              aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                       <  1:851589/139‑1 (MQ=255)
                                                              aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                       >  2:851589/1‑139 (MQ=255)
                                                              aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                       <  1:228948/139‑1 (MQ=255)
                                                              aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                       >  1:825746/1‑139 (MQ=255)
                                                              aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                       <  2:178972/139‑1 (MQ=255)
                                                              aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                       <  2:825746/139‑1 (MQ=255)
                                                                  tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTg                                                                                                                 <  1:833167/94‑1 (MQ=255)
                                                                  tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTg                                                                                                                 >  2:833167/1‑94 (MQ=255)
                                                                   cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                  <  2:423748/139‑1 (MQ=255)
                                                                   cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCTCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                  >  1:539603/1‑139 (MQ=255)
                                                                          ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                                                           >  2:676945/1‑139 (MQ=255)
                                                                          ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                                                           >  2:410264/1‑139 (MQ=255)
                                                                                  tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                   <  1:319082/139‑1 (MQ=255)
                                                                                  tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                   <  1:676945/139‑1 (MQ=255)
                                                                                     gaTCAGGAAGAGGTTTCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                <  2:833538/139‑1 (MQ=255)
                                                                                         aGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                             >  1:543260/1‑122 (MQ=255)
                                                                                         aGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                             <  2:543260/122‑1 (MQ=255)
                                                                                                ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                   >  1:493126/1‑125 (MQ=255)
                                                                                                ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                   <  2:493126/125‑1 (MQ=255)
                                                                                                ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTag                                     >  1:662168/1‑139 (MQ=255)
                                                                                                ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTag                                     >  2:135248/1‑139 (MQ=255)
                                                                                                   tACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                                  <  2:662168/139‑1 (MQ=255)
                                                                                                     ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                                                   <  2:200250/73‑1 (MQ=255)
                                                                                                     ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                                                   >  1:200250/1‑73 (MQ=255)
                                                                                                          aCGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCa                           <  1:171966/139‑1 (MQ=255)
                                                                                                           cGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                                                                    >  1:109802/1‑50 (MQ=255)
                                                                                                           cGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                                                                    <  2:109802/50‑1 (MQ=255)
                                                                                                               ggACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                             >  2:553032/1‑100 (MQ=255)
                                                                                                               ggACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                             <  1:553032/100‑1 (MQ=255)
                                                                                                               ggACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGt                      <  2:539603/139‑1 (MQ=255)
                                                                                                                      atcCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCATATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg               <  2:876866/139‑1 (MQ=255)
                                                                                                                       tcCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                                                   <  1:922191/55‑1 (MQ=255)
                                                                                                                       tcCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                                                   >  2:922191/1‑55 (MQ=255)
                                                                                                                           ccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGccc                                                                                                              >  2:696978/1‑40 (MQ=255)
                                                                                                                           ccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGccc                                                                                                              <  1:696978/40‑1 (MQ=255)
                                                                                                                           ccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATc          <  2:330483/139‑1 (MQ=255)
                                                                                                                                  gtgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                       <  1:935345/71‑1 (MQ=255)
                                                                                                                                  gtgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                       >  2:935345/1‑71 (MQ=255)
                                                                                                                                   tgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCCGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTCTTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAAt  >  2:170889/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                                     
GGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAAT  >  NC_000913/2765650‑2765919

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