Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,563,518:1 (A)5→6 100% intergenic (+401/‑74) yffP → / → yffQ 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,563,5131.A98.5% 228.4 / ‑3.4 67intergenic (+396/‑79)yffP/yffQCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (36/30);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTA  >  NC_000913/2563415‑2563603
                                                                                                   |                                                                                          
ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTa                                                                                           >  1:366553/1‑100 (MQ=255)
ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTa                                                                                           >  1:1927382/1‑100 (MQ=255)
  tGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaa                                                                                         >  1:1718386/1‑100 (MQ=255)
  ngAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaa                                                                                         >  1:72072/2‑100 (MQ=255)
    aaTTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaa                                                                                       >  1:480374/1‑100 (MQ=255)
       tttGTGATGCAGATCGCGTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaaatg                                                                                    <  2:1507911/100‑4 (MQ=255)
        ttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaaatga                                                                                   >  2:742154/1‑96 (MQ=255)
          gtgATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGaaa                                                                                  >  2:124147/1‑100 (MQ=255)
           tgATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAAc                                                                                 <  1:2443156/100‑1 (MQ=255)
            gATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACg                                                                                >  1:1688898/1‑100 (MQ=255)
            gATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACg                                                                                >  1:2197857/1‑100 (MQ=255)
            gATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACg                                                                                >  1:2226491/1‑100 (MQ=255)
              tGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACTAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACgg                                                                               >  2:1148199/1‑99 (MQ=255)
              tGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGACAGGCAGTTCCTAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACgg                                                                               >  2:2495488/1‑99 (MQ=255)
              tGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGt                                                                              >  2:30380/1‑100 (MQ=255)
                cAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTgg                                                                            >  1:1484546/1‑100 (MQ=255)
                  gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGa                                                                          >  1:1625489/1‑100 (MQ=255)
                  gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGa                                                                          >  1:1989476/1‑100 (MQ=255)
                  gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGa                                                                          <  2:521549/100‑1 (MQ=255)
                  gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGa                                                                          >  1:340429/1‑100 (MQ=255)
                    tCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAg                                                                        >  1:943229/1‑100 (MQ=255)
                        tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGc                                                                    <  2:101653/100‑1 (MQ=255)
                           ttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCttt                                                                 >  2:955109/1‑100 (MQ=255)
                            tCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTa                                                                <  2:1939110/100‑1 (MQ=255)
                             cccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTAc                                                               <  2:2070087/100‑1 (MQ=255)
                             cccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTAc                                                               >  2:2285799/1‑100 (MQ=255)
                              ccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCACAACAACAAAATAAAGATAGAAAGTAAAAAATGAAACGGTGGGAAGTGGCTTTACt                                                              <  2:683429/100‑1 (MQ=255)
                               cTTCAGGCATTTTTGTTTAATTCAATGAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTg                                                             <  2:505899/100‑1 (MQ=255)
                                  cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                          <  2:2418631/100‑1 (MQ=255)
                                  cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                          <  1:858460/100‑1 (MQ=255)
                                  cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGAAAAAAAAAAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGt                                                           <  2:693716/99‑1 (MQ=255)
                                   aGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGa                                                         >  1:1406912/1‑100 (MQ=255)
                                       aTTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGt                                                     <  2:1059207/100‑1 (MQ=255)
                                             gTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                                               >  1:2154307/1‑100 (MQ=255)
                                             gTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                                               >  2:2191391/1‑100 (MQ=255)
                                              tttAATAAAATCAACCAAAGGCAGTTCCCATCAACGGAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTCGTGaaagcattg                                                 <  2:894339/97‑9 (MQ=255)
                                               ttAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAGTAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTgg                                             >  2:1368050/1‑100 (MQ=255)
                                               ttAACTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTgg                                             <  1:2398906/100‑1 (MQ=255)
                                                tAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTgg                                             >  2:1512385/1‑99 (MQ=255)
                                                tAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGt                                            >  1:1223526/1‑100 (MQ=255)
                                                tAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGt                                            <  1:2125790/100‑1 (MQ=255)
                                                  aTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGttt                                          >  1:1121595/1‑100 (MQ=255)
                                                   ttCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTg                                         >  2:2396462/1‑100 (MQ=255)
                                                      aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGGGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                      <  1:278801/100‑1 (MQ=255)
                                                            cGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGa                                <  2:130624/100‑1 (MQ=255)
                                                               aaGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTNTGCTGTATTGAACt                             <  1:349128/100‑1 (MQ=255)
                                                                  gCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACtgtg                          >  1:1497336/1‑100 (MQ=255)
                                                                    aGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGaa                        >  2:66578/1‑100 (MQ=255)
                                                                    aGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGGGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACtgtg                          <  2:545564/98‑1 (MQ=255)
                                                                     gTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAag                       <  2:1917960/100‑1 (MQ=255)
                                                                     gTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTGACTGGTGATGGTAATGTTGGTTTGCTGGATTGAACTGTGAag                       >  2:468040/1‑100 (MQ=255)
                                                                        cccAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggag                    <  2:1650315/100‑1 (MQ=255)
                                                                         ccAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggagg                   >  1:432529/1‑100 (MQ=255)
                                                                          caacaacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggagg                   >  2:858504/1‑99 (MQ=255)
                                                                          caacaacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAggagga                  <  1:1805061/100‑1 (MQ=255)
                                                                          caacaacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAggagga                  >  1:1216016/1‑100 (MQ=255)
                                                                             caacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGc               <  2:943221/100‑1 (MQ=255)
                                                                              aacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGGTGGTAATGTTGGTTTGCTGTAGTGAAGTGTGAAGGAGGACGcc              >  1:573310/1‑100 (MQ=255)
                                                                              aacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGcc              >  1:962671/1‑100 (MQ=255)
                                                                              aacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGcc              <  1:720793/100‑1 (MQ=255)
                                                                               acaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa             >  2:2198514/1‑100 (MQ=255)
                                                                               acaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa             >  1:603173/1‑100 (MQ=255)
                                                                                caaAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATAGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATg           >  2:2242179/1‑100 (MQ=255)
                                                                                 aaaaTAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATg           <  2:2154293/100‑1 (MQ=255)
                                                                                    aTAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGaaa        <  2:2168191/100‑1 (MQ=255)
                                                                                     tAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAAc       >  2:1820283/1‑100 (MQ=255)
                                                                                      aaaGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAGACt      >  2:136133/1‑100 (MQ=255)
                                                                                      aaaGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACt      <  1:66579/100‑1 (MQ=255)
                                                                                        aGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACttt    <  2:2514596/100‑1 (MQ=255)
                                                                                        aGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACttt    >  1:736822/1‑100 (MQ=255)
                                                                                         gATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt   <  1:1178054/100‑1 (MQ=255)
                                                                                         gATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt   >  1:372836/1‑100 (MQ=255)
                                                                                         gATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt   <  1:2056406/100‑1 (MQ=255)
                                                                                          aTGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTa  <  1:2289241/100‑1 (MQ=255)
                                                                                                   |                                                                                          
GGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTA  >  NC_000913/2563415‑2563603

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: