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

GGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑A‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCG  >  NC_000913/2563415‑2563607
                                                                                                   |                                                                                               
ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑a                                                                                                >  1:793624/1‑100 (MQ=255)
  tGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑a‑aa                                                                                             >  2:823231/1‑100 (MQ=255)
      ttttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑a‑aaaaat                                                                                         >  1:1027275/1‑98 (MQ=255)
       tttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑a‑aaaaatg                                                                                        <  2:1129036/100‑4 (MQ=255)
        ttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑a‑aaaaatga                                                                                       <  1:1735060/100‑5 (MQ=255)
        ttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑a‑aaaaatga                                                                                       <  1:1543967/100‑5 (MQ=255)
            gATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACg                                                                                    >  2:1478385/1‑100 (MQ=255)
              tGCAGATCGCTTTTTTTTTTTTTTTATTTTTGTTTAATTCAATCAACGATATGCAGTTTTCATCAACAATATAAAGATGGGAGGT‑A‑AAAATGAAAAGGt                                                                                  >  1:373878/1‑99 (MQ=255)
                 aGATCGCTTTTTCTTTTTAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTggg                                                                               >  1:1469581/1‑100 (MQ=255)
                 aGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑A‑AAAATGAAACGGTGGGa                                                                              >  2:1256180/1‑100 (MQ=255)
                  gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATAAAAAGTAA‑AAAATGAAACGGTGGGa                                                                              <  1:1196841/100‑1 (MQ=255)
                   aTCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGa                                                                              <  1:1050206/99‑1 (MQ=255)
                    tCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑A‑AAAATGAAACGGTGGGAAGt                                                                           <  1:586890/100‑1 (MQ=255)
                      gCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAAGGGTGGGaa                                                                             >  2:1163873/1‑97 (MQ=255)
                       cTTTTTCCCTTCAGGCATTTTTTTTTTTTTTTTTCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGTTCGGa                                                                              >  1:974848/1‑95 (MQ=255)
                       cTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTgg                                                                         <  2:137945/100‑1 (MQ=255)
                       cTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTgg                                                                         <  1:942553/100‑1 (MQ=255)
                        tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGc                                                                        >  2:1473545/1‑100 (MQ=255)
                         ttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑A‑AAAATGAAACGGTGGGAAGTGGCtt                                                                      <  1:809374/100‑1 (MQ=255)
                         ttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑A‑AAAATGAAACGGTGGGAAGTGGCtt                                                                      >  2:1057071/1‑100 (MQ=255)
                           ttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCttt                                                                     <  1:863147/100‑1 (MQ=255)
                           ttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCttt                                                                     <  2:1027276/100‑1 (MQ=255)
                               cTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑A‑AAAATGAAACGGTGGGAAGTGGCTTTACTgg                                                                <  1:1057070/100‑1 (MQ=255)
                                ttCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTgg                                                                <  1:748266/100‑1 (MQ=255)
                                 tCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGt                                                               >  2:214272/1‑100 (MQ=255)
                                 tCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGt                                                               <  1:1669263/100‑1 (MQ=255)
                                  cAGGCATTTTTGTTTAATTCAATCAACGACAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑ATAAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              <  1:1167548/100‑1 (MQ=255)
                                  cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              <  1:1819215/100‑1 (MQ=255)
                                  cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              <  2:27667/100‑1 (MQ=255)
                                  cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              >  2:1409931/1‑100 (MQ=255)
                                   aGGCATTTTTTTTTTTTTTTTTCAACGAAAGGCAGTTCCCAACAACAAATTAAAGATGGGAGGTAA‑AAAATGAAACTGTGGGAAGTGGCTTTACTGGt                                                               >  1:739373/1‑98 (MQ=255)
                                    ggCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGa                                                             >  2:620529/1‑99 (MQ=255)
                                    ggCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGAt                                                            >  2:1598539/1‑100 (MQ=255)
                                      cATTTTTGTTTAAGTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATgg                                                          <  1:414164/100‑1 (MQ=255)
                                        tttttGTTTAATTCAATCAACGAAAGGCAGTTCCCAGCAACAAAATAAAGATGGGAGGT‑A‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTaa                                                       <  2:1364906/100‑1 (MQ=255)
                                         ttttGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTaa                                                       <  1:19617/100‑1 (MQ=255)
                                             gTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                                                   >  2:315367/1‑100 (MQ=255)
                                              tttAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTg                                                  <  1:1555085/100‑1 (MQ=255)
                                               ttAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTgg                                                 <  1:703645/100‑1 (MQ=255)
                                                 aaTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGtt                                               >  2:1550109/1‑100 (MQ=255)
                                                  aTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGttt                                              >  2:409659/1‑100 (MQ=255)
                                                  aTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGttt                                              >  2:1822217/1‑100 (MQ=255)
                                                  aTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGttt                                              >  2:70774/1‑100 (MQ=255)
                                                     caatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCt                                           >  2:1664173/1‑100 (MQ=255)
                                                     caatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑A‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          >  2:1207310/1‑100 (MQ=255)
                                                      aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          >  2:494453/1‑100 (MQ=255)
                                                        tcaaCGAAAGGCAGTACCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTa                                        <  2:1914054/100‑1 (MQ=255)
                                                         caaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt                                       >  2:1017256/1‑100 (MQ=255)
                                                           aCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTg                                     >  1:1661039/1‑100 (MQ=255)
                                                           aCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTg                                     <  2:434465/100‑1 (MQ=255)
                                                            cGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGa                                    <  2:1460538/100‑1 (MQ=255)
                                                              aaaGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAAc                                  <  1:1478377/100‑1 (MQ=255)
                                                               aaGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑A‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACtg                                >  1:190262/1‑100 (MQ=255)
                                                                 ggCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACtgt                               >  1:234773/1‑100 (MQ=255)
                                                                  gCAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACtgtg                              <  1:189513/100‑1 (MQ=255)
                                                                   cAGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGa                             <  1:250700/100‑1 (MQ=255)
                                                                    aGTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGaa                            >  2:158773/1‑100 (MQ=255)
                                                                     gTTCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAag                           <  1:823230/100‑1 (MQ=255)
                                                                      ttCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagg                          <  1:1215070/100‑1 (MQ=255)
                                                                       tCCCAACAACAAAATAAAGATGGGAGGTAA‑AAAACGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagga                         <  1:1030220/100‑1 (MQ=255)
                                                                        cccAACAACAAAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagga                         <  2:793626/99‑1 (MQ=255)
                                                                        cccAACAACAAAATAAAGATGGGAGGT‑A‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggagg                       >  2:1432939/1‑100 (MQ=255)
                                                                           aacaacaaAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGAc                     >  1:1167019/1‑100 (MQ=255)
                                                                            acaacaaAATAAAGATGGGAGGT‑A‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGc                   >  1:869587/1‑100 (MQ=255)
                                                                             caacaaAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGc                   >  1:283068/1‑100 (MQ=255)
                                                                                caaAATAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGTTGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCAt                >  2:787416/1‑100 (MQ=255)
                                                                                  aaaTAAAGATGGGAGGT‑A‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGaa             <  2:883148/100‑1 (MQ=255)
                                                                                   aaTAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGaa             >  2:1795418/1‑100 (MQ=255)
                                                                                     tAAAGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAAc           <  1:259488/100‑1 (MQ=255)
                                                                                       aaGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtt         <  1:44074/100‑1 (MQ=255)
                                                                                        aGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACttt        <  2:449037/100‑1 (MQ=255)
                                                                                        aGATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACttt        <  1:214271/100‑1 (MQ=255)
                                                                                         gATGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt       <  1:1441021/100‑1 (MQ=255)
                                                                                          aTGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTa      <  1:1724253/100‑1 (MQ=255)
                                                                                          aTGGGAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTa      >  2:1663160/1‑100 (MQ=255)
                                                                                            gggAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAAt    <  2:974849/100‑1 (MQ=255)
                                                                                              gAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGATCTGTGAAGGAGGACGCCATGAAACTTTTAATCg  <  1:1167365/100‑1 (MQ=255)
                                                                                              gAGGTAA‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCg  <  2:1600562/100‑1 (MQ=255)
                                                                                                   |                                                                                               
GGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑A‑AAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCG  >  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: ‑