Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 284,666 C→A 11.5% E70* (GAG→TAG)  yagI ← CP4‑6 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273284,6660CA11.5% 79.2 / 2.5 35E70* (GAG→TAG) yagICP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (10/21);  new base A (2/2);  total (12/23)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

GATAGCCCTCCAGCAGGGCGTTCAGCTCGGCCTCGCCCAGCCAGGCAATCAACACCTTGCCGATGGCGGTGGCGTGCACCGGCAGGCGGCGGCCGATGCGTGAATAGGCGATGGCGGCCAGCTTGCCTTCAATCTTCTCGATATAGACCCCTTCACGCCCGTCCAGGATCCCCAGATGGGTGGTCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCA  >  CP009273/284396‑284935
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                             
gATAGCCCTCCAGCAGGGCGTTCAGCTCGGCCTCGCCCAGCCAGGCAATCAACACCTTGCCGATGGCGGTGGCGTGCACCGGCAGGCGGCGGCCGATGCGTGAATAGGCGATGGCGGCCAGCTTGCCTTCAATCTTCTCGATATAGACCCCTTCACGCACGGCCAGGATCCCCAGATGGGTGGTCTGCCCGGTCCGCCGGGATAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCACCGCTcaccc                                                                                                                                                                                                                                                                           >  1:101994/1‑271 (MQ=255)
    gCCCTCCAGCAGGGCGTTCAGCTCGGCCTCGCCCAGCCAGGCAATCAACACCTTGCCGATGGCGGTGGCGTGCACCGGCAGGCGGCGGCCGATGCGTGAATAGGCGATGGCGGCCAGCTTGCCTTCAATCTTCTCGATATAGACCCGTTCACGCCCGTNCAGGATCCCCAGATGGGTGGTCTGCCCGGTACGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCACGAACAGCt                                                                                                                                                                                                                                                                       >  1:747454/1‑275 (MQ=255)
         ccAGCAGGGCGTTCAGCTCGGCCTCGCCCAGCCAGGCAATCAACACCTTGCCGATGGCGGTGCCGTGCACCGGCAGGCGGCGGCCGATGCGTGAATAGGCGATGGCGGCCCGCTTGCCTTCAATCTTCTCGATATAGACCCCGTCACGCCCGTCCAGGATCCCCAGATGGATGGTCTGCCCGGNCCGCCGGGNCGGCTCCGTCAGCCAGCCCTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCCCGCGACCAGCTTCATg                                                                                                                                                                                                                                                                  >  1:560909/1‑275 (MQ=255)
                                            gcacgcgctcccTGGCCGGTGGCGGTCGCGGGCGCGGGCTGGCGGCGGGGGTTGCGTGAATGGGGGATGGCAGCCACCTTGCTTTCAATCTTCGCGATATAGACCCAGTCACGCCCGTCCAGGAATCCCAGATGGGTGGTCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTTCCTTATGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGGTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGAt                                                                                                                                                                                                                                <  2:747448/264‑1 (MQ=255)
                                                                                ggggtggggcgggcgATGCGTGAATAGGGGATGCGGGCCAGCTTGCCTTCTATCTTTTCGATATACAGCCCTGGTCGCCCGTCCAGTCACCCCGGGTGGGGGGTCTGCCCGGTCCGCCGGGACGGCTCCGTCATCCAGCCGGTTGCCTTCTGCGGAATATCGATGGAGCCCACGACAGAATGGCCGCGCTCGACCAGCTGCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGa                                                                                                                                                                                           <  2:171743/268‑1 (MQ=255)
                                                                                 gCAGGCGGCGGCCGATGCGTGAATAGGCGATGGCGGCCAGCTTGCCTTCAATCTTCTCGATATAGACCCCTTCACGCCCGTCCAGGATCCCCAGATGGGTGGTCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATTTAGCCGTGAAGCTGCCGGGCTTTTACCAACGAg                                                                                                                                                                                          >  1:171744/1‑275 (MQ=255)
                                                                                   tgggaccggcCTATGCTTTATTAGGCGTTGCCGGCCAGCTTGCCTTCATTCTTCTCTATATAGACCCCTTCACGCCCGTCCCGGATCCCCAGTTGGGTGGTCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGAAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTg                                                                                                                                                                                        <  1:351585/269‑1 (MQ=255)
                                                                                     gcggcggcCGTTGGTTGAATAGGCGATGGTGGCCAGCTTGCCGTCAATCTTCTCGATATAGACCCCGTCACGCCCGTCCCGGATCCCGAGATGGGTGGTCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTGCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGa                                                                                                                                                                                      <  1:274028/275‑1 (MQ=255)
                                                                                      cggcggcCGATGCGTGAATAGGCGATGGCGGCCAGCTTGCCTTCAATCTTCTCGATATAGACCCCTTCACGCCCGTCCAGGATCCCCAGATGGGTGGTCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGa                                                                                                                                                                                      <  1:204948/274‑1 (MQ=255)
                                                                                         gcgcgcggtcGTGCATGGGCGGTGGCGGCCACCTTGCCTTCCATCTTCTCGATGTAGACCCCTTCACGCACGTCCAGGATCCCCAGTTGGGTGGTCTGCCCGGTCAGCCGAGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAggg                                                                                                                                                                                   <  2:101994/265‑1 (MQ=255)
                                                                                            ccGATGCGTGAATAGGCGATGGCGGCCAGCTTGCCTTCAATCTTCTCGATATAGACCCCTTCACGCCCGTCCAGGATCCCCAGATGGGTGGTCTGCTCGGTCAGCCGGGACCGCTCCGTCACCCAGCCTTTTGCCTTCTGCCGAATATCGATGGCGCCCACGCCAAAAGGGCCGCGCGCGACCAGCTGCATGCCCAGGCGAAACTAGCCCATCACCGGGCGCGGAGCGATATAGCCGTGATGAGGCAGGGTTTTTTTGAGCGAGTGggggagata                                                                                                                                                                               >  2:204945/1‑270 (MQ=255)
                                                                                                    tGAATAGGCGATTGCGGCCAGCTTGCCTTCAATCTTCTCGCTATAGACCCCTTCACGCCCGTCCAGGATCCCCAGTTGGGTGGTCTGCCCGGTCCGACGCGACAGCTCCGGCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTc                                                                                                                                                                        <  1:428188/274‑1 (MQ=255)
                                                                                                    ggAATAGGCGATGGCGGCCAGCCTGCCTTCAATCTTCTCGATATAGACCCCTGCACGCCCGTCCACGATCCCCAGATGGGTGGTCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGAt                                                                                                                                                                                                                                <  1:778668/217‑1 (MQ=255)
                                                                                                                                                       ttCACGCCCGTCCAGGATCCCCAGATGGGTGGTCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATCTCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATCCCGGGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATCTAGCCGTGAAGCTCCAGGGCTTTTAGCAGCGATTGTAGGGTACTCTCGCTCAGCCCCAGCAGTTTGCTGAGTGCGGGGATCTTCGGCTCGGTGTCCtccgggt                                                                                                                    >  2:222051/1‑269 (MQ=255)
                                                                                                                                                       ttCACGCCCGTCCAGGATCCCCAGATGGATGGTCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTATCCGGGTTCTGATCGATATAGCCGTGAAGNGGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTAGGTGATCTTAAGGTCGGTGGCCTGCTCGt                                                                                                                    >  1:613569/1‑275 (MQ=255)
                                                                                                                                                           cGCCCGTACAGGCTCCCCAGATGGGTGGTCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGaa                                                                                                                <  1:222054/275‑1 (MQ=255)
                                                                                                                                                           cGCCCGGTCAGGATCCCCAGATGGGTGGTCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTGCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGaa                                                                                                                <  1:332130/275‑1 (MQ=255)
                                                                                                                                                              ccGTCCTGGGTCCCCAGATGTGTGTTCGGCCCGTTCCGCCGGGACAGCTCCGTCAGCCAGACTTTTCCCTTCGGCCGAAGAGCGATGGAGCCCACGACAAAATGGCCTCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACa                                                                                                              <  1:418983/274‑1 (MQ=255)
                                                                                                                                                                 tCCAGGATCCCCAGATGGGTGGTCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTc                                                                                                                                                                                                                                                                     >  2:489266/1‑120 (MQ=255)
                                                                                                                                                                     ggATCCCCAGATGGGTGGTCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGAc                                                                                                                                                                                                                                                                            <  1:275062/109‑1 (MQ=255)
                                                                                                                                                                     ggATCCCCAGATGGGTGGTCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGc                                                                                                                                                                                                                                                                 <  1:489269/120‑1 (MQ=255)
                                                                                                                                                                                  ggtggtCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGTGGGTACTCtt                                                                                                                                                                            >  2:375848/1‑192 (MQ=255)
                                                                                                                                                                                      gtCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTc                                                                                                                                                                        <  1:375851/192‑1 (MQ=255)
                                                                                                                                                                                        cTGCCCGGTCCGCCGGTACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCt                                                                                                                                                       <  1:319312/207‑1 (MQ=255)
                                                                                                                                                                                           cccGGTCAGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCGGAATATCGAAGAAGCCAACGACAAAATGGCCGCGATAGACCAGCATCATGcc                                                                                                                                                                                                                                                                >  2:657921/1‑99 (MQ=255)
                                                                                                                                                                                               gTCCTCCGTGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGTGCCCACGACAAAATGGCCGCGCTAGACCAGCTTCATGCCGAgg                                                                                                                                                                                                                                                            <  1:657927/99‑1 (MQ=255)
                                                                                                                                                                                                         aCAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGAc                                                                                                                                                                                                                                                                            <  2:273360/73‑1 (MQ=255)
                                                                                                                                                                                                         aCAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGAc                                                                                                                                                                                                                                                                            <  2:458690/73‑1 (MQ=255)
                                                                                                                                                                                                                  tCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATg                                                                                                                                                   >  1:116502/1‑185 (MQ=255)
                                                                                                                                                                                                                      ccagccTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCgg                                                                                                                                               <  2:116502/185‑1 (MQ=255)
                                                                                                                                                                                                                      ccaactttttGTNTTCTGCCGAATAGCGCGGGAGCCCACGACCATATGCGCGCGCTCGACTAGCGTCATGCCGAGGCGATACTTGCCGTTATCCGGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGc                                                       <  2:613563/267‑1 (MQ=255)
                                                                                                                                                                                                                          ccTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATAc                                                                                                                                                                                                                                                      >  1:387323/1‑78 (MQ=255)
                                                                                                                                                                                                                              ttGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGc                                                                                                                                                                                                                                                  <  2:387320/78‑1 (MQ=255)
                                                                                                                                                                                                                                                cGATGGAGCCCACGACAAAATGGCCGCGCTAGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGc                                                                                                                                                                                                               >  1:281686/1‑95 (MQ=255)
                                                                                                                                                                                                                                                    ggAGCCCACGACAAAATGGCCGCGCTAGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGCTGCa                                                                                                                                                                                                           <  2:281683/95‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gcgcTCGACCAGTTTCATGCCGAGGGGGTACTTGCCGTTCTCCCGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCANCAGTTTGCTGACGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCa  <  1:737203/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                             
GATAGCCCTCCAGCAGGGCGTTCAGCTCGGCCTCGCCCAGCCAGGCAATCAACACCTTGCCGATGGCGGTGGCGTGCACCGGCAGGCGGCGGCCGATGCGTGAATAGGCGATGGCGGCCAGCTTGCCTTCAATCTTCTCGATATAGACCCCTTCACGCCCGTCCAGGATCCCCAGATGGGTGGTCTGCCCGGTCCGCCGGGACAGCTCCGTCAGCCAGCCTTTTGCCTTCTGCCGAATATCGATGGAGCCCACGACAAAATGGCCGCGCTCGACCAGCTTCATGCCGAGGCGATACTTGCCGTTCTCCGGGTTCTGATCGATATAGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCA  >  CP009273/284396‑284935

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