Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,991 T→A 18.1% intergenic (+93/‑377) 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,9910TA18.1% 100.2 / 5.1 52intergenic (+93/‑377)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (24/18);  new base A (3/6);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAG  >  NC_000913/2560845‑2561129
                                                                                                                                                   |                                                                                                                                           
gAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGG‑TCAAAATTTg                                                                                                                                           >  2:195879/1‑149 (MQ=255)
gAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGGGGG‑TCAAAAGTTg                                                                                                                                           <  6:206373/149‑1 (MQ=255)
   aGGGAGTGCTACGAATACCTGTTATTATTCAATCCGTACAACATGAATTAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCGTTGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAATGTTGcc                                                                                                                                         <  5:85351/148‑1 (MQ=255)
     ggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGTTGGG‑TCAAAAGTTGCAGt                                                                                                                                       >  5:87320/1‑148 (MQ=255)
         tGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑Agg                                                                                                                                 >  6:362613/1‑149 (MQ=255)
           cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGt                                                                                                                                <  5:362613/148‑1 (MQ=255)
           cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑GCAAAAGTTGCCGTTA‑AGGt                                                                                                                                <  1:399641/148‑1 (MQ=255)
            aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTgg                                                                                                                              >  1:356036/1‑149 (MQ=255)
              cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGt                                                                                                                            <  8:18518/149‑1 (MQ=255)
              cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGt                                                                                                                            <  7:147270/149‑1 (MQ=255)
                   aCCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGTGCTATGTTCGCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATTGTGATGGG‑GCAAAAGATGGCGTTA‑AGGGGTGTGaaat                                                                                                                       >  4:282561/1‑148 (MQ=255)
                   aCCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTACATGGATTCAGGTGTGTCAAAAGTTGACGTTATTTCTGTG‑TAAAAAGTTGAAGTTA‑AGgtgtgtaatat                                                                                                                       >  2:86000/1‑144 (MQ=255)
                     cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCaaaatt                                                                                                                     >  5:288752/1‑147 (MQ=255)
                        ttattaTTCAAGCCGTACAACATGAAATAACGGAATCAGGCGTGCAATGTTCTAGCCTTTTCTTTTCCTTAATTACATGAATTCAGGTGGGTCAAAAGTTGACGTTAGTGGTGTT‑ACAAAAGTTGCCGTCA‑AGGTGGGTCTACAGTTg                                                                                                                   >  7:350182/1‑148 (MQ=255)
                        ttattaTTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGc                                                                                                                  <  7:11988/149‑1 (MQ=255)
                        ttattaTTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTAAGGTGGGTCAAAAGTTGCCGTTTATGGTTTG‑TAATATTTTGACTTTA‑AGGTGGGTGaatttttt                                                                                                                   >  3:308487/1‑142 (MQ=255)
                        ttattaTTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTAAGGTGGGTCAAAAGTTGCCGTTAGTGATGTG‑TAAAAATTTGCAGTTA‑AGGAGTGTCAAGAGTTg                                                                                                                   >  5:336177/1‑148 (MQ=255)
                          attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCaaaa                                                                                                                       <  1:11061/142‑1 (MQ=255)
                          attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAACAGTTGcgg                                                                                                                >  5:211267/1‑147 (MQ=255)
                          attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTACCCTTAATTCCATGAATTCAGGTGGCTCAAAAGTTTCCGTAAGTTGTGGG‑TCAAAAGATGCCGTTA‑AGGTGGGTCaaaa                                                                                                                       >  2:11061/1‑142 (MQ=255)
                           ttattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAATGTTGCCTTAA‑AGGTGGGTCAAAATTTg                                                                                                                   >  1:205995/1‑145 (MQ=255)
                           ttattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCATAAATTACATGAATTCAGGTCGGTCAAAAGTTGCCGTTAGTGGTGTG‑TCAAATTATGCCGTAA‑TGGGGGGTCAATAGTTgactg                                                                                                               >  5:334095/1‑145 (MQ=255)
                                 aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TAAAAATTTGCAGTTA‑AGGTGGGTCAAAAGTTGCAGTTTAggg                                                                                                         >  6:107979/1‑148 (MQ=255)
                                 aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTGAGTGGTGGG‑TCAAAAGTTGCCGTTT‑AGGTGGGTCAAAAGTTGCCgagtaggg                                                                                                         >  5:327987/1‑142 (MQ=255)
                                   gcagTACACCATGAAATAAAGGAATCAGTAGTGCAATTTTCACGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                       <  4:328964/146‑1 (MQ=37)
                                                       ggAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAATTTGCCGTTA‑AGGTGGGTCAAAAGTTGCAGTTAAGGTGGGTCAAAAGTTGCCGAGTAAg                                                                                   >  8:45963/1‑149 (MQ=255)
                                                              ggcggcCACTTTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGTGGCGTCAAAAGTTGCCGTTAGTGGTGGG‑GCAAAAGTTGCCTTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTcta                                                                            <  5:108021/149‑1 (MQ=37)
                                                              ggCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TAAAAAGTTGCGGTGA‑AGGTGGCGCAAAATTTGCCTTTAATGTGGGGGAAAATTTTCCTAGTTAGAATCctc                                                                            >  3:23727/1‑148 (MQ=255)
                                                              ggCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTACATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGGGGG‑TAAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTcta                                                                            <  7:377550/149‑1 (MQ=255)
                                                                                      tttGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGCTGAAAATTTGACTAGTAAGAATTCTACAATTAAGAATTCAAGGCGACgg                                                     >  1:120079/1‑148 (MQ=37)
                                                                                        tGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAATTTTGCAGTTA‑AGGTGTGTCAAAAGTTGCCGTTTAGGTGGGTGTAAAGTTGCCGAGTAAGGATTCTACTATTAATAATTCTAGGCGACGGGat                                                  >  1:344178/1‑148 (MQ=37)
                                                                                         gCCTTAATTCCATGAATTCAGGTGAGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGt                                                 >  7:227476/1‑149 (MQ=255)
                                                                                               aTTCCATGAATTCAGCTGGGTAAAAATGTGCCTTTATTGGTGGG‑TCAAAAGTTGCAGTTA‑ATGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                           <  5:176133/149‑1 (MQ=37)
                                                                                                   cATGAATTCAGGGGGCTCAAAATAAGGCCTTAGTGGGGGG‑TCAAAAGATGCCGTTA‑CGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCt                                       <  1:195879/149‑1 (MQ=25)
                                                                                                    aTGAATTCAGTTTGGTAAAAATAAGCCTTTATGGGTGGG‑TCATATCTTGCCGTTT‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTc                                      <  8:353517/149‑1 (MQ=25)
                                                                                                             ggTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGCGGAGTAAGAATTGTACTATTAAGAATTGTAGGCGACGGGAGTGTTTTAGCGTCCcggtcgg                              >  1:124521/1‑148 (MQ=37)
                                                                                                                   ctaaaaGATGCCGTTATTGTTGCC‑TAAAAAGCAGCCGTTA‑ATGTGGGGCAAAAGTAGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATaga                       <  1:86000/147‑1 (MQ=255)
                                                                                                                   ccAAAAGGAGCCTTTAGTGGTGGG‑TAAAAAGTTGCCTTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATaga                       <  5:385477/148‑1 (MQ=38)
                                                                                                                      cacgTTGCCGTTAGTGGGGGG‑TCAAAACTTGCCGTCA‑CGGGGGGTCAAAAGTTGCCTTTAAGGTGCGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagaga                     <  7:340837/145‑1 (MQ=255)
                                                                                                                      aaaGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagagag                    >  1:316488/1‑149 (MQ=255)
                                                                                                                        aGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGGCGTGAAGGGGGGTCAAAAGGTGCCGAGATAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTc                  >  2:371509/1‑149 (MQ=38)
                                                                                                                           aggccttACTGGGGGG‑TAAAAACAGGCCTTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGcc               <  3:368017/144‑1 (MQ=38)
                                                                                                                             gcGTTACGGGTGGG‑TCACAACACGCCGTTA‑CGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTa             <  7:364276/148‑1 (MQ=38)
                                                                                                                             cccttATTGTTGGG‑TAAATACATGCCTTTA‑AGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTa             <  2:356036/146‑1 (MQ=38)
                                                                                                                                tttctggtgcG‑GCAAAAGTGTCCGTTA‑AGTTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTagag          <  7:323/145‑1 (MQ=38)
                                                                                                                                ttAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAAGTGGGTCAAAAGTTGACGAGGAAGAATTCTAATATTAAGAATTCTAGCCGACTGGACTTTATAAGCATCAAGGTCCGTCATATAGAGTCGCATagat          >  6:178068/1‑148 (MQ=21)
                                                                                                                                taaGTGGTGCG‑TCAAAAGTTCACGATA‑ACGTGGGTCAAAAGATGCCGTAAACGAGGCTCAAAAGTTCCCCACTAAGAATACTACAATAAAGAATTCTCCACGACAGCACTTTATTAGCCTCCAGTACAGTCATAGAGAGTCGCCTagag          >  6:183829/3‑149 (MQ=255)
                                                                                                                                cttttggtggGCTCATATGATGCCTTTT‑ATGTGGGTCAAAAGTTGCACTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga           <  2:46073/145‑1 (MQ=38)
                                                                                                                                attcaggtggG‑TCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga           <  7:40183/144‑1 (MQ=255)
                                                                                                                                       tggG‑TCAAAAGTTTCAGTTA‑AGGTGGGTCAATAGTTGCCGTTAAGGTTGGTCAAATGTTGCCGAGTAATTTTTCTACTATTAAGAATTCTAGGCGACGGGAGGTTTTTAGCCTCCCGGTCGGTCAtt                         <  1:389939/127‑2 (MQ=255)
                                                                                                                                       gggg‑TCATAAGATGCCTTTA‑AGGTGGGTCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACa   <  3:282561/148‑1 (MQ=255)
                                                                                                                                            tCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCACGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag  >  2:103357/1‑146 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
GAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAG  >  NC_000913/2560845‑2561129

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