Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 580,851 G→A 23.2% K6K (AAG→AAA)  nohD → DLP12 prophage; DNA packaging protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,8510GA23.2% 115.8 / 30.1 64K6K (AAG→AAA) nohDDLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (18/31);  new base A (7/8);  total (25/39)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGC  >  NC_000913/580756‑580929
                                                                                               |                                                                              
aagGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagcagc                                                                            >  1:1283993/1‑100 (MQ=255)
   gAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctgg                                                                         <  2:1514672/100‑1 (MQ=255)
    aaaCGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAgctggc                                                                        <  2:2225480/100‑1 (MQ=255)
      aCGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggctg                                                                      <  2:2028895/100‑1 (MQ=255)
       cGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTTTTTGTTGTTTGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGa                                                                     >  1:3075032/1‑100 (MQ=17)
        gACAGGTGTTGAAAGCGAGCTTTTTGGGCTCTGTCGTTTCCTTTCTCTGATTTTGTCCGTGGAATGTGTTATGGAAGTCAACAAAAAGCAgctggctg                                                                      <  1:2694838/98‑1 (MQ=17)
         aCAGGTGATGAAAGCGAGCTTTTTGGCCTCTGTCGTTTACTTTCTCTATGTTTGGCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGa                                                                     <  1:2965454/98‑1 (MQ=14)
            ggTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAttt                                                                >  2:1285582/1‑100 (MQ=255)
             gTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTCTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAt                                                                  <  1:1825288/97‑1 (MQ=255)
              tGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTc                                                              <  2:2084692/100‑1 (MQ=255)
              tGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTc                                                              <  1:1883252/100‑1 (MQ=255)
              tGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAACGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTc                                                              <  1:1852448/100‑1 (MQ=255)
               gCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCg                                                             <  2:1283997/100‑1 (MQ=255)
                 tGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGt                                                           <  1:3083117/100‑1 (MQ=255)
                  gAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTg                                                          >  1:756248/1‑100 (MQ=255)
                  gAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTg                                                          <  1:569147/100‑1 (MQ=255)
                   aaaGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGc                                                         <  2:458529/100‑1 (MQ=255)
                      gCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTGCGTGGAGTGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTGGGTGGGAg                                                      >  1:2955871/1‑100 (MQ=14)
                      gCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAg                                                      <  2:2943567/100‑1 (MQ=255)
                        gAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAg                                                      <  1:518002/98‑1 (MQ=255)
                          gCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATc                                                  >  2:1569222/1‑100 (MQ=255)
                          gCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATc                                                  >  1:1612004/1‑100 (MQ=255)
                           cTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATcc                                                 >  1:938115/1‑100 (MQ=255)
                             ctttGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCTGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGt                                               <  2:3227943/99‑1 (MQ=255)
                                tGGGCTCTGCTGTTTCCTTTCTCTGTTGTTGTCCGTGGAATGTGGAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTAc                                             <  1:2937179/99‑1 (MQ=14)
                                tGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTAcc                                            >  1:3335595/1‑100 (MQ=255)
                                tGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTAcc                                            <  2:3302362/100‑1 (MQ=255)
                                   gctctGTCGTTTCCTTTCTCTGTTTTTGTCCGTGCAATTTTCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCAt                                          <  1:1082394/98‑1 (MQ=11)
                                   ccTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCAtt                                         >  1:3139315/1‑100 (MQ=255)
                                   ccTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCAtt                                         >  2:2075824/1‑100 (MQ=255)
                                   ccTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCAtt                                         >  2:2043744/1‑100 (MQ=255)
                                   ccTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCAtt                                         <  2:3020267/100‑1 (MQ=255)
                                    ctctGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTc                                        <  1:2043730/100‑1 (MQ=255)
                                     tctGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCa                                       <  2:1059516/100‑1 (MQ=255)
                                     tctGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCa                                       >  1:2999308/1‑100 (MQ=255)
                                      ctGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAg                                      <  2:3241176/100‑1 (MQ=255)
                                      ctGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAg                                      >  1:736744/1‑100 (MQ=255)
                                      ctGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAg                                      >  1:3309802/1‑100 (MQ=255)
                                       tGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGa                                     >  1:1344694/1‑100 (MQ=255)
                                          cGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACt                                  >  2:2938917/1‑100 (MQ=255)
                                           gTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTg                                 >  2:2824085/1‑100 (MQ=255)
                                             ttCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGc                               >  1:1730641/1‑100 (MQ=255)
                                              tCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCa                              <  1:2075811/100‑1 (MQ=255)
                                               ccTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAg                             >  2:82113/1‑100 (MQ=255)
                                                 tttattaactttaGTACGTGCAATATGCTTTGGATTTCTACAAAAAGCAGCTGGCTGACATTTTCGGTCCGAGTATCCGTACCATTCAGCACTGGCCgg                            >  1:75286/10‑99 (MQ=2)
                                                 tttCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGa                           <  2:3364726/100‑1 (MQ=40)
                                                    ctctGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACa                        >  1:2492927/1‑100 (MQ=255)
                                                     tctGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAg                       <  2:867030/100‑1 (MQ=255)
                                                     tctGTTTTTGTCAGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAg                       >  1:2337734/1‑100 (MQ=255)
                                                     tctGTTTGTGTGCGTGGAATGTGCGATGGAGGTCAACAAAAAGCAGCTGGCTGAGCTTTGGGGTGCGAGTATCGGTACCATTCAGATCTGGCAg                             <  1:617960/94‑1 (MQ=14)
                                                       tGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAggg                     <  2:83285/100‑1 (MQ=255)
                                                       tGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAggg                     <  2:87457/100‑1 (MQ=255)
                                                       tGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAggg                     <  1:1745793/100‑1 (MQ=255)
                                                       tGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAggg                     <  2:1231498/100‑1 (MQ=40)
                                                       tGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAgg                      >  1:770019/1‑99 (MQ=40)
                                                         tttttGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGaa                   >  2:1236432/1‑100 (MQ=255)
                                                          ttttGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAAt                  <  2:589505/100‑1 (MQ=40)
                                                                gagcggAATGTGCAATTGAAGTCAAAAAACCGCTTCTGGCTGACCTTTTCGGTGAGAGTATCCGTACCATTCAGAACTGGCAGGGACAggg                     <  1:2019327/87‑1 (MQ=11)
                                                                 cGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGtt           <  2:1219629/100‑1 (MQ=38)
                                                                  gTGGACTGTGCAATGCAAGTAAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTc          >  1:1583938/1‑100 (MQ=255)
                                                                  gTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTc          <  2:2520068/100‑1 (MQ=38)
                                                                  gTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTc          <  2:1424007/100‑1 (MQ=35)
                                                                   tGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCt         <  1:963480/100‑1 (MQ=35)
                                                                       aTGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGa     <  1:2250198/100‑1 (MQ=35)
                                                                          tgCAATGGAAGTCAACAAAAAGCAGCGGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGc  <  1:1641695/100‑1 (MQ=14)
                                                                                               |                                                                              
AAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGC  >  NC_000913/580756‑580929

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