Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,421,638 (C)7→8 coding (520/858 nt) ydaU → Rac prophage; DUF1376 domain‑containing protein YdaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,421,6311.C94.1% 171.0 / 3.2 51coding (513/858 nt)ydaURac prophage; DUF1376 domain‑containing protein YdaU
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (29/19);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 5.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GAAGCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACC  >  NC_000913/1421492‑1421762
                                                                                                                                            |                                                                                                                                   
gaagCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGAt                                                                                                                                      <  2:1246376/140‑1 (MQ=255)
   gCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATccc                                                                                                                                  >  2:693293/1‑140 (MQ=255)
   gCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATccc                                                                                                                                  <  2:1196659/140‑1 (MQ=255)
     aGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATccccc                                                                                                                                <  1:754671/140‑1 (MQ=255)
           gAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATcccccccctaa                                                                                                                          >  1:1251181/1‑136 (MQ=255)
                   tGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATcccccc                                                                                                                   >  1:919113/1‑140 (MQ=255)
                   tGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATcccccc                                                                                                                   <  1:805864/140‑1 (MQ=255)
                      tCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCgg                                                                                                                <  2:1251181/140‑1 (MQ=255)
                                        aaaTTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAg                                                                                                      <  2:289610/132‑1 (MQ=255)
                                        aaaTTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAg                                                                                                      >  1:289610/1‑132 (MQ=255)
                                        aaaTTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCaaa                                                                                              <  1:1391214/140‑1 (MQ=255)
                                              aaCAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTg                                                                                       <  2:201392/140‑1 (MQ=255)
                                               aCAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTg                                                                                       <  2:319542/140‑1 (MQ=255)
                                               aCAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTg                                                                                       <  1:233149/140‑1 (MQ=255)
                                                  ggTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGAcc                                                                                    <  2:278294/140‑1 (MQ=255)
                                                   gtgtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGAccc                                                                                   >  2:110718/1‑140 (MQ=255)
                                                    tgtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                                  <  1:51151/140‑1 (MQ=255)
                                                    tgtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                                  <  1:1355309/140‑1 (MQ=255)
                                                    tgtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                                  <  1:693293/140‑1 (MQ=255)
                                                     gtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCctc                                                                                 >  2:1140710/1‑140 (MQ=255)
                                                      ttCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCctct                                                                                >  1:909853/1‑140 (MQ=255)
                                                      ttCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCctct                                                                                <  2:919113/140‑1 (MQ=255)
                                                      ttCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCctct                                                                                <  2:361502/140‑1 (MQ=255)
                                                      ttCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCctct                                                                                >  1:537933/1‑140 (MQ=255)
                                                      ttCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCctct                                                                                <  2:435529/140‑1 (MQ=255)
                                                      ttCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCctct                                                                                >  1:361502/1‑140 (MQ=255)
                                                       tCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTg                                                                               <  2:537933/140‑1 (MQ=255)
                                                           aCAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGata                                                                           <  2:666819/140‑1 (MQ=255)
                                                             aGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATAtt                                                                         >  1:1126909/1‑140 (MQ=255)
                                                               aTGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTAc                                                                       >  2:491205/1‑140 (MQ=255)
                                                                    gAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGc                                                                  >  2:88409/1‑140 (MQ=255)
                                                                     aaGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCg                                                                <  2:689711/140‑1 (MQ=255)
                                                                          gtgtTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACt                                                            >  1:1388992/1‑140 (MQ=255)
                                                                          gtgtTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACt                                                            >  2:663819/1‑140 (MQ=255)
                                                                               tGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGAtt                                                       >  1:923650/1‑140 (MQ=255)
                                                                                     tgatgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGtc                                                 <  1:663819/140‑1 (MQ=255)
                                                                                     tgatgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGtc                                                 <  2:1042443/140‑1 (MQ=255)
                                                                                         gTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGc                                             >  2:594050/1‑140 (MQ=255)
                                                                                                 aaaGGCAACTAATAAAGATACTGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTg                                    >  2:133628/1‑140 (MQ=255)
                                                                                                  aaGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTg                                    >  2:58072/1‑140 (MQ=255)
                                                                                                              aaaGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAAt                        >  2:1146694/1‑140 (MQ=255)
                                                                                                                   tACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGc                   >  2:942126/1‑140 (MQ=255)
                                                                                                                    aCAGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGcc                  >  2:347336/1‑140 (MQ=255)
                                                                                                                      aGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGt                             <  1:388243/127‑1 (MQ=255)
                                                                                                                      aGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGt                             >  2:388243/1‑127 (MQ=255)
                                                                                                                      aGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAg                >  2:440867/1‑140 (MQ=255)
                                                                                                                      aGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAg                >  2:387623/1‑140 (MQ=255)
                                                                                                                      aGATAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAg                >  1:794832/1‑140 (MQ=255)
                                                                                                                        aTAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTggg                              >  2:895657/1‑124 (MQ=255)
                                                                                                                        aTAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTggg                              <  1:895657/124‑1 (MQ=255)
                                                                                                                        aTAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGGCAGGc              >  1:911901/1‑140 (MQ=255)
                                                                                                                        aTAAAGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGc              >  2:421032/1‑140 (MQ=255)
                                                                                                                           aaGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCAtt           >  2:878260/1‑140 (MQ=255)
                                                                                                                            aGATCTAAAAACAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTg          >  2:804137/1‑140 (MQ=255)
                                                                                                                                  aaaaaCAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGaaaa    >  1:692098/1‑140 (MQ=255)
                                                                                                                                   aaaaCAGATCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAAcc  >  2:437835/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                   
GAAGCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACC  >  NC_000913/1421492‑1421762

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: