Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,418,386 G→A intergenic (‑157/‑285) kilR ← / → sieB killing protein, Rac prophage, FtsZ inhibitor protein/phage superinfection exclusion protein, Rac prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,418,3860GA98.4% 198.2 / NA 64intergenic (‑157/‑285)kilR/sieBkilling protein, Rac prophage, FtsZ inhibitor protein/phage superinfection exclusion protein, Rac prophage
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (34/29);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCGGTGTCTGGTGCCTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCGCATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAGGCGCACTCCCACCATTTAAATTTAA  >  NC_000913/1418263‑1418515
                                                                                                                           |                                                                                                                                 
aTGCGGTGTCTGGTGCCTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTTCGCTTAGCCGCATTCACCACATCACAAAATTCAc                                                                                                                    >  1:421186/1‑139 (MQ=255)
      tgtCTGGTGCCTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTaaa                                                                                                              <  2:170293/139‑1 (MQ=255)
       gtCTGGTGCCTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACCTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTaaaa                                                                                                             <  2:421186/139‑1 (MQ=255)
           ggTGCCTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAggg                                                                                                         <  1:614121/139‑1 (MQ=255)
                cTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAAtt                                                                                                                       >  1:1128266/1‑120 (MQ=255)
                cTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAAtt                                                                                                                       <  2:1128266/120‑1 (MQ=255)
                      gTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCTCATTCACCACATCACAAAAtt                                                                                                                       <  1:1549883/114‑1 (MQ=255)
                      gTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAAtt                                                                                                                       >  2:1549883/1‑114 (MQ=255)
                                          aTTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACt                                                                          >  2:1546788/1‑139 (MQ=255)
                                           atACAGTCGGCTTTCCCACCCAAACCAATACGGACTAACATGACTTTTAACTGTGCCCCGTGCGCTTAGCCGCACTCTCCGCATCACAAAATTTACTTTTCAAACGGCGGCCATCAGCCGAACTTTAAGAAAAAAACTg                                                                         >  1:880168/2‑139 (MQ=255)
                                                 tCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATgccg                                                                   >  2:946045/1‑139 (MQ=255)
                                                  cGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATgccgc                                                                  <  1:1546788/139‑1 (MQ=255)
                                                   ggCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATgccgcc                                                                 >  2:1501318/1‑139 (MQ=255)
                                                             cccaAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTacac                                                       >  2:149190/1‑139 (MQ=255)
                                                                  aCCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACAt                                                                                                  <  1:359237/91‑1 (MQ=255)
                                                                  aCCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACAt                                                                                                  >  2:359237/1‑91 (MQ=255)
                                                                  aCCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATgccgc                                                                  <  2:1070032/123‑1 (MQ=255)
                                                                  aCCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATgccgc                                                                  >  1:1070032/1‑123 (MQ=255)
                                                                      aTAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATg                                              >  1:1017475/1‑139 (MQ=255)
                                                                       tAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGt                                             >  1:884250/1‑139 (MQ=255)
                                                                       tAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGt                                             >  1:1367288/1‑139 (MQ=255)
                                                                            aCTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTaaaa                                                                                                             <  1:1592034/70‑1 (MQ=255)
                                                                            aCTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTaaaa                                                                                                             >  2:1592034/1‑70 (MQ=255)
                                                                            aCTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATgccgc                                                                  <  1:200139/113‑1 (MQ=255)
                                                                            aCTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATgccgc                                                                  >  2:200139/1‑113 (MQ=255)
                                                                              tAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATc                                                                                                                               <  2:1306389/50‑1 (MQ=255)
                                                                              tAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATc                                                                                                                               >  1:1306389/1‑50 (MQ=255)
                                                                              tAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTAtt                                      >  1:190597/1‑139 (MQ=255)
                                                                                aCATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGAt                                                                       <  2:832805/104‑1 (MQ=255)
                                                                                aCATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGAt                                                                       >  1:832805/1‑104 (MQ=255)
                                                                                aCATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGAt                                                                       <  2:1606748/104‑1 (MQ=255)
                                                                                aCATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGAt                                                                       >  1:1606748/1‑104 (MQ=255)
                                                                                     aCTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATgcc                                                                    <  2:1461528/102‑1 (MQ=255)
                                                                                     aCTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATgcc                                                                    >  1:1461528/1‑102 (MQ=255)
                                                                                     aCTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAAcc                               >  1:73267/1‑139 (MQ=255)
                                                                                         ttAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCa                                                                <  2:588903/102‑1 (MQ=255)
                                                                                         ttAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCa                                                                >  1:588903/1‑102 (MQ=255)
                                                                                          tAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATc                                                                                                                               >  2:351738/1‑38 (MQ=255)
                                                                                          tAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATc                                                                                                                               <  1:351738/38‑1 (MQ=255)
                                                                                          tAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAAtt                                                                                                                       >  2:960292/1‑46 (MQ=255)
                                                                                          tAACTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAAtt                                                                                                                       <  1:960292/46‑1 (MQ=255)
                                                                                            aCTGTGCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAGgcg                        <  1:149190/139‑1 (MQ=255)
                                                                                                tgCCACGTGCGCTTAGCCGCATTCACCACATCACAAAAt                                                                                                                        <  2:720950/39‑1 (MQ=14)
                                                                                                tgCCACGTGCGCTTAGCCGCATTCACCACATCACAAAAt                                                                                                                        >  1:720950/1‑39 (MQ=14)
                                                                                                tgCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGa                                                             >  2:845873/1‑98 (MQ=255)
                                                                                                tgCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGa                                                             <  1:845873/98‑1 (MQ=255)
                                                                                                tgCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTaca                                                        >  2:797465/1‑103 (MQ=255)
                                                                                                tgCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTaca                                                        <  1:797465/103‑1 (MQ=255)
                                                                                                tgCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAGgcg                        <  2:1019583/135‑1 (MQ=255)
                                                                                                tgCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAGgcg                        >  1:1019583/1‑135 (MQ=255)
                                                                                                 gCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAgg                          >  1:422670/1‑132 (MQ=255)
                                                                                                 gCCACGTGCGCTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAgg                          <  2:422670/132‑1 (MQ=255)
                                                                                                        gcgcTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATgccgc                                                                  <  1:1600334/85‑1 (MQ=255)
                                                                                                        gcgcTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATgccgc                                                                  >  2:1600334/1‑85 (MQ=255)
                                                                                                          gcTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATgccgc                                                                  <  2:921609/83‑1 (MQ=255)
                                                                                                          gcTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATgccgc                                                                  >  1:921609/1‑83 (MQ=255)
                                                                                                          gcTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCa                                                 >  1:1590858/1‑100 (MQ=255)
                                                                                                          gcTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCa                                                 <  2:1590858/100‑1 (MQ=255)
                                                                                                           cTTAGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAGGCGCACTCCCACCattta         >  1:1129837/1‑139 (MQ=255)
                                                                                                              aGCCGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAGGCGCACTCCCACCatttaaat      >  1:1012755/1‑139 (MQ=255)
                                                                                                                 cGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAGGCGCACTccc                 >  2:829721/1‑125 (MQ=255)
                                                                                                                 cGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAGGCGCACTccc                 <  1:829721/125‑1 (MQ=255)
                                                                                                                 cGCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAGGCGCACTCCCACCatttaaattta   <  1:1501318/139‑1 (MQ=255)
                                                                                                                  gCATTCACCACATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAGGCGCACTCCCACCatttaaatttaa  <  2:1367288/139‑1 (MQ=255)
                                                                                                                           |                                                                                                                                 
ATGCGGTGTCTGGTGCCTCCAGGTGACTGCAACCAGTTAACAATTACAGTCGGCTTTCCCACCCAAACCAATAAGGACTAACATGACTTTTAACTGTGCCACGTGCGCTTAGCCGCATTCACCGCATCACAAAATTCACTTTAAAAAGGGCGGACATCAGCCGAACTTCAAGAAAAAAACTGATGCCGCCAGGACTACACACAGCAATGTCGTTATTTACAACCGGAGGCGCACTCCCACCATTTAAATTTAA  >  NC_000913/1418263‑1418515

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