Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,422,015 A→G E11E (GAA→GAG ydaV → Rac prophage; putative ATP‑binding protein YdaV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,422,0150AG98.3% 214.5 / ‑4.4 59E11E (GAA→GAGydaVRac prophage; putative ATP‑binding protein YdaV
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (34/24);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

GTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAACGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTGAGGAGATCAACCGTCAGAATCGTCAGTTGCGGGTGGAAA  >  NC_000913/1421880‑1422151
                                                                                                                                       |                                                                                                                                        
gTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTa                                                                                                                                      >  2:1241640/1‑140 (MQ=255)
 ttCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTAt                                                                                                                                     >  1:710658/1‑140 (MQ=255)
  tCGCGCCGAAAGGTGTTCGCCCTGAGACGTTCCTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATc                                                                                                                                    <  1:181014/140‑1 (MQ=255)
     cgcCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGc                                                                                                                                 >  2:994991/1‑140 (MQ=255)
        cGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGa                                                                                                                              <  2:287303/140‑1 (MQ=255)
             ggTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGc                                                                                                                         <  2:78566/140‑1 (MQ=255)
              gtgtTCGACCTGAGACGTTCCTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGcc                                                                                                                        <  1:326986/140‑1 (MQ=255)
              gtgtTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGcc                                                                                                                        <  2:87865/140‑1 (MQ=255)
                  tCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGt                                                                                                                    <  1:5022/140‑1 (MQ=255)
                  tCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGt                                                                                                                    >  1:651234/1‑140 (MQ=255)
                   cGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTc                                                                                                                   <  2:1167458/140‑1 (MQ=255)
                         gagaCGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGt                                                                                                             >  1:654630/1‑140 (MQ=255)
                                tACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACcgcg                                                                                                      <  2:1228271/140‑1 (MQ=255)
                                 aCTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCgcgc                                                                                                     >  1:136942/1‑140 (MQ=255)
                                  ccccGACAGGTTAACCCCGTCTCGGTACCGGATAGTGGGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGGCTGGCCCCGTCACATGTAACCGCGcc                                                                                                    <  2:762097/138‑1 (MQ=255)
                                        cAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCaa                                                                                              <  1:1112609/140‑1 (MQ=255)
                                          ggTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTGAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGa                                                                                            >  1:1386108/1‑140 (MQ=255)
                                          ggTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGa                                                                                            <  2:405724/140‑1 (MQ=255)
                                          ggTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGa                                                                                            <  1:994991/140‑1 (MQ=255)
                                          ggTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGa                                                                                            <  1:1095963/140‑1 (MQ=255)
                                           gTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGAc                                                                                           >  2:1559750/1‑140 (MQ=255)
                                                     tctcGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGGCCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAg                                                                                 >  1:892119/1‑140 (MQ=255)
                                                     tctcGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAg                                                                                 >  2:719825/1‑140 (MQ=255)
                                                     tctcGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGCAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAg                                                                                 >  1:892466/1‑140 (MQ=255)
                                                     tctcGTTACCGGATAGTGCGATCCAGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAg                                                                                 >  2:1192059/1‑140 (MQ=255)
                                                            aCCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGcgcg                                                                          >  1:1055871/1‑140 (MQ=255)
                                                              cGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAg                                                                        >  2:1101937/1‑140 (MQ=255)
                                                                gATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTg                                                                      >  1:690723/1‑140 (MQ=255)
                                                                 aTAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTgg                                                                     >  1:453612/1‑140 (MQ=255)
                                                                  tAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGc                                                                    <  2:213749/140‑1 (MQ=255)
                                                                  tAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGc                                                                    <  2:892119/140‑1 (MQ=255)
                                                                    gTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCaa                                                                  >  1:25417/1‑140 (MQ=255)
                                                                     tGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAAc                                                                 >  2:741334/1‑140 (MQ=255)
                                                                     tGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAAc                                                                 >  1:723584/1‑140 (MQ=255)
                                                                       cGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACtt                                                               >  1:1327849/1‑140 (MQ=255)
                                                                           cccGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCg                                                           <  2:1180019/140‑1 (MQ=255)
                                                                            ccGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGa                                                          >  1:1581206/1‑140 (MQ=255)
                                                                             cGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGaa                                                         <  2:105906/140‑1 (MQ=255)
                                                                                cAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGc                                                      <  2:651234/140‑1 (MQ=255)
                                                                                        aGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGTCAGAAACg                                              >  1:1096285/1‑140 (MQ=255)
                                                                                          ggggTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGtt                                            <  2:1386108/140‑1 (MQ=255)
                                                                                          ggggTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGtt                                            <  2:25417/140‑1 (MQ=255)
                                                                                            ggTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTtgt                                          >  1:838619/1‑140 (MQ=255)
                                                                                              tAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTga                                        >  1:820734/1‑140 (MQ=255)
                                                                                               aaCTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTgag                                       >  1:1106299/1‑140 (MQ=255)
                                                                                                    aCCATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTGAGGAGAt                                  >  1:606043/1‑140 (MQ=255)
                                                                                                      cATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTGAGGAGATCa                                >  2:1412742/1‑140 (MQ=255)
                                                                                                      cATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTGAGGAGATCa                                <  2:136942/140‑1 (MQ=255)
                                                                                                      cATGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTGAGGAGATCa                                >  2:1030385/1‑140 (MQ=255)
                                                                                                        tGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTGAGGAGATCAAc                              >  2:638829/1‑140 (MQ=255)
                                                                                                        tGAAAAATATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTGAGGAGATCAAc                              <  2:723584/140‑1 (MQ=255)
                                                                                                          aaaaaTATTGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTGAGGAGATCAACCg                            >  2:1059613/1‑140 (MQ=255)
                                                                                                                 ttGCGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTGAGGAGATCAACCGTCAGAAt                     <  2:892466/140‑1 (MQ=255)
                                                                                                                    cGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTGAGGAGATCAACCGTCAGAATCGt                  >  1:1454950/1‑140 (MQ=255)
                                                                                                                    cGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTGAGGAGATCAACCGTCAGAATCGt                  >  2:1575994/1‑140 (MQ=255)
                                                                                                                    cGACAGGCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTGAGGAGATCAACCGTCAGAATCGt                  >  1:257115/1‑140 (MQ=255)
                                                                                                                    cGACAGGCGATGTTCTTGAACGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTGAGGAGATCAACCGTCAGAATCGt                  >  1:312192/1‑140 (MQ=255)
                                                                                                                         ggCGATGTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTGAGGAGATCAACCGTCAGAATCGTCAGtt             >  2:62128/1‑140 (MQ=255)
                                                                                                                               gTTCTTGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGTCAGAAACGTTGTGAGGAGATCAACCGTCAGAATCGTCAGTTGCGGGt       <  2:1096285/140‑1 (MQ=255)
                                                                                                                                    tGAGCGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTGAGGAGATCAACCGTCAGAATCGTCAGTTGCGGGTGGaaa  <  2:654630/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
GTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAATATTGCGACAGGCGATGTTCTTGAACGTATCCGCAGACTGGCCCCGTCACATGTAACCGCGCCATTCAAGACGGTAGCGGAGTGGCGCGAGTGGCAACTTTCCGAAGGCCAGAAACGTTGTGAGGAGATCAACCGTCAGAATCGTCAGTTGCGGGTGGAAA  >  NC_000913/1421880‑1422151

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