Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,207,942 (G)7→8 coding (442/630 nt) stfP → e14 prophage, uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,207,9351.G95.8% 164.2 / ‑0.9 48coding (435/630 nt)stfPe14 prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (27/19);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGT  >  NC_000913/1207799‑1208065
                                                                                                                                          |                                                                                                                                   
GAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg                                                                                                                                  <  1:426825‑M1/139‑1 (MQ=255)
 AGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggg                                                                                                                                 >  2:698453‑M1/1‑139 (MQ=255)
        gCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggggggat                                                                                                                          <  2:653308/139‑4 (MQ=255)
               tCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATa                                                                                                                    <  2:218454/139‑1 (MQ=255)
                    gcagcaACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGtt                                                                                                               <  2:39033/139‑1 (MQ=255)
                      agccACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGCGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggggg                                                                                                                             >  2:38167/1‑122 (MQ=255)
                      agcaACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGt                                                                                                             >  2:474264/1‑139 (MQ=255)
                      agcaACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggggg                                                                                                                             <  1:38167/122‑1 (MQ=255)
                        caACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGACTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTAt                                                                                                           >  2:666549/1‑139 (MQ=255)
                                    tAGGACGCACCGATTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATAGAAATACAGTTGTATTCCCGGTTGCGt                                                                                               <  2:459869/139‑1 (MQ=255)
                                       gACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggggg                                                                                                                             <  2:376671/105‑1 (MQ=255)
                                       gACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggggg                                                                                                                             >  1:376671/1‑105 (MQ=255)
                                        aCGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGAGGGGGGATGGAAATAAAGTTGTATTCCCGGTTGCGTTtaa                                                                                           >  1:104849/1‑139 (MQ=255)
                                           cACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaataa                                                                                        >  2:93816/1‑139 (MQ=255)
                                            aCCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTaataat                                                                                       >  1:623311/1‑139 (MQ=255)
                                                               cGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTTGGATTGGAATTCCGCTAAAACACCTGCCTCTGGGGGGGGATGGAAATACAGTTGAAGTCTCGGTTTCGTTTAATAATGTCTGTGTTGCCGttgttg                                                                    <  2:104849/139‑1 (MQ=255)
                                                                 gctATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGttgttg                                                                    >  2:644505/3‑137 (MQ=255)
                                                                 gctATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGttgttg                                                                    <  1:644505/135‑1 (MQ=255)
                                                                   tATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGttgttg                                                                    >  2:281584/1‑135 (MQ=255)
                                                                   tATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGttgttg                                                                    <  1:281584/135‑1 (MQ=255)
                                                                      gATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTa                                                             >  2:53637/1‑139 (MQ=255)
                                                                         tggtggTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGa                                                          >  2:562045/1‑139 (MQ=255)
                                                                            tggtATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACaa                                                       <  2:426825/139‑1 (MQ=255)
                                                                            tggtATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACaa                                                       <  1:3267/139‑1 (MQ=255)
                                                                              gtATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg                                                                                                                                  >  1:105837/1‑61 (MQ=255)
                                                                              gtATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg                                                                                                                                  <  2:105837/61‑1 (MQ=255)
                                                                                     gCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgcgc                                              >  1:411180/1‑139 (MQ=255)
                                                                                      cAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgggg                                                                                                                                >  2:176424/1‑55 (MQ=255)
                                                                                      cAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgggg                                                                                                                                <  1:176424/55‑1 (MQ=255)
                                                                                        aTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCAcc                                           >  1:93931/1‑139 (MQ=255)
                                                                                         tCAGTAATCATTCAGTTGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCa                                                                  <  2:150188/115‑1 (MQ=255)
                                                                                         tCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCa                                                                  >  1:150188/1‑115 (MQ=255)
                                                                                         tCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCt                                          >  2:85883/1‑139 (MQ=255)
                                                                                         tCAGTAATCAGTCAGAGGGGGTGGAATGTCGCAAAAGCACCTGCCTCGGGGGGGGGATGGAAATACAGGGTTATTCCCGGTTGCGTTTAGTAATGCCTGGGTTGCCATTGATGCAATTTATGACAATTTTACCGCActt                                          >  2:630588/1‑137 (MQ=255)
                                                                                           ggTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTAt                                        >  2:481694/2‑139 (MQ=255)
                                                                                              aaTCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggggg                                                                                                                             <  2:95482/50‑1 (MQ=255)
                                                                                              aaTCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggggg                                                                                                                             >  1:95482/1‑50 (MQ=255)
                                                                                                   ttCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCgg                                                                                                     <  1:121070/70‑1 (MQ=255)
                                                                                                   ttCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCgg                                                                                                     >  2:121070/1‑70 (MQ=255)
                                                                                                   ttCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAg                                >  1:104034/1‑139 (MQ=255)
                                                                                                   ttCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggggg                                                                                                                             >  2:544373/1‑45 (MQ=255)
                                                                                                   ttCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggggg                                                                                                                             <  1:544373/45‑1 (MQ=255)
                                                                                                     cAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTg                              >  2:32340/1‑139 (MQ=255)
                                                                                                     cAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTg                              >  2:521160/1‑139 (MQ=255)
                                                                                                        tGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGTGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCa                           >  1:479568/1‑139 (MQ=255)
                                                                                                          ggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAAc                         <  1:93816/139‑1 (MQ=255)
                                                                                                          ggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAAc                         >  2:448595/1‑139 (MQ=255)
                                                                                                           gggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACg                        <  1:303576/139‑1 (MQ=255)
                                                                                                             gTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaat                                                                                          <  2:279625/71‑1 (MQ=255)
                                                                                                             gTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaat                                                                                          >  1:279625/1‑71 (MQ=255)
                                                                                                             gTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACggg                      <  1:562045/139‑1 (MQ=255)
                                                                                                             gTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACggg                      >  2:385761/1‑139 (MQ=255)
                                                                                                                     ccGCAAAAGCATCTGCCTCTGGGGGGGGATGGaaa                                                                                                                      <  2:478065/35‑1 (MQ=255)
                                                                                                                     ccGCAAAAGCATCTGCCTCTGGGGGGGGATGGaaa                                                                                                                      >  1:478065/1‑35 (MQ=255)
                                                                                                                     ccGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgc                                                <  1:26379/105‑1 (MQ=255)
                                                                                                                     ccGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgc                                                >  2:26379/1‑105 (MQ=255)
                                                                                                                     ccGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGAtata              >  1:61981/1‑139 (MQ=255)
                                                                                                                      cGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCaaa                                                                <  1:455705/88‑1 (MQ=255)
                                                                                                                      cGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATAc             >  1:405736/1‑139 (MQ=255)
                                                                                                                      cGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGATGCaaa                                                                >  2:455705/1‑88 (MQ=255)
                                                                                                                                tCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGc                            <  1:132727/114‑1 (MQ=255)
                                                                                                                                tCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGc                            >  2:132727/1‑114 (MQ=255)
                                                                                                                                 cTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGt  <  1:474264/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                   
GAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGT  >  NC_000913/1207799‑1208065

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