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

AGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑AAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCG  >  NC_000913/2563414‑2563607
                                                                                                    |                                                                                               
aGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGt                                                                                                  <  1:444630/100‑1 (MQ=255)
  gTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑aa                                                                                               >  2:1046749/1‑100 (MQ=255)
  gTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑aa                                                                                               >  2:580289/1‑100 (MQ=255)
   tGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑aaa                                                                                              <  1:1031909/100‑1 (MQ=255)
   tGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑aaa                                                                                              >  2:1394505/1‑100 (MQ=255)
      aTTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑aaaaaa                                                                                           <  2:1509288/100‑2 (MQ=255)
        tttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑aaaaaatg                                                                                         <  1:1263650/100‑4 (MQ=255)
        tttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑aaaaaatg                                                                                         <  1:580291/100‑4 (MQ=255)
        tttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑aaaaaatg                                                                                         <  1:1441234/100‑4 (MQ=255)
           gtgATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGaaa                                                                                      >  2:1925898/1‑100 (MQ=255)
               tGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGt                                                                                  >  1:446310/1‑100 (MQ=255)
               tGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGt                                                                                  >  1:1774258/1‑100 (MQ=255)
                gCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTg                                                                                 <  1:992420/100‑1 (MQ=255)
                  aGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTggg                                                                               <  2:1774244/100‑1 (MQ=255)
                   gATCGCTTTTTCCTTTTAGTCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTTCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGa                                                                              >  2:818172/1‑100 (MQ=255)
                   gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGa                                                                              >  2:547369/1‑100 (MQ=255)
                   gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGa                                                                              >  1:708800/1‑100 (MQ=255)
                   gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGa                                                                              >  2:269873/1‑100 (MQ=255)
                   gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGa                                                                              >  1:1266412/1‑100 (MQ=255)
                    aTCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGaa                                                                             <  1:890780/100‑1 (MQ=255)
                    aTCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGaa                                                                             >  1:867786/1‑100 (MQ=255)
                    aTCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGaa                                                                             >  1:867527/1‑100 (MQ=255)
                       gCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTg                                                                          >  1:1738019/1‑100 (MQ=255)
                        cTTTTTCCCTTCAGGCATTTTTTTTTTTTTTTTTCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGATAAAAAA‑TGAAACGGTGGGAAGt                                                                           >  2:1300693/1‑98 (MQ=255)
                         tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑AAAAA‑TGAAACGGTGGGAAGTGGc                                                                        <  1:1768745/99‑1 (MQ=255)
                            ttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCttt                                                                     <  2:55631/100‑1 (MQ=255)
                            ttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCttt                                                                     >  2:1082062/1‑100 (MQ=255)
                               ccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACt                                                                  <  1:848527/100‑1 (MQ=255)
                                cTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTg                                                                 >  2:1560154/1‑100 (MQ=255)
                                 ttCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGATAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTgg                                                                <  1:1300700/100‑1 (MQ=255)
                                   cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              <  2:421228/100‑1 (MQ=255)
                                     ggCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGAt                                                            >  1:1940815/1‑100 (MQ=255)
                                     ggCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑AAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATg                                                           <  1:315265/100‑1 (MQ=255)
                                       cATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATgg                                                          <  2:910353/100‑1 (MQ=255)
                                        aTTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGt                                                         >  2:2005962/1‑100 (MQ=255)
                                        aTTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGt                                                         >  2:2027757/1‑100 (MQ=255)
                                        aTTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑AAAAACTGAAACGGTGGGAAGTGGCTTTACTGGTGATGGt                                                         >  2:728084/1‑100 (MQ=255)
                                         tttttGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTa                                                        >  1:1876064/1‑100 (MQ=255)
                                          ttttGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTaa                                                       >  2:897920/1‑100 (MQ=255)
                                          ttttGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTaa                                                       <  1:1649405/100‑1 (MQ=255)
                                               tttAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTg                                                  >  1:1724178/1‑100 (MQ=255)
                                                ttAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTgg                                                 <  1:2042573/100‑1 (MQ=255)
                                                ttAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTgg                                                 <  2:1648770/100‑1 (MQ=255)
                                                 tAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGt                                                >  2:434556/1‑100 (MQ=255)
                                                 tAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGt                                                >  1:661894/1‑100 (MQ=255)
                                                   aTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGttt                                              >  2:728133/1‑100 (MQ=255)
                                                   aTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGttt                                              >  2:1982907/1‑100 (MQ=255)
                                                     tcaatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGc                                            <  1:434557/100‑1 (MQ=255)
                                                     tcaatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGc                                            <  2:1538095/100‑1 (MQ=255)
                                                       aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          <  1:1330665/100‑1 (MQ=255)
                                                       aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          <  2:1947302/100‑1 (MQ=255)
                                                        atcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGt                                         >  1:1716166/1‑100 (MQ=255)
                                                          caaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt                                       <  1:1042929/100‑1 (MQ=255)
                                                          caaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt                                       >  2:212560/1‑100 (MQ=255)
                                                          caaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt                                       <  1:1394514/100‑1 (MQ=255)
                                                               aaaGGCAGTTAACAACAATAAAAGAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAAc                                  <  2:13343/100‑1 (MQ=255)
                                                                  ggCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACtgt                               >  1:1173548/1‑100 (MQ=255)
                                                                       cttccgatctcaaAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagg                          <  2:14558/90‑1 (MQ=255)
                                                                        tCCCAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagga                         <  1:693396/100‑1 (MQ=255)
                                                                         cccAACAACAAAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggag                        >  2:848027/1‑100 (MQ=255)
                                                                                acaaAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa                 >  2:438024/1‑100 (MQ=255)
                                                                                acaaAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa                 >  1:1820726/1‑100 (MQ=255)
                                                                                acaaAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa                 >  1:1281726/1‑100 (MQ=255)
                                                                                 caaAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggagatcggaag                >  1:14558/1‑92 (MQ=255)
                                                                                 caaAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCAt                <  2:986036/100‑1 (MQ=255)
                                                                                 caaAATAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCAt                <  1:547370/100‑1 (MQ=255)
                                                                                    aaTAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGaa             <  1:736098/100‑1 (MQ=255)
                                                                                      tAAAGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAAc           <  2:1851098/100‑1 (MQ=255)
                                                                                       aaaGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACt          >  1:1807999/1‑100 (MQ=255)
                                                                                        aaGATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtt         >  1:1456260/1‑100 (MQ=255)
                                                                                          gATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt       <  1:725445/100‑1 (MQ=255)
                                                                                          gATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt       <  1:269874/100‑1 (MQ=255)
                                                                                          gATGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt       >  1:1523728/1‑100 (MQ=255)
                                                                                            tGGGAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTaa     >  1:1461781/1‑100 (MQ=255)
                                                                                              ggAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAAt    >  1:1836944/1‑99 (MQ=255)
                                                                                               gAGGTAAAAAA‑TGAATCGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCg  <  1:2009191/100‑1 (MQ=255)
                                                                                               gAGGTAAAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTaa     <  2:872019/97‑1 (MQ=255)
                                                                                                    |                                                                                               
AGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑AAAAA‑TGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCG  >  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: ‑