Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 574,486:1 +C 9.5% intergenic (+147/+43) quuD → / ← nmpC DLP12 prophage; putative antitermination protein/DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913574,4861.C 9.5% 175.2 / 14.4 42intergenic (+147/+43)quuD/nmpCDLP12 prophage; putative antitermination protein/DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage
Reads supporting (aligned to +/- strand):  ref base . (24/14);  new base C (2/2);  total (26/16)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

TTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCA  >  NC_000913/574348‑574632
                                                                                                                                           |                                                                                                                                                  
ttGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAAc                                                                                                                                          <  3:497950/149‑1 (MQ=255)
     tACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTc                                                                                                                                     <  2:31498/149‑1 (MQ=255)
        tAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGa                                                                                                                                                <  1:543110/135‑1 (MQ=255)
        tAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGa                                                                                                                                                >  2:543110/1‑135 (MQ=255)
            aaTTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTg                                                                                                                                                     <  3:660733/127‑1 (MQ=255)
                                  agaCGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGAtt                                                                                                        >  4:149813/1‑149 (MQ=255)
                                    aCGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAg                                                                                                      >  1:441041/1‑149 (MQ=255)
                                       cccccAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAAc                                                                                                   >  1:238228/1‑149 (MQ=255)
                                               aaTAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTa                                                                                              >  3:802375/1‑146 (MQ=255)
                                               aaTAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTa                                                                                              <  4:802375/146‑1 (MQ=255)
                                                           tatataACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGc                                                                               >  1:637739/1‑149 (MQ=255)
                                                              ataACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAAc                                                                            <  3:149813/149‑1 (MQ=255)
                                                              ataACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAAc                                                                            >  4:597246/1‑149 (MQ=255)
                                                                aaCAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTGATATGTTGCGGGAGGCGTTGCCTCCC‑CAACatataa                                                                                                                        <  1:545778/103‑3 (MQ=25)
                                                                aaCAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTGATATGTTGCGGGAGGCGTTGCCTCCC‑CAACatataa                                                                                                                        >  2:545778/1‑101 (MQ=25)
                                                                   aGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGt                                                                       >  1:341386/1‑149 (MQ=255)
                                                                            atataGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGc                                                                                                                                             <  4:4235/70‑1 (MQ=255)
                                                                            atataGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGc                                                                                                                                             >  3:4235/1‑70 (MQ=255)
                                                                                  ttAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTacac                                                        >  3:469969/1‑149 (MQ=255)
                                                                                  ttAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTacac                                                        <  2:637739/149‑1 (MQ=255)
                                                                                         cAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGc                                                 <  2:341386/149‑1 (MQ=255)
                                                                                              ttGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAg                                            <  3:597246/149‑1 (MQ=255)
                                                                                                  gCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGAt                                                                         <  2:670075/116‑1 (MQ=255)
                                                                                                  gCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGAt                                                                         >  1:670075/1‑116 (MQ=255)
                                                                                                    ttCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAgg                                      >  4:706250/1‑149 (MQ=255)
                                                                                                    ttCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAgg                                      >  2:672546/1‑149 (MQ=255)
                                                                                                           ggAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCtttt                                                                                                                     <  3:766362/63‑1 (MQ=255)
                                                                                                           ggAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCtttt                                                                                                                     >  4:766362/1‑63 (MQ=255)
                                                                                                                   cTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCgggg                       >  1:163102/1‑149 (MQ=255)
                                                                                                                          gAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCggg                        >  3:276634/1‑141 (MQ=255)
                                                                                                                          gAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCggg                        <  4:276634/141‑1 (MQ=255)
                                                                                                                              cACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCgg                         >  2:579961/1‑136 (MQ=255)
                                                                                                                              cACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCgg                         <  1:579961/136‑1 (MQ=255)
                                                                                                                              cACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTtctc            >  4:682780/1‑149 (MQ=255)
                                                                                                                               aCTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCg           >  3:511879/1‑149 (MQ=255)
                                                                                                                               aCTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCg           >  2:260721/1‑149 (MQ=255)
                                                                                                                                 agatatGTTGCGGGAGGCGTTGCCTCCC‑CAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAAtt                  <  3:467245/138‑1 (MQ=38)
                                                                                                                                 agatatGTTGCGGGAGGCGTTGCCTCCC‑CAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAAtt                  >  4:467245/3‑140 (MQ=38)
                                                                                                                                      tGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGCTGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACt    >  1:227604/1‑149 (MQ=255)
                                                                                                                                       gTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAgg                                           <  2:457898/109‑1 (MQ=255)
                                                                                                                                       gTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAgg                                           >  1:457898/1‑109 (MQ=255)
                                                                                                                                       gTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTc   >  4:519728/1‑149 (MQ=255)
                                                                                                                                        ttG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCa  >  3:661431/1‑149 (MQ=255)
                                                                                                                                          g‑gggAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAAc                                                                                           >  4:220028/1‑58 (MQ=255)
                                                                                                                                          g‑gggAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAAc                                                                                           <  3:220028/58‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
TTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCA  >  NC_000913/574348‑574632

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