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.A95.6% 234.1 / 3.7 68intergenic (+396/‑79)yffP/yffQCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base A (37/28);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

AGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCG  >  NC_000913/2563414‑2563607
                                                                                                    |                                                                                               
aGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGt                                                                                                  <  1:444630/100‑1 (MQ=255)
  gTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaa                                                                                               >  2:1046749/1‑100 (MQ=255)
  gTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaa                                                                                               >  2:580289/1‑100 (MQ=255)
   tGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaa                                                                                              <  1:1031909/100‑1 (MQ=255)
   tGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaa                                                                                              >  2:1394505/1‑100 (MQ=255)
      aTTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaaa                                                                                           <  2:1509288/100‑2 (MQ=255)
        tttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaaatg                                                                                         <  1:1263650/100‑4 (MQ=255)
        tttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaaatg                                                                                         <  1:580291/100‑4 (MQ=255)
        tttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaaatg                                                                                         <  1:1441234/100‑4 (MQ=255)
           gtgATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGaaa                                                                                      >  2:1925898/1‑100 (MQ=255)
               tGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGt                                                                                  >  1:446310/1‑100 (MQ=255)
               tGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGt                                                                                  >  1:1774258/1‑100 (MQ=255)
                gCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTg                                                                                 <  1:992420/100‑1 (MQ=255)
                  aGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTggg                                                                               <  2:1774244/100‑1 (MQ=255)
                   gATCGCTTTTTCCTTTTAGTCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTTCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGa                                                                              >  2:818172/1‑100 (MQ=255)
                   gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGa                                                                              >  2:547369/1‑100 (MQ=255)
                   gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGa                                                                              >  1:708800/1‑100 (MQ=255)
                   gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGa                                                                              >  2:269873/1‑100 (MQ=255)
                   gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGa                                                                              >  1:1266412/1‑100 (MQ=255)
                    aTCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGaa                                                                             <  1:890780/100‑1 (MQ=255)
                    aTCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGaa                                                                             >  1:867786/1‑100 (MQ=255)
                    aTCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGaa                                                                             >  1:867527/1‑100 (MQ=255)
                       gCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTg                                                                          >  1:1738019/1‑100 (MQ=255)
                        cTTTTTCCCTTCAGGCATTTTTTTTTTTTTTTTTCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGATAAAAAATGAAACGGTGGGAAGt                                                                           >  2:1300693/1‑98 (MQ=255)
                         tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGc                                                                        <  1:1768745/99‑1 (MQ=255)
                            ttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCttt                                                                     <  2:55631/100‑1 (MQ=255)
                            ttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCttt                                                                     >  2:1082062/1‑100 (MQ=255)
                               ccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACt                                                                  <  1:848527/100‑1 (MQ=255)
                                cTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTg                                                                 >  2:1560154/1‑100 (MQ=255)
                                 ttCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGATAAAAAATGAAACGGTGGGAAGTGGCTTTACTgg                                                                <  1:1300700/100‑1 (MQ=255)
                                   cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              <  2:421228/100‑1 (MQ=255)
                                     ggCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGAt                                                            >  1:1940815/1‑100 (MQ=255)
                                     ggCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATg                                                           <  1:315265/100‑1 (MQ=255)
                                       cATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATgg                                                          <  2:910353/100‑1 (MQ=255)
                                        aTTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGt                                                         >  2:2005962/1‑100 (MQ=255)
                                        aTTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGt                                                         >  2:2027757/1‑100 (MQ=255)
                                        aTTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAACTGAAACGGTGGGAAGTGGCTTTACTGGTGATGGt                                                         >  2:728084/1‑100 (MQ=255)
                                         tttttGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTa                                                        >  1:1876064/1‑100 (MQ=255)
                                          ttttGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTaa                                                       >  2:897920/1‑100 (MQ=255)
                                          ttttGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTaa                                                       <  1:1649405/100‑1 (MQ=255)
                                               tttAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTg                                                  >  1:1724178/1‑100 (MQ=255)
                                                ttAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTgg                                                 <  1:2042573/100‑1 (MQ=255)
                                                ttAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTgg                                                 <  2:1648770/100‑1 (MQ=255)
                                                 tAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGt                                                >  2:434556/1‑100 (MQ=255)
                                                 tAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGt                                                >  1:661894/1‑100 (MQ=255)
                                                   aTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGttt                                              >  2:728133/1‑100 (MQ=255)
                                                   aTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGttt                                              >  2:1982907/1‑100 (MQ=255)
                                                     tcaatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGc                                            <  1:434557/100‑1 (MQ=255)
                                                     tcaatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGc                                            <  2:1538095/100‑1 (MQ=255)
                                                       aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          <  1:1330665/100‑1 (MQ=255)
                                                       aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          <  2:1947302/100‑1 (MQ=255)
                                                        atcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGt                                         >  1:1716166/1‑100 (MQ=255)
                                                          caaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt                                       <  1:1042929/100‑1 (MQ=255)
                                                          caaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt                                       >  2:212560/1‑100 (MQ=255)
                                                          caaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt                                       <  1:1394514/100‑1 (MQ=255)
                                                               aaaGGCAGTTAACAACAATAAAAGAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAAc                                  <  2:13343/100‑1 (MQ=255)
                                                                  ggCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACtgt                               >  1:1173548/1‑100 (MQ=255)
                                                                       cttccgatctcaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagg                          <  2:14558/90‑1 (MQ=255)
                                                                        tCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagga                         <  1:693396/100‑1 (MQ=255)
                                                                         cccAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggag                        >  2:848027/1‑100 (MQ=255)
                                                                                acaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa                 >  2:438024/1‑100 (MQ=255)
                                                                                acaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa                 >  1:1820726/1‑100 (MQ=255)
                                                                                acaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa                 >  1:1281726/1‑100 (MQ=255)
                                                                                 caaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggagatcggaag                >  1:14558/1‑92 (MQ=255)
                                                                                 caaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCAt                <  2:986036/100‑1 (MQ=255)
                                                                                 caaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCAt                <  1:547370/100‑1 (MQ=255)
                                                                                    aaTAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGaa             <  1:736098/100‑1 (MQ=255)
                                                                                      tAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAAc           <  2:1851098/100‑1 (MQ=255)
                                                                                       aaaGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACt          >  1:1807999/1‑100 (MQ=255)
                                                                                        aaGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtt         >  1:1456260/1‑100 (MQ=255)
                                                                                          gATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt       <  1:725445/100‑1 (MQ=255)
                                                                                          gATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt       <  1:269874/100‑1 (MQ=255)
                                                                                          gATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt       >  1:1523728/1‑100 (MQ=255)
                                                                                            tGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTaa     >  1:1461781/1‑100 (MQ=255)
                                                                                              ggAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAAt    >  1:1836944/1‑99 (MQ=255)
                                                                                               gAGGTAAAAAATGAATCGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCg  <  1:2009191/100‑1 (MQ=255)
                                                                                               gAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTaa     <  2:872019/97‑1 (MQ=255)
                                                                                                    |                                                                                               
AGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCG  >  NC_000913/2563414‑2563607

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