Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 266,505 A→G F17L (TTT→CTT)  yafY ← CP4‑6 prophage; inner membrane lipoprotein YafY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913266,5050AG98.1% 181.2 / ‑4.1 52F17L (TTT→CTT) yafYCP4‑6 prophage; inner membrane lipoprotein YafY
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (30/21);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.97e-01

TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGT  >  NC_000913/266368‑266637
                                                                                                                                         |                                                                                                                                    
tAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACa                                                                                                                                    <  2:462251/140‑1 (MQ=255)
   gtagtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCagag                                                                                                                                 <  2:628185/140‑1 (MQ=255)
      gtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCa                                                                                                                                    <  1:381543/134‑1 (MQ=255)
      gtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCa                                                                                                                                    >  2:381543/1‑134 (MQ=255)
         aaCCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCagag                                                                                                                                 <  2:64515/134‑1 (MQ=255)
         aaCCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCagag                                                                                                                                 >  1:64515/1‑134 (MQ=255)
           ccTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCa                                                                                                                         >  1:933050/1‑140 (MQ=255)
               tgACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCa                                                                                                                     >  1:827642/1‑140 (MQ=255)
                gACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAg                                                                                                                    <  2:133515/140‑1 (MQ=255)
                  ccATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCACCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAgcg                                                                                                                  <  1:731247/140‑1 (MQ=255)
                     ttcttcGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCca                                                                                                               >  2:157354/1‑140 (MQ=255)
                                catcatCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAgt                                                                                                    >  2:183052/1‑140 (MQ=255)
                                  tcatcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAgtgt                                                                                                  >  2:758236/1‑140 (MQ=255)
                                    atcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGc                                                                                                       >  1:634420/1‑133 (MQ=255)
                                    atcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGc                                                                                                       <  2:634420/133‑1 (MQ=255)
                                             gCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCa                                                                                       >  1:846883/1‑140 (MQ=255)
                                              cTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCat                                                                                      >  1:709578/1‑140 (MQ=255)
                                              cTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCat                                                                                      >  1:479217/1‑140 (MQ=255)
                                              cTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCat                                                                                      >  1:162663/1‑140 (MQ=255)
                                                gACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCatta                                                                                    >  2:121396/1‑140 (MQ=255)
                                                 aCCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCattat                                                                                   >  2:468506/1‑140 (MQ=255)
                                                                cgcgGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATc                                                                    >  1:689858/1‑140 (MQ=255)
                                                                cgcgGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATc                                                                    >  1:821306/1‑140 (MQ=255)
                                                                cgcgGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATc                                                                    >  1:361856/1‑140 (MQ=255)
                                                                  cgGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAg                                                                  >  2:391260/1‑140 (MQ=255)
                                                                   ggTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGa                                                                 >  1:424351/1‑140 (MQ=255)
                                                                        ggTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAAc                                                            <  1:121396/140‑1 (MQ=255)
                                                                                 cTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGc                                                   <  1:758236/140‑1 (MQ=255)
                                                                                  ttAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCa                                                  <  2:306206/140‑1 (MQ=255)
                                                                                   tAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCaa                                                 <  2:827642/140‑1 (MQ=255)
                                                                                    aaaTTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCaaa                                                <  2:844370/140‑1 (MQ=255)
                                                                                     aaTTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAg                                               <  2:479217/140‑1 (MQ=255)
                                                                                         tACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACAt                                           >  2:718359/1‑140 (MQ=255)
                                                                                         tACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACAt                                           <  1:468506/140‑1 (MQ=255)
                                                                                              aTGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATg                                      >  1:867108/1‑140 (MQ=255)
                                                                                                ggCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGcc                                    <  2:424351/140‑1 (MQ=255)
                                                                                                     tCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGa                               <  1:718359/140‑1 (MQ=255)
                                                                                                            atcaCTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATAcc                        >  2:114882/1‑140 (MQ=255)
                                                                                                            atcaCTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATAcc                        >  1:946822/1‑140 (MQ=255)
                                                                                                            atcaCTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATAcc                        >  1:704137/1‑140 (MQ=255)
                                                                                                            atcaCTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATAcc                        >  1:33376/1‑140 (MQ=255)
                                                                                                               aCTTCTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCAc                     <  1:157354/140‑1 (MQ=255)
                                                                                                                   cTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGc                 >  1:558771/1‑140 (MQ=255)
                                                                                                                   cTGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGc                 >  2:733048/1‑140 (MQ=255)
                                                                                                                    tGTCATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCa                <  2:33376/140‑1 (MQ=255)
                                                                                                                       cATCGCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGc             >  2:407466/1‑140 (MQ=255)
                                                                                                                           gCATGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCcgac         >  2:441201/1‑140 (MQ=255)
                                                                                                                              tGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGaa      <  2:704137/140‑1 (MQ=255)
                                                                                                                              tGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGaa      <  1:183785/140‑1 (MQ=255)
                                                                                                                              tGCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGaa      <  2:933050/140‑1 (MQ=255)
                                                                                                                               gCGATGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAAc     <  1:114882/140‑1 (MQ=255)
                                                                                                                                  aTGAGAAGCAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGgt  >  2:536954/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                    
TAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGT  >  NC_000913/266368‑266637

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