Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,433,698 T→C intergenic (‑24/+293) pinR ← / ← ynaE Rac prophage; putative site‑specific recombinase/Rac prophage; uncharacterized protein YnaE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,433,6980TC98.3% 128.7 / ‑4.3 59intergenic (‑24/+293)pinR/ynaERac prophage; putative site‑specific recombinase/Rac prophage; uncharacterized protein YnaE
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (20/38);  total (20/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

TTCGATTATTTGCTGAGGTTTGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAATGCGTATGCAACCCCGTAATTT  >  NC_000913/1433561‑1433801
                                                                                                                                         |                                                                                                       
ttcGATTATTTGCTGAGGTTTGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGgcgt                                                                                                       >  2:91947/1‑140 (MQ=17)
 tcGATTATTTGCTGAGGTTTGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGtt                                                                                                      <  2:534742/140‑1 (MQ=17)
   gATTATTTGCTGAGGTTTGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTtca                                                                                                    <  2:489089/140‑1 (MQ=18)
   gATTATTTGCTGAGGTTTGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTtca                                                                                                    <  1:1163296/140‑1 (MQ=18)
   gATTATTTGCTGAGGTTTGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTtca                                                                                                    <  1:96120/140‑1 (MQ=17)
   gATTATTTGCTGAGGTTTGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTtca                                                                                                    <  2:130071/140‑1 (MQ=18)
     ttattTGCTGAGGTTTGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTtcata                                                                                                  >  1:213223/1‑140 (MQ=17)
        tttGCTGAGGTTTGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTtcatatca                                                                                               <  1:1154854/140‑1 (MQ=18)
         ttGCTGAGGTTTGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTtcatatcat                                                                                              >  1:526574/1‑140 (MQ=17)
          tGCTGAGGTTTGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTtcatatcata                                                                                             <  1:18260/140‑1 (MQ=18)
              gAGGTTTGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGc                                                                                         >  1:1134070/1‑140 (MQ=17)
               aGGTTTGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCa                                                                                        >  2:1093162/1‑140 (MQ=17)
               aGGTTTGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCa                                                                                        >  2:222875/1‑140 (MQ=17)
               aGGTTTGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCa                                                                                        >  2:714754/1‑140 (MQ=17)
                  tttGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCAtct                                                                                     <  1:455524/140‑1 (MQ=17)
                   ttGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCAtctc                                                                                    >  1:1018635/1‑140 (MQ=17)
                    tGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCa                                                                                        >  2:892475/1‑135 (MQ=17)
                    tGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCa                                                                                        <  1:892475/135‑1 (MQ=18)
                    tGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCa                                                                                   <  2:1005277/140‑1 (MQ=18)
                       ttttAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGaa                                                                                <  1:410778/140‑1 (MQ=18)
                           aaaaCCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAaata                                                                            >  1:458712/1‑140 (MQ=17)
                           aaaaCCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAaata                                                                            >  2:320342/1‑140 (MQ=17)
                                cTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTAt                                                                       <  2:1114883/140‑1 (MQ=18)
                                cTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTAt                                                                       <  2:526574/140‑1 (MQ=18)
                                 tGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATg                                                                      >  2:874850/1‑140 (MQ=17)
                                   cACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGat                                                                    <  1:714754/140‑1 (MQ=17)
                                   cACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGat                                                                    <  1:1093162/140‑1 (MQ=17)
                                            tttCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAATAAGCAAAAATTCGAAACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCAACTCAGAAAATAATTATGATTATTTTtgt                                                           <  2:480353/140‑1 (MQ=255)
                                             ttCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTtgtg                                                          <  1:874850/140‑1 (MQ=17)
                                             ttCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTtgtg                                                          <  2:1001615/140‑1 (MQ=17)
                                             ttCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTtgtg                                                          <  2:846265/140‑1 (MQ=17)
                                              tCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGc                                                         <  2:566685/140‑1 (MQ=17)
                                               cccGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCa                                                        <  2:302277/140‑1 (MQ=17)
                                                ccGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCat                                                       <  1:91947/140‑1 (MQ=17)
                                                 cGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCata                                                      <  2:985667/140‑1 (MQ=17)
                                                 cGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCata                                                      <  1:681238/140‑1 (MQ=17)
                                                 cGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCata                                                      <  1:324775/140‑1 (MQ=17)
                                                  ggCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCatat                                                     <  1:325729/140‑1 (MQ=17)
                                                  ggCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCatat                                                     <  2:1134070/140‑1 (MQ=17)
                                                  ggCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCatat                                                     >  2:325729/1‑140 (MQ=17)
                                                   gCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATAtt                                                    <  2:1019152/140‑1 (MQ=17)
                                                    cGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATAttt                                                   <  2:156011/140‑1 (MQ=32)
                                                                 tggtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGtt                                      <  2:1018635/140‑1 (MQ=17)
                                                                        aTCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAAt                               >  1:864466/1‑140 (MQ=18)
                                                                             gCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACg                          <  2:1163225/140‑1 (MQ=17)
                                                                              cGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGa                         >  1:136216/1‑140 (MQ=17)
                                                                              cGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGa                         >  2:190214/1‑140 (MQ=18)
                                                                                ttGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAAt                       <  2:265204/140‑1 (MQ=17)
                                                                                ttGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAAt                       <  2:213223/140‑1 (MQ=17)
                                                                                  gATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAATGc                     <  1:190214/140‑1 (MQ=17)
                                                                                  gATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAATGc                     >  1:991854/1‑140 (MQ=18)
                                                                                  gATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAATGc                     >  2:1028546/1‑140 (MQ=18)
                                                                                  gATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAATGc                     >  1:47709/1‑140 (MQ=18)
                                                                                   atatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAAt                       >  1:1005307/1‑137 (MQ=18)
                                                                                   atatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAAt                       <  2:1005307/137‑1 (MQ=17)
                                                                                         gACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTATATGTACACGTTCGAAAATAAACGAATGCGTATGCa              <  2:458712/140‑1 (MQ=17)
                                                                                          aCAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAATGCGTATGCaa             <  2:45571/140‑1 (MQ=17)
                                                                                               aaGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAATGCGTATGCAACCCCg        <  1:913491/140‑1 (MQ=17)
                                                                                               aaGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAATGCGTATGCAACCCCg        <  2:261443/140‑1 (MQ=17)
                                                                                                     aaaTTCGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAATGCGTATGCAACCCCGTAAttt  >  1:648180/1‑140 (MQ=17)
                                                                                                                                         |                                                                                                       
TTCGATTATTTGCTGAGGTTTGATTTTAAAACCTGCACTTTCGATTTCCCGGCGTTGATTTTCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTCGAAAATAAACGAATGCGTATGCAACCCCGTAATTT  >  NC_000913/1433561‑1433801

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