Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,563,518 (A)5→6 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.A81.9% 214.4 / 42.9 72intergenic (+396/‑79)yffP/yffQCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (7/6);  new base A (26/33);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

GGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCG  >  NC_000913/2563415‑2563607
                                                                                                   |                                                                                               
ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTa                                                                                                >  1:793624/1‑100 (MQ=255)
  tGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaa                                                                                             >  2:823231/1‑100 (MQ=255)
      ttttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaaat                                                                                         >  1:1027275/1‑98 (MQ=255)
       tttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaaatg                                                                                        <  2:1129036/100‑4 (MQ=255)
        ttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaaatga                                                                                       <  1:1735060/100‑5 (MQ=255)
        ttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaaatga                                                                                       <  1:1543967/100‑5 (MQ=255)
            gATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACg                                                                                    >  2:1478385/1‑100 (MQ=255)
              tGCAGATCGCTTTTTTTTTTTTTTTATTTTTGTTTAATTCAATCAACGATATGCAGTTTTCATCAACAATATAAAGATGGGAGGTAAAAATGAAAAGGt                                                                                  >  1:373878/1‑99 (MQ=255)
                 aGATCGCTTTTTCTTTTTAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTggg                                                                               >  1:1469581/1‑100 (MQ=255)
                 aGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGa                                                                              >  2:1256180/1‑100 (MQ=255)
                  gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATAAAAAGTAAAAAATGAAACGGTGGGa                                                                              <  1:1196841/100‑1 (MQ=255)
                   aTCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGa                                                                              <  1:1050206/99‑1 (MQ=255)
                    tCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGt                                                                           <  1:586890/100‑1 (MQ=255)
                      gCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAAGGGTGGGaa                                                                             >  2:1163873/1‑97 (MQ=255)
                       cTTTTTCCCTTCAGGCATTTTTTTTTTTTTTTTTCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGTTCGGa                                                                              >  1:974848/1‑95 (MQ=255)
                       cTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTgg                                                                         <  2:137945/100‑1 (MQ=255)
                       cTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTgg                                                                         <  1:942553/100‑1 (MQ=255)
                        tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGc                                                                        >  2:1473545/1‑100 (MQ=255)
                         ttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCtt                                                                      <  1:809374/100‑1 (MQ=255)
                         ttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCtt                                                                      >  2:1057071/1‑100 (MQ=255)
                           ttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCttt                                                                     <  1:863147/100‑1 (MQ=255)
                           ttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCttt                                                                     <  2:1027276/100‑1 (MQ=255)
                               cTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTgg                                                                <  1:1057070/100‑1 (MQ=255)
                                ttCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTgg                                                                <  1:748266/100‑1 (MQ=255)
                                 tCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGt                                                               >  2:214272/1‑100 (MQ=255)
                                 tCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGt                                                               <  1:1669263/100‑1 (MQ=255)
                                  cAGGCATTTTTGTTTAATTCAATCAACGACAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTATAAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              <  1:1167548/100‑1 (MQ=255)
                                  cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              <  1:1819215/100‑1 (MQ=255)
                                  cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              <  2:27667/100‑1 (MQ=255)
                                  cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              >  2:1409931/1‑100 (MQ=255)
                                   aGGCATTTTTTTTTTTTTTTTTCAACGAAAGGCAGTTCCCAACAACAAATTAAAGATGGGAGGTAAAAAATGAAACTGTGGGAAGTGGCTTTACTGGt                                                               >  1:739373/1‑98 (MQ=255)
                                    ggCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGa                                                             >  2:620529/1‑99 (MQ=255)
                                    ggCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGAt                                                            >  2:1598539/1‑100 (MQ=255)
                                      cATTTTTGTTTAAGTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATgg                                                          <  1:414164/100‑1 (MQ=255)
                                        tttttGTTTAATTCAATCAACGAAAGGCAGTTCCCAGCAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTaa                                                       <  2:1364906/100‑1 (MQ=255)
                                         ttttGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTaa                                                       <  1:19617/100‑1 (MQ=255)
                                             gTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                                                   >  2:315367/1‑100 (MQ=255)
                                              tttAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTg                                                  <  1:1555085/100‑1 (MQ=255)
                                               ttAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTgg                                                 <  1:703645/100‑1 (MQ=255)
                                                 aaTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGtt                                               >  2:1550109/1‑100 (MQ=255)
                                                  aTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGttt                                              >  2:409659/1‑100 (MQ=255)
                                                  aTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGttt                                              >  2:1822217/1‑100 (MQ=255)
                                                  aTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGttt                                              >  2:70774/1‑100 (MQ=255)
                                                     caatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCt                                           >  2:1664173/1‑100 (MQ=255)
                                                     caatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          >  2:1207310/1‑100 (MQ=255)
                                                      aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          >  2:494453/1‑100 (MQ=255)
                                                        tcaaCGAAAGGCAGTACCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTa                                        <  2:1914054/100‑1 (MQ=255)
                                                         caaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt                                       >  2:1017256/1‑100 (MQ=255)
                                                           aCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTg                                     >  1:1661039/1‑100 (MQ=255)
                                                           aCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTg                                     <  2:434465/100‑1 (MQ=255)
                                                            cGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGa                                    <  2:1460538/100‑1 (MQ=255)
                                                              aaaGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAAc                                  <  1:1478377/100‑1 (MQ=255)
                                                               aaGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACtg                                >  1:190262/1‑100 (MQ=255)
                                                                 ggCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACtgt                               >  1:234773/1‑100 (MQ=255)
                                                                  gCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACtgtg                              <  1:189513/100‑1 (MQ=255)
                                                                   cAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGa                             <  1:250700/100‑1 (MQ=255)
                                                                    aGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGaa                            >  2:158773/1‑100 (MQ=255)
                                                                     gTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAag                           <  1:823230/100‑1 (MQ=255)
                                                                      ttCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagg                          <  1:1215070/100‑1 (MQ=255)
                                                                       tCCCAACAACAAAATAAAGATGGGAGGTAAAAAACGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagga                         <  1:1030220/100‑1 (MQ=255)
                                                                        cccAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagga                         <  2:793626/99‑1 (MQ=255)
                                                                        cccAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggagg                       >  2:1432939/1‑100 (MQ=255)
                                                                           aacaacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGAc                     >  1:1167019/1‑100 (MQ=255)
                                                                            acaacaaAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGc                   >  1:869587/1‑100 (MQ=255)
                                                                             caacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGc                   >  1:283068/1‑100 (MQ=255)
                                                                                caaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGTTGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCAt                >  2:787416/1‑100 (MQ=255)
                                                                                  aaaTAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGaa             <  2:883148/100‑1 (MQ=255)
                                                                                   aaTAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGaa             >  2:1795418/1‑100 (MQ=255)
                                                                                     tAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAAc           <  1:259488/100‑1 (MQ=255)
                                                                                       aaGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtt         <  1:44074/100‑1 (MQ=255)
                                                                                        aGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACttt        <  2:449037/100‑1 (MQ=255)
                                                                                        aGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACttt        <  1:214271/100‑1 (MQ=255)
                                                                                         gATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt       <  1:1441021/100‑1 (MQ=255)
                                                                                          aTGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTa      <  1:1724253/100‑1 (MQ=255)
                                                                                          aTGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTa      >  2:1663160/1‑100 (MQ=255)
                                                                                            gggAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAAt    <  2:974849/100‑1 (MQ=255)
                                                                                              gAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGATCTGTGAAGGAGGACGCCATGAAACTTTTAATCg  <  1:1167365/100‑1 (MQ=255)
                                                                                              gAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCg  <  2:1600562/100‑1 (MQ=255)
                                                                                                   |                                                                                               
GGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCG  >  NC_000913/2563415‑2563607

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