Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,996 G→T 19.3% intergenic (+98/‑372) yffL → / → yffM CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,9960GT19.3% 97.2 / 4.9 48intergenic (+98/‑372)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (21/17);  new base T (6/4);  total (27/21)
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

AGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCG  >  NC_000913/2560848‑2561145
                                                                                                                                                    |                                                                                                                                                     
aGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCg                                                                                                                                                       <  6:263754/149‑1 (MQ=255)
    aGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTaa                                                                                                                                                   <  4:266444/149‑1 (MQ=255)
      tGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTCAAAAGTTGCCGTTAAgg                                                                                                                                                 >  3:235021/1‑149 (MQ=255)
           cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                                                            <  4:280547/149‑1 (MQ=255)
               tACCAGTTATTATTAAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCAGCCTTTTCTTTTCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGACTTTAAGGTGGGTCaa                                                                                                                                         >  8:281739/1‑148 (MQ=255)
                  cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATAGGTGGGTCAAAAGTAGCCTTTAAGGAGTGTCaaaatt                                                                                                                                     >  4:115004/1‑147 (MQ=255)
                  cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                                                                                                                             >  1:199995/1‑141 (MQ=255)
                  cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                                                                                                                             <  2:199995/141‑1 (MQ=255)
                   aGTTATTATTCAAGCCTTACAACATGAAATAACGGAATCAGGCTTGCAATGTTACAGACTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAATTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                                                                                                                    <  7:281739/149‑1 (MQ=255)
                   aGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                                                                                                                    <  1:123748/149‑1 (MQ=255)
                   aGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTGAATGGTGGGTCAAAAGTTGAGGTTAAGGAGGGTCAAAAttt                                                                                                                                    >  6:106767/1‑149 (MQ=255)
                        ttattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                                       <  2:8106/141‑1 (MQ=255)
                        ttattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                                       >  1:8106/1‑141 (MQ=255)
                          attCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTAAAATGTTGCCTTTTAGGTGGCTAATAAGTTGCCGtg                                                                                                                              >  5:335349/1‑147 (MQ=255)
                              aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGTGTCAAAAGTTGCAGTTAAGATGGGGCAAAAGTTGCCTATAAg                                                                                                                           >  7:101582/1‑147 (MQ=255)
                                gCCGTACAACATGAAATACAGGCATCAGTCTGGCAATGTTCCTGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTTTTGGGTCAAAAGTTGCAGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                       <  2:233671/149‑1 (MQ=38)
                                  cGTACATCATGAAATAAAGTTATAAGTAGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCAATGCATTCAGGTGGGTCAAAAGATGCCGTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                                     <  1:240774/149‑1 (MQ=25)
                                   gTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTTCTGGAGGGTCAAAAGATGCCTTTAAGGTGGGTCAAAATTTGTAGATAACGTggcg                                                                                                                     >  8:276467/1‑146 (MQ=255)
                                         caTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGGGGGGCAAAAGTTGCCGTTAAGGTGGGTGAAATGt                                                                                                              >  7:179608/1‑149 (MQ=255)
                                         caTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGACTTTAGTGGTGGGTAAAAATTTGACTTTTAGGATTGTCAAAATTTTGAGATAAGGTGGGTCAAATGt                                                                                                              >  3:252013/1‑149 (MQ=255)
                                                                         ccGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGAGGGTCAAAAGTTGGCGTGAAGGTGTGGGTAAATTTGAATAGGAAGGATTGTACAATTAGGAAtta                                                                               >  8:160678/1‑147 (MQ=21)
                                                                              ttttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCCGTTAAGGAGGATAAATATTTTCCGATAAGGTCGAGCAAAATTTTACGAGATAAAATTCTACAATAAATATTACTAgac                                                                          >  8:99686/1‑146 (MQ=16)
                                                                                         ttAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTTTTGTTGGGTAAAAATTTGCCTTTAAGGAGGCTAAAAATTTTAAGTTAAAGGGGCTGAAAATTTTACTAGTAAGAATTCTACAATTAAGATTTCTACGCGACGTGATtttt                                                              >  8:165228/1‑149 (MQ=255)
                                                                                            aTTCCATCAATGCAGGCGTGATAAAATCTGCCCGTAGTGCTGGGTAAAAAGTGGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                                           <  1:266199/149‑1 (MQ=25)
                                                                                                         aGGTGGGTCAAAAGCTGCCTTTCTTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtcg                                               <  7:74574/148‑1 (MQ=37)
                                                                                                               gTCAAAAGTTGCCGTTAGGGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCGACAATTAAGAATTCTAGGCTACGGTAGTTTATTAGCCTCCCGGTCGGTCa                                           >  4:205414/1‑146 (MQ=37)
                                                                                                                  caaaGTAGCCTTTAGTGGTGGGTCACAAGCTGCCTTTAAGGTGGGTCAAAAGTTGCCCTTAAGGTGGGGCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagaga                                     <  3:174359/148‑1 (MQ=255)
                                                                                                                         gCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTCGGTCAAAAGTTTCAGAGAAAGAATTCTGCTATTAAGAATTCGAGGGGACGTGAGTTATTTAGCCTCGCGGTAGGTGATAGAGAGTCggc                               >  2:43625/1‑146 (MQ=25)
                                                                                                                         gCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTAAAAAGTTGCAGTTAATGTCGATCAAAAGTTGCATAGTAAGCATTCGAGAATTAATAATTCTAGGCAACGGGAGTTTTTTATACTCGCGGTCGGTAATAGAtcgtcggc                               >  6:139685/1‑146 (MQ=19)
                                                                                                                          cccatcctgggggGTAAAAAGCTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGc                                <  7:134418/139‑1 (MQ=38)
                                                                                                                             ttAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTTACGAGTAAGAATTCGAGTATTAAGAATTCTAGGCGACGGGAGTTTTTTATACTCCCGGTCGGTCATAGATAGTCGCCTagcg                          >  2:218214/1‑147 (MQ=37)
                                                                                                                                 gggtggGTCAATAGTTGCCTTTAAGGTGGGTGAATTGTTGCGGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCTACGTGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaa                      <  3:38342/148‑1 (MQ=38)
                                                                                                                                      ggTCAAAAGTTGCCGTTAAGGTGTGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCg                                                                         <  7:163227/93‑1 (MQ=255)
                                                                                                                                        tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTATAGTTAAACag                  >  2:182384/1‑146 (MQ=255)
                                                                                                                                        tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag                  >  2:215280/1‑146 (MQ=255)
                                                                                                                                        tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag                  >  7:293912/1‑146 (MQ=255)
                                                                                                                                        tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag                  >  1:328489/1‑146 (MQ=255)
                                                                                                                                        tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGATCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag                  >  3:151173/1‑146 (MQ=255)
                                                                                                                                          aaaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGATTTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAgaagaa             >  1:139367/1‑149 (MQ=255)
                                                                                                                                          aaaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAgaagaa             >  6:360880/1‑149 (MQ=255)
                                                                                                                                            tacatGCCGTTTAGTTGGGTCAAAAGTTTCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAatt           <  7:276467/145‑1 (MQ=255)
                                                                                                                                            aaGTTCCCGTAAAGGTGGGTCAAAAGATTCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAatt           <  7:305264/149‑1 (MQ=255)
                                                                                                                                              gctGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAattat         <  7:160678/147‑1 (MQ=255)
                                                                                                                                              gTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCggtcggtc                                            <  7:111514/114‑1 (MQ=255)
                                                                                                                                              gTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCggtcggtc                                            >  8:111514/1‑114 (MQ=255)
                                                                                                                                               ttGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAattatt        <  5:360880/149‑1 (MQ=255)
                                                                                                                                               atGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAattatt        <  7:87052/148‑1 (MQ=255)
                                                                                                                                                  ccGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTc      <  2:328489/148‑1 (MQ=255)
                                                                                                                                                   cGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgc    >  7:41469/1‑149 (MQ=255)
                                                                                                                                                    gTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgccg  >  3:199766/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCG  >  NC_000913/2560848‑2561145

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: