Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 574,494 A→G 12.2% intergenic (+155/+35) 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,4940AG12.2% 139.6 / 9.1 49intergenic (+155/+35)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 A (28/15);  new base G (3/3);  total (31/18)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCAT  >  NC_000913/574348‑574642
                                                                                                                                                   |                                                                                                                                                    
ttGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAAc                                                                                                                                                    <  3:497950/149‑1 (MQ=255)
     tACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTc                                                                                                                                               <  2:31498/149‑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)
                                                                                  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                                                                                   >  1:670075/1‑116 (MQ=255)
                                                                                                  gCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGAt                                                                                   <  2:670075/116‑1 (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                     >  2:260721/1‑149 (MQ=255)
                                                                                                                               aCTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCg                     >  3:511879/1‑149 (MQ=255)
                                                                                                                                 agatatGTTGCGGGAGGCGTTGCCTCCC‑CAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAAtt                            >  4:467245/3‑140 (MQ=38)
                                                                                                                                 agatatGTTGCGGGAGGCGTTGCCTCCC‑CAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAAtt                            <  3:467245/138‑1 (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)
                                                                                                                                        ttgcgggAGGCGTTGCCTCCC‑CAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGc                                                           >  3:657233/5‑102 (MQ=37)
                                                                                                                                        ttgcgggAGGCGTTGCCTCCC‑CAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGc                                                           <  4:657233/98‑1 (MQ=37)
                                                                                                                                        ttG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCa            >  3:661431/1‑149 (MQ=255)
                                                                                                                                          g‑gggAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAAc                                                                                                     >  4:220028/1‑58 (MQ=255)
                                                                                                                                          g‑gggAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAAc                                                                                                     <  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     <  3:393767/148‑1 (MQ=255)
                                                                                                                                                 gCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcattt     >  4:393767/1‑148 (MQ=255)
                                                                                                                                                   aaCGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttcat  >  2:785228/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTG‑GGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCAT  >  NC_000913/574348‑574642

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: ‑