Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 574,496 C→T 12.0% intergenic (+157/+33) 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,4960CT12.0% 144.0 / 8.9 50intergenic (+157/+33)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 C (29/15);  new base T (3/3);  total (32/18)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATT  >  NC_000913/574348‑574643
                                                                                                                                                     |                                                                                                                                                   
ttGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAAc                                                                                                                                                     <  3:497950/149‑1 (MQ=255)
     tACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTc                                                                                                                                                <  2:31498/149‑1 (MQ=255)
                                  agaCGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGAtt                                                                                                                   >  4:149813/1‑149 (MQ=255)
                                    aCGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAg                                                                                                                 >  1:441041/1‑149 (MQ=255)
                                       cccccAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAAc                                                                                                              >  1:238228/1‑149 (MQ=255)
                                               aaTAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTa                                                                                                         <  4:802375/146‑1 (MQ=255)
                                               aaTAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTa                                                                                                         >  3:802375/1‑146 (MQ=255)
                                                           tatataACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGc                                                                                          >  1:637739/1‑149 (MQ=255)
                                                              ataACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAAc                                                                                       >  4:597246/1‑149 (MQ=255)
                                                              ataACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAAc                                                                                       <  3:149813/149‑1 (MQ=255)
                                                                aaCAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTGATATGTTGCGGGAGGCGTTGCCTCCCCAACatataa                                                                                                                                   >  2:545778/1‑101 (MQ=25)
                                                                aaCAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTGATATGTTGCGGGAGGCGTTGCCTCCCCAACatataa                                                                                                                                   <  1:545778/103‑3 (MQ=25)
                                                                   aGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGt                                                                                  >  1:341386/1‑149 (MQ=255)
                                                                                  ttAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTacac                                                                   <  2:637739/149‑1 (MQ=255)
                                                                                  ttAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTacac                                                                   >  3:469969/1‑149 (MQ=255)
                                                                                         cAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGc                                                            <  2:341386/149‑1 (MQ=255)
                                                                                              ttGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAg                                                       <  3:597246/149‑1 (MQ=255)
                                                                                                  gCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGAt                                                                                    <  2:670075/116‑1 (MQ=255)
                                                                                                  gCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGAt                                                                                    >  1:670075/1‑116 (MQ=255)
                                                                                                    ttCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAgg                                                 >  4:706250/1‑149 (MQ=255)
                                                                                                    ttCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAgg                                                 >  2:672546/1‑149 (MQ=255)
                                                                                                           ggAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCtttt                                                                                                                                >  4:766362/1‑63 (MQ=255)
                                                                                                           ggAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCtttt                                                                                                                                <  3:766362/63‑1 (MQ=255)
                                                                                                                   cTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCgggg                                  >  1:163102/1‑149 (MQ=255)
                                                                                                                          gAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCggg                                   <  4:276634/141‑1 (MQ=255)
                                                                                                                          gAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCggg                                   >  3:276634/1‑141 (MQ=255)
                                                                                                                              cACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCgg                                    <  1:579961/136‑1 (MQ=255)
                                                                                                                              cACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCgg                                    >  2:579961/1‑136 (MQ=255)
                                                                                                                              cACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTtctc                       >  4:682780/1‑149 (MQ=255)
                                                                                                                               aCTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCg                      >  3:511879/1‑149 (MQ=255)
                                                                                                                               aCTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCg                      >  2:260721/1‑149 (MQ=255)
                                                                                                                                 agatatGTTGCGGGAGGCGTTGCCTCCCCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAAtt                             <  3:467245/138‑1 (MQ=38)
                                                                                                                                 agatatGTTGCGGGAGGCGTTGCCTCCCCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAAtt                             >  4:467245/3‑140 (MQ=38)
                                                                                                                                      tGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGCTGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACt               >  1:227604/1‑149 (MQ=255)
                                                                                                                                       gTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAgg                                                      <  2:457898/109‑1 (MQ=255)
                                                                                                                                       gTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAgg                                                      >  1:457898/1‑109 (MQ=255)
                                                                                                                                       gTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTc              >  4:519728/1‑149 (MQ=255)
                                                                                                                                        ttgcgggAGGCGTTGCCTCCCCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGc                                                            >  3:657233/5‑102 (MQ=37)
                                                                                                                                        ttgcgggAGGCGTTGCCTCCCCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGc                                                            <  4:657233/98‑1 (MQ=37)
                                                                                                                                        ttGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCa             >  3:661431/1‑149 (MQ=255)
                                                                                                                                          ggggAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAAc                                                                                                      >  4:220028/1‑58 (MQ=255)
                                                                                                                                          ggggAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAAc                                                                                                      <  3:220028/58‑1 (MQ=255)
                                                                                                                                            gggAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCa             <  3:502878/146‑1 (MQ=255)
                                                                                                                                            gggAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCa             >  4:502878/1‑146 (MQ=255)
                                                                                                                                              gAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcat        >  2:681793/1‑149 (MQ=255)
                                                                                                                                                 gCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttc     >  4:272140/1‑149 (MQ=255)
                                                                                                                                                 gCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttc     >  1:391946/1‑149 (MQ=255)
                                                                                                                                                 gCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcattt      >  4:393767/1‑148 (MQ=255)
                                                                                                                                                 gCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcattt      <  3:393767/148‑1 (MQ=255)
                                                                                                                                                   aaCGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttcat   >  2:785228/1‑149 (MQ=255)
                                                                                                                                                    aCGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttcatt  <  2:238228/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATT  >  NC_000913/574348‑574643

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: