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

GAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAG  >  NC_000913/2560845‑2561129
                                                                                                                                                   |                                                                                                                                           
gAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTCAAAATTTg                                                                                                                                           >  2:195879/1‑149 (MQ=255)
gAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGGGGGTCAAAAGTTg                                                                                                                                           <  6:206373/149‑1 (MQ=255)
   aGGGAGTGCTACGAATACCTGTTATTATTCAATCCGTACAACATGAATTAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCGTTGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGcc                                                                                                                                         <  5:85351/148‑1 (MQ=255)
     ggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGTTGGGTCAAAAGTTGCAGt                                                                                                                                       >  5:87320/1‑148 (MQ=255)
         tGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAgg                                                                                                                                 >  6:362613/1‑149 (MQ=255)
           cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                                                <  5:362613/148‑1 (MQ=255)
           cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGGCAAAAGTTGCCGTTAAGGt                                                                                                                                <  1:399641/148‑1 (MQ=255)
            aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                              >  1:356036/1‑149 (MQ=255)
              cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                                            <  8:18518/149‑1 (MQ=255)
              cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                                            <  7:147270/149‑1 (MQ=255)
                   aCCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGTGCTATGTTCGCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATTGTGATGGGGCAAAAGATGGCGTTAAGGGGTGTGaaat                                                                                                                       >  4:282561/1‑148 (MQ=255)
                   aCCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTACATGGATTCAGGTGTGTCAAAAGTTGACGTTATTTCTGTGTAAAAAGTTGAAGTTAAGgtgtgtaatat                                                                                                                       >  2:86000/1‑144 (MQ=255)
                     cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaatt                                                                                                                     >  5:288752/1‑147 (MQ=255)
                        ttattaTTCAAGCCGTACAACATGAAATAACGGAATCAGGCGTGCAATGTTCTAGCCTTTTCTTTTCCTTAATTACATGAATTCAGGTGGGTCAAAAGTTGACGTTAGTGGTGTTACAAAAGTTGCCGTCAAGGTGGGTCTACAGTTg                                                                                                                   >  7:350182/1‑148 (MQ=255)
                        ttattaTTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGc                                                                                                                  <  7:11988/149‑1 (MQ=255)
                        ttattaTTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTAAGGTGGGTCAAAAGTTGCCGTTTATGGTTTGTAATATTTTGACTTTAAGGTGGGTGaatttttt                                                                                                                   >  3:308487/1‑142 (MQ=255)
                        ttattaTTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTAAGGTGGGTCAAAAGTTGCCGTTAGTGATGTGTAAAAATTTGCAGTTAAGGAGTGTCAAGAGTTg                                                                                                                   >  5:336177/1‑148 (MQ=255)
                          attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                       <  1:11061/142‑1 (MQ=255)
                          attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAACAGTTGcgg                                                                                                                >  5:211267/1‑147 (MQ=255)
                          attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTACCCTTAATTCCATGAATTCAGGTGGCTCAAAAGTTTCCGTAAGTTGTGGGTCAAAAGATGCCGTTAAGGTGGGTCaaaa                                                                                                                       >  2:11061/1‑142 (MQ=255)
                           ttattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGCCTTAAAGGTGGGTCAAAATTTg                                                                                                                   >  1:205995/1‑145 (MQ=255)
                           ttattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCATAAATTACATGAATTCAGGTCGGTCAAAAGTTGCCGTTAGTGGTGTGTCAAATTATGCCGTAATGGGGGGTCAATAGTTgactg                                                                                                               >  5:334095/1‑145 (MQ=255)
                                 aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAATTTGCAGTTAAGGTGGGTCAAAAGTTGCAGTTTAggg                                                                                                         >  6:107979/1‑148 (MQ=255)
                                 aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTGAGTGGTGGGTCAAAAGTTGCCGTTTAGGTGGGTCAAAAGTTGCCgagtaggg                                                                                                         >  5:327987/1‑142 (MQ=255)
                                   gcagTACACCATGAAATAAAGGAATCAGTAGTGCAATTTTCACGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                       <  4:328964/146‑1 (MQ=37)
                                                       ggAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCCGTTAAGGTGGGTCAAAAGTTGCAGTTAAGGTGGGTCAAAAGTTGCCGAGTAAg                                                                                   >  8:45963/1‑149 (MQ=255)
                                                              ggcggcCACTTTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGTGGCGTCAAAAGTTGCCGTTAGTGGTGGGGCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTcta                                                                            <  5:108021/149‑1 (MQ=37)
                                                              ggCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCGGTGAAGGTGGCGCAAAATTTGCCTTTAATGTGGGGGAAAATTTTCCTAGTTAGAATCctc                                                                            >  3:23727/1‑148 (MQ=255)
                                                              ggCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTACATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGGGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTcta                                                                            <  7:377550/149‑1 (MQ=255)
                                                                                      tttGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGCTGAAAATTTGACTAGTAAGAATTCTACAATTAAGAATTCAAGGCGACgg                                                     >  1:120079/1‑148 (MQ=37)
                                                                                        tGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATTTTGCAGTTAAGGTGTGTCAAAAGTTGCCGTTTAGGTGGGTGTAAAGTTGCCGAGTAAGGATTCTACTATTAATAATTCTAGGCGACGGGat                                                  >  1:344178/1‑148 (MQ=37)
                                                                                         gCCTTAATTCCATGAATTCAGGTGAGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGt                                                 >  7:227476/1‑149 (MQ=255)
                                                                                               aTTCCATGAATTCAGCTGGGTAAAAATGTGCCTTTATTGGTGGGTCAAAAGTTGCAGTTAATGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                           <  5:176133/149‑1 (MQ=37)
                                                                                                   cATGAATTCAGGGGGCTCAAAATAAGGCCTTAGTGGGGGGTCAAAAGATGCCGTTACGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCt                                       <  1:195879/149‑1 (MQ=25)
                                                                                                    aTGAATTCAGTTTGGTAAAAATAAGCCTTTATGGGTGGGTCATATCTTGCCGTTTAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTc                                      <  8:353517/149‑1 (MQ=25)
                                                                                                             ggTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGCGGAGTAAGAATTGTACTATTAAGAATTGTAGGCGACGGGAGTGTTTTAGCGTCCcggtcgg                              >  1:124521/1‑148 (MQ=37)
                                                                                                                   ctaaaaGATGCCGTTATTGTTGCCTAAAAAGCAGCCGTTAATGTGGGGCAAAAGTAGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATaga                       <  1:86000/147‑1 (MQ=255)
                                                                                                                   ccAAAAGGAGCCTTTAGTGGTGGGTAAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATaga                       <  5:385477/148‑1 (MQ=38)
                                                                                                                      cacgTTGCCGTTAGTGGGGGGTCAAAACTTGCCGTCACGGGGGGTCAAAAGTTGCCTTTAAGGTGCGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagaga                     <  7:340837/145‑1 (MQ=255)
                                                                                                                      aaaGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagagag                    >  1:316488/1‑149 (MQ=255)
                                                                                                                        aGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGGCGTGAAGGGGGGTCAAAAGGTGCCGAGATAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTc                  >  2:371509/1‑149 (MQ=38)
                                                                                                                           aggccttACTGGGGGGTAAAAACAGGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGcc               <  3:368017/144‑1 (MQ=38)
                                                                                                                             gcGTTACGGGTGGGTCACAACACGCCGTTACGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTa             <  7:364276/148‑1 (MQ=38)
                                                                                                                             cccttATTGTTGGGTAAATACATGCCTTTAAGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTa             <  2:356036/146‑1 (MQ=38)
                                                                                                                                tttctggtgcGGCAAAAGTGTCCGTTAAGTTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTagag          <  7:323/145‑1 (MQ=38)
                                                                                                                                ttAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAAGTGGGTCAAAAGTTGACGAGGAAGAATTCTAATATTAAGAATTCTAGCCGACTGGACTTTATAAGCATCAAGGTCCGTCATATAGAGTCGCATagat          >  6:178068/1‑148 (MQ=21)
                                                                                                                                taaGTGGTGCGTCAAAAGTTCACGATAACGTGGGTCAAAAGATGCCGTAAACGAGGCTCAAAAGTTCCCCACTAAGAATACTACAATAAAGAATTCTCCACGACAGCACTTTATTAGCCTCCAGTACAGTCATAGAGAGTCGCCTagag          >  6:183829/3‑149 (MQ=255)
                                                                                                                                cttttggtggGCTCATATGATGCCTTTTATGTGGGTCAAAAGTTGCACTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga           <  2:46073/145‑1 (MQ=38)
                                                                                                                                attcaggtggGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga           <  7:40183/144‑1 (MQ=255)
                                                                                                                                       tggGTCAAAAGTTTCAGTTAAGGTGGGTCAATAGTTGCCGTTAAGGTTGGTCAAATGTTGCCGAGTAATTTTTCTACTATTAAGAATTCTAGGCGACGGGAGGTTTTTAGCCTCCCGGTCGGTCAtt                         <  1:389939/127‑2 (MQ=255)
                                                                                                                                       ggggTCATAAGATGCCTTTAAGGTGGGTCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACa   <  3:282561/148‑1 (MQ=255)
                                                                                                                                            tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCACGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag  >  2:103357/1‑146 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
GAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAG  >  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: