Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,426,178 A→T 15.6% intergenic (+96/‑276) ydaY → / → tmpR pseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,426,1780AT15.6% 199.1 / 26.8 77intergenic (+96/‑276)ydaY/tmpRpseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base A (34/31);  new base T (7/5);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

ATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGA  >  NC_000913/1426079‑1426275
                                                                                                   |                                                                                                 
aTAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGa                                                                                                   >  2:1763759/1‑100 (MQ=255)
aTAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGa                                                                                                   >  2:1764009/1‑100 (MQ=255)
 tAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGAc                                                                                                  <  1:783427/100‑1 (MQ=255)
  aGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGAc                                                                                                  <  2:779412/99‑1 (MQ=255)
   gACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACag                                                                                                >  2:1088398/1‑100 (MQ=255)
   gACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACag                                                                                                <  1:1891185/100‑1 (MQ=255)
     cccGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACagag                                                                                              <  1:1850581/100‑1 (MQ=255)
       cGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGAt                                                                                            >  1:1805515/1‑100 (MQ=255)
        ggTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATa                                                                                           <  1:820326/100‑1 (MQ=255)
        ggTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATa                                                                                           <  1:2066890/100‑1 (MQ=255)
           gCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAAt                                                                                        <  1:2239207/100‑1 (MQ=255)
           gCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTAGAATGTCTGGGGTTCCTTTGCCTTTATCTCTGACAGAGATaaa                                                                                         >  2:1645907/1‑99 (MQ=255)
            cTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaagg                                                                                       <  1:26162/100‑3 (MQ=255)
             ttAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaagg                                                                                       <  1:1743464/99‑3 (MQ=255)
             ttAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCCGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaaggc                                                                                      <  2:2113269/100‑4 (MQ=255)
              tAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaaggc                                                                                      <  2:2368264/99‑4 (MQ=255)
               aaGGAGTTGTTGTACGCCTATTCGGTAATATGCCGTGCCCGAGGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAat                                                                                    >  1:1608785/1‑100 (MQ=255)
               aaGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaaggcaa                                                                                    <  1:734886/100‑6 (MQ=255)
                aGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAata                                                                                   <  1:1176801/100‑1 (MQ=255)
                 ggAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGATGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAatat                                                                                  >  2:104653/1‑100 (MQ=255)
                  gAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATAtt                                                                                 >  2:2338554/1‑100 (MQ=255)
                    gttgttGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTa                                                                               >  1:582217/1‑100 (MQ=255)
                     ttgttgTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAg                                                                              >  2:294701/1‑100 (MQ=255)
                       gttgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGcc                                                                            >  2:17059/1‑100 (MQ=255)
                         tgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCAc                                                                          >  2:280626/1‑100 (MQ=255)
                          gtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACt                                                                         <  1:283758/100‑1 (MQ=255)
                          gtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACt                                                                         >  2:1638937/1‑100 (MQ=255)
                               ccTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATcc                                                                    >  1:2015096/1‑100 (MQ=255)
                               ccTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAAAGAATATTTAGCCACTCATcc                                                                    >  2:562335/1‑100 (MQ=255)
                                cTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCg                                                                   >  1:765457/1‑100 (MQ=255)
                                 tATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCgg                                                                  <  1:2210759/100‑1 (MQ=255)
                                  aTTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGt                                                                 >  1:1396569/1‑100 (MQ=255)
                                   ttCGGTAATATCCCGTGTCTGACGTTATGCTGGATTGTCTGTGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTa                                                                <  2:1185106/100‑1 (MQ=255)
                                     nggTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTatt                                                              >  2:726593/2‑100 (MQ=255)
                                     cGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTatt                                                              >  1:423401/1‑100 (MQ=255)
                                     cGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTatt                                                              <  1:2332340/100‑1 (MQ=255)
                                       gTAATATCCCGTGCCCGACGTTATGCTGGAATGGCGGGGGTGCCTTTGCCGTTATCTCGGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattga                                                            <  2:1029204/100‑1 (MQ=255)
                                        tAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgat                                                           <  1:2259735/100‑1 (MQ=255)
                                         aaTATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgatt                                                          >  2:1815073/1‑100 (MQ=255)
                                          atatCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgattg                                                         <  1:1597993/100‑1 (MQ=255)
                                             tCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAgc                                                      >  1:705816/1‑100 (MQ=255)
                                                tgTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGcgcg                                                   >  1:363596/2‑100 (MQ=255)
                                                   gCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATg                                                <  1:424177/100‑1 (MQ=255)
                                                   gCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATg                                                >  1:436849/1‑100 (MQ=255)
                                                       cgatctATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAtt                                            <  2:1480039/95‑1 (MQ=255)
                                                            tATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAg                                       <  1:2079531/100‑1 (MQ=255)
                                                             aTGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAttagatcg                                      >  1:1480048/1‑94 (MQ=255)
                                                              tGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCa                                     >  1:91554/1‑100 (MQ=255)
                                                                cTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGt                                   >  2:3808/1‑100 (MQ=255)
                                                                  ggAATGGCTGGGGTGCCTTTGCCTTTTTCTCTGACAGGGATAAATGAATATTAAGCCACTCATCCGGTATTGATTGAGCGCGTTGAATTTGAAGcgg                                    >  2:1499066/1‑95 (MQ=255)
                                                                  ggAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGa                                 >  1:1260654/1‑100 (MQ=255)
                                                                    aaTGGCTGGGGTGCGTTTGCCTTTCTCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATc                               <  2:1074299/100‑1 (MQ=255)
                                                                    aaTGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATc                               <  2:2279552/100‑1 (MQ=255)
                                                                    aaTGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATc                               <  2:18456/100‑1 (MQ=255)
                                                                    aaTGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATc                               <  2:290684/100‑1 (MQ=255)
                                                                       ggCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTATCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCttt                            <  2:355734/100‑1 (MQ=255)
                                                                       ggCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCttt                            <  1:294703/100‑1 (MQ=255)
                                                                         cTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGc                          <  2:945271/100‑1 (MQ=255)
                                                                             ggTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGTGATGAATTTGAAGCAGTGATCTTTGCACTg                      >  2:2197231/1‑100 (MQ=255)
                                                                              gTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTgg                     >  1:434969/1‑100 (MQ=255)
                                                                               tgcctttgccttTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGa                    <  2:2124090/100‑1 (MQ=255)
                                                                                 cctttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATg                  >  1:1771975/1‑100 (MQ=255)
                                                                                   tttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAc                <  1:1088403/100‑1 (MQ=255)
                                                                                   tttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAc                <  1:310792/100‑1 (MQ=255)
                                                                                   tattcatttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAc                >  1:194578/7‑100 (MQ=255)
                                                                                    attcatttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAcc               >  1:614926/6‑100 (MQ=255)
                                                                                     ttcatttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCa              >  2:58093/5‑100 (MQ=255)
                                                                                     ttcatttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTAAGCGCGACGAATTTGAAGCAGTGATCTTTGCACTGGATGACCa              >  2:322065/5‑100 (MQ=255)
                                                                                      tcatttATCGCTGTCAGCGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGCTGACCAg             >  2:1617144/4‑100 (MQ=255)
                                                                                      gcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGTTCTTTGCACTGGATGACCAg             <  2:1452943/100‑1 (MQ=255)
                                                                                       catttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGt            >  2:1237851/3‑100 (MQ=255)
                                                                                        atttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTa           >  1:101416/2‑100 (MQ=255)
                                                                                          ttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTAtt         >  1:1209487/1‑100 (MQ=255)
                                                                                           tATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTAttt        >  1:442639/1‑100 (MQ=255)
                                                                                             tctcTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTc      >  2:267878/1‑100 (MQ=255)
                                                                                              ctctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCa     <  1:2337657/100‑1 (MQ=255)
                                                                                              ctctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCa     <  1:1826273/100‑1 (MQ=255)
                                                                                                ctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAgg   >  2:1693116/1‑100 (MQ=255)
                                                                                                 tGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGa  >  2:1171307/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
ATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGA  >  NC_000913/1426079‑1426275

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