Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,635,362 G→A H162H (CAC→CAT ydfN ← Qin prophage; putative side tail fiber assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,635,3620GA94.5% 83.3 / NA 55H162H (CAC→CATydfNQin prophage; putative side tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (19/33);  total (22/33)
Fisher's exact test for biased strand distribution p-value = 5.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCT  >  NC_000913/1635225‑1635387
                                                                                                                                         |                         
tGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGtatg                         >  1:5191870/1‑140 (MQ=17)
tGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGtatg                         <  2:719768/140‑1 (MQ=17)
 gcgcACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTTGTTTTCGTCCCCAAATGAGTACTGGTTGCGGGGGCGCTGGGGGTATGc                        >  2:3693171/1‑140 (MQ=17)
 gcgcACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGc                        >  2:4482450/1‑140 (MQ=17)
 gcgcACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGc                        >  1:245256/1‑140 (MQ=17)
    cACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGAt                     >  1:6652150/1‑140 (MQ=17)
    cACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGAt                     >  1:4483409/1‑140 (MQ=17)
    cACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGAt                     >  2:720508/1‑140 (MQ=17)
    cACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCAAAAGCTGTACTGGATGGGCTGGCGCTGTGGGTATGCGAt                     >  2:2624343/1‑140 (MQ=17)
     aCCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGAtt                    <  2:5770760/140‑1 (MQ=14)
       cGGCACTATTTGCGGTTCCGCTTATACTATGTGGATGCGGCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTa                  <  2:3489054/140‑1 (MQ=11)
       cGGCACTATTTGCGGTCCCGCTTATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTa                  <  1:2695309/140‑1 (MQ=14)
         gCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGCCCCCAAATGTGTCCTGCATGCGCTGCCGCTATGGGTATGCGATTTAAt                >  2:3081228/1‑140 (MQ=14)
          cACTATTTGCGGTCCCGCTAATTCTATGGGTCTGCGCCCAGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATg               <  2:3026228/140‑1 (MQ=14)
          cACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATg               <  2:6154715/140‑1 (MQ=14)
          cACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATg               <  1:5340056/140‑1 (MQ=14)
          cACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATg               <  2:6997225/140‑1 (MQ=14)
          agcTATTTGCGGTCCCGCTAATCCTATGGGTCTGCTCCCCGGTGTTATCCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATg               <  2:1927277/138‑1 (MQ=11)
           aCTATTTGCGGTCCCTCTAATACTAAGGGCATGCGCCTCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGc              <  1:2259121/140‑1 (MQ=11)
              aTTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCTAATCTGTAATGGATGCGCTGGCGCTGTGGGAATGCGATTTAATGCCGt           >  2:2259121/1‑140 (MQ=14)
               tttGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAGATCTGTACGGGAAGCGGTGGCGCTGAGGGTGTGCGATTTAATGCCGTc          >  2:2003362/1‑140 (MQ=14)
               tttGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTc          <  1:4233861/140‑1 (MQ=11)
                 tGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTCTTCGTGGATTTAGTGCCGTAATCAAACGACTATGTGGTTTTCGTCCCCAAACCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCt        <  2:4410867/140‑1 (MQ=11)
                 tGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCt        <  1:720508/140‑1 (MQ=11)
                 tGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCt        <  1:7460122/140‑1 (MQ=11)
                 tGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCt        >  2:7499931/1‑140 (MQ=14)
                 tGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCt        >  2:7460122/1‑140 (MQ=14)
                 tGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCt        <  1:5519829/140‑1 (MQ=11)
                 tGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCt        <  1:621257/140‑1 (MQ=11)
                 tGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCt        <  1:4038983/140‑1 (MQ=11)
                 tGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGAAGCGCTGGAGCTATGGGTATGCGATTTAATGCCGTCCt        >  1:2050211/1‑140 (MQ=14)
                 tGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCACCAAATCTGTACTGGATGGGCTTGCGCTGTGGGTATGCGATTTAATGCCGTCCt        >  2:2622529/1‑140 (MQ=11)
                 tGCGGTCCCGCTAATACGATGGGTATGCGCCCCGGCGTTTTTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCt        <  2:2687617/140‑1 (MQ=11)
                  gCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTg       <  2:4255012/140‑1 (MQ=11)
                   cGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGt      >  1:5730703/1‑140 (MQ=11)
                   cGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGt      <  1:254937/140‑1 (MQ=11)
                     gTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATg               <  2:4532783/129‑1 (MQ=11)
                     gTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATg               >  1:4532783/1‑129 (MQ=11)
                      tCCCGCTAATACTATTGGTATGCGCCCCGCTGTTATTAGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGTTcc   <  1:2371255/140‑1 (MQ=11)
                      tCCCGCTAATACTATGGGTATGCGCTCCGGTGTTATTAGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGTTcc   <  1:2481101/140‑1 (MQ=11)
                      tCCCGCTAATACTATGGGTATGCGCCCCGGTTTTACTCGAGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGTTcc   <  2:2588179/140‑1 (MQ=11)
                      tCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGTTcc   <  1:4600771/140‑1 (MQ=11)
                      tCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGTTcc   <  1:5259378/140‑1 (MQ=11)
                      tCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGTTcc   <  2:6652150/140‑1 (MQ=11)
                      tCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGTTcc   <  2:127307/140‑1 (MQ=11)
                      tCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGTTcc   <  2:529439/140‑1 (MQ=11)
                      tCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCATGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGTTcc   <  1:3955327/140‑1 (MQ=11)
                      tCCCGCTAATACTATCGGTATGCGCACCGATGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGTTcc   <  1:3176717/140‑1 (MQ=11)
                       cccGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCt  >  1:2588179/1‑140 (MQ=255)
                       cccGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCt        >  2:5420478/1‑134 (MQ=11)
                       cccGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCt        >  2:4339227/1‑134 (MQ=11)
                       cccGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCt        <  1:4339227/134‑1 (MQ=1)
                       cccGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCt        <  1:5420478/134‑1 (MQ=1)
                       cccGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGTTcc   >  2:1857786/1‑139 (MQ=11)
                       cccGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGTTcc   <  1:1857786/139‑1 (MQ=1)
                       cccGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGTTCCt  >  1:1548130/1‑140 (MQ=11)
                       cccGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTCGGTGTGTGATTTAATGCCGTCCTGTTCCt  >  1:2282461/1‑140 (MQ=14)
                       cccGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTAGCGCTGTGGGTATGAGATTTAGTGCCGTCCTGTTCCt  >  1:1927277/1‑140 (MQ=11)
                                                                                                                                         |                         
TGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCT  >  NC_000913/1635225‑1635387

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