Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,767,445 (C)7→9 intergenic (+90/‑265) rnlB → / → yfjP CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,767,4381.C83.3% 116.0 / 16.4 42intergenic (+83/‑272)rnlB/yfjPCP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base . (5/2);  new base C (14/21);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
*NC_0009132,767,4382.C82.9% 115.0 / 16.5 41intergenic (+83/‑272)rnlB/yfjPCP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base . (5/2);  new base C (13/21);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑CCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCG  >  NC_000913/2767305‑2767563
                                                                                                                                       ||                                                                                                                                       
tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAa‑‑ccac‑‑c                                                                                                                                >  2:538825/1‑139 (MQ=255)
 ctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc‑‑cc                                                                                                                               >  1:1182665/1‑139 (MQ=255)
 ctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑accc‑‑cc                                                                                                                               >  2:1457008/1‑139 (MQ=255)
    ttttGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc‑‑cc‑‑ccc                                                                                                                          <  2:1182665/139‑3 (MQ=255)
    ttttGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc‑‑cc‑‑ccc                                                                                                                          >  2:1303059/1‑137 (MQ=255)
    ttttGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc‑‑cc‑‑ccc                                                                                                                          <  2:1417033/139‑3 (MQ=255)
       tGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑CCCC‑‑CC‑‑C‑‑T‑‑caca                                                                                                                   <  2:1220730/139‑1 (MQ=255)
         tctcCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc‑‑cc‑‑ccctcaca                                                                                                                     <  1:1278360/139‑8 (MQ=255)
            cccGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCACGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACg                                                                                                                >  2:166299/1‑139 (MQ=255)
              cGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTc                                                                                                              >  1:1391263/1‑139 (MQ=255)
                     gTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATa                                                                                                       <  2:1329542/139‑1 (MQ=255)
                     gTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc‑‑c                                                                                                                                <  2:589466/118‑1 (MQ=255)
                     gTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc‑‑c                                                                                                                                >  1:589466/1‑118 (MQ=255)
                       aCGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGt                                                                                                               >  2:307983/1‑129 (MQ=255)
                       aCGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGt                                                                                                               <  1:307983/129‑1 (MQ=255)
                        cGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc                                                                                                                                   <  2:232935/114‑1 (MQ=255)
                        cGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc                                                                                                                                   >  1:232935/1‑114 (MQ=255)
                           cAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAac                                                                                                 <  1:1303059/139‑1 (MQ=255)
                           cAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAac                                                                                                 <  1:620359/139‑1 (MQ=255)
                              gTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc‑‑cc‑‑ccctcaca                                                                                                                     <  2:957870/118‑8 (MQ=255)
                              gTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc‑‑cc‑‑ccctcaca                                                                                                                     >  1:957870/1‑111 (MQ=255)
                              gTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc‑‑cc                                                                                                                               <  1:898449/110‑1 (MQ=255)
                              gTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc‑‑cc                                                                                                                               >  2:898449/1‑110 (MQ=255)
                                   gTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑CCACC‑‑CCAC‑‑T‑‑CAAATCGTCTTGCATAGGGTACACGTTTTc                                                                                         >  1:1431920/1‑139 (MQ=255)
                                         ggATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGt                                                                                                               <  1:679133/111‑1 (MQ=255)
                                         ggATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGt                                                                                                               >  2:679133/1‑111 (MQ=255)
                                         ggATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       >  2:189758/1‑135 (MQ=255)
                                         ggATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       <  1:189758/135‑1 (MQ=255)
                                                     ctctTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCACGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACt                                                                       <  1:166299/139‑1 (MQ=255)
                                                              ctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                              <  1:1188348/139‑1 (MQ=255)
                                                                        aggTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTaaa                                                    <  2:688549/137‑1 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAA‑‑CCCAC‑‑CC‑‑C‑‑T‑‑CACAACGTCACGCATAGGGAACACGATTTACCTTGTGCACCCACCACTAAATTTCCCCCCCTCAAAACGTCCCATTTATg                                       >  1:515637/1‑139 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CCCAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTtcc                                                              >  1:861358/1‑116 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTc                                                                                                              <  2:861049/68‑1 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTc                                                                                                              >  1:1194319/1‑68 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTc                                                                                                              <  2:1194319/68‑1 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTc                                                                                                              >  1:861049/1‑68 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGttt                                                                                           <  1:1561684/87‑1 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGttt                                                                                           >  2:1561684/1‑87 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                              >  1:312641/1‑100 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                              <  2:312641/100‑1 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                              <  2:861358/116‑1 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑TCCC‑‑CCAGT‑‑T‑‑CCCAACTCCCTGCCTAGGCAACAAGTGTTTCCCTGTGCTCCCAAGACAAAATTTCCCCCCCTTTCACCTTCTCAATAATg                                       >  2:184506/1‑139 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑CCCC‑‑CC‑‑C‑‑TTCCATAACGTACTCCTTAGTGAAAACGGTTTCCACTGTCCACCCACGCCTAAATTCCCCCCCCTTTCAACTCCTCcataata                                       >  1:859514/1‑138 (MQ=255)
                                                                                       ccTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc‑‑cc                                                                                                                               >  1:236859/1‑53 (MQ=255)
                                                                                       ccTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc‑‑cc                                                                                                                               <  2:236859/53‑1 (MQ=255)
                                                                                                 tGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTc                           <  2:859514/139‑1 (MQ=255)
                                                                                                     atTAACTTGTCCAATATGATCAAATAGCATTAA‑‑CCCC‑‑CC‑‑CCTT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGa                       <  1:184506/139‑1 (MQ=255)
                                                                                                       tAACTTGTCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGaaa                     <  2:1191681/139‑1 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                              <  1:1456403/75‑1 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                              >  2:1456403/1‑75 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGACGTGCATAGGGAACACGTTTTCCCCTCTGCACCCACGACTAGATTTCCCCCCCTTAAACTTCTCAATAATGTCACCAATTTTCTCGAAATACTatga             >  2:836839/1‑139 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑C‑‑CCCAACGTCCTGCATAGAGAAAACGTTTTCCACCGTCCACCCACGACTAAACATCCCACCCTTCTCCCTCCCCACTAATCTCACCAATTTCCTGGCAATCCTatt              >  1:621266/1‑138 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑A‑‑CACAACGTCCTCCATAGGGAAAAACTTTTACCCCGAGCACCCAAGACTAAAATTCACCCCCCTTAAACCTTTCAAAAATGACAACAATTTTTAGGCAAACCTata              >  1:1281577/1‑138 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAA‑‑CCCCCACC‑‑C‑‑T‑‑AACAACATCCAGCATAAGGAACACGATATTCTCTGGGCAACCACGACCAAATTTGACCCCCTTCGAAATTCGCAATAATATAACCAACTTTCGGGAAATACTata              >  1:997578/1‑138 (MQ=255)
                                                                                                                 actattatCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACCACCTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTCAACTTCTCACTAATGTCACCAATTTTCTGGAAATCCTatgatg           <  1:538825/137‑1 (MQ=255)
                                                                                                                   tatGATCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CAAAACGTCCTGCTTAGGGCACACGTTTTCCCCTGTGCACCCACGACTCAATTTCCCCCCCCTTAAACTTCTCAATAATGTCCCCCATTTTCTGGAAATCCAATGATGcc         >  1:109959/1‑139 (MQ=255)
                                                                                                                       aTCAAATAGCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc     <  2:621266/139‑1 (MQ=255)
                                                                                                                         cAAATAGCAGTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg  <  1:836839/139‑1 (MQ=255)
                                                                                                                          aaaTATCATTAACCCCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg  <  1:1457008/139‑1 (MQ=255)
                                                                                                                                       ||                                                                                                                                       
TCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑CCCC‑‑CC‑‑C‑‑T‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCG  >  NC_000913/2767305‑2767563

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