Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,473,413 C→T 7.0% pseudogene (173/708 nt) oweS ← pseudogene, CPS‑53 (KpLE1) prophage; bacteriophage replication protein O family;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,473,4130CT 7.0% 177.3 / 3.9 57pseudogene (173/708 nt)oweSpseudogene, CPS‑53 (KpLE1) prophage; bacteriophage replication protein O family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (27/26);  new base T (2/2);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

TGGCGCTTTTTAGGGCTGAATACAACCGCATCAGGATGGCGAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCAGAGTCTGACAGTTGAGAAAATGCCGCTGCCTGAAGCTTCGCAACGTTAAGCTG  >  NC_000913/2473266‑2473561
                                                                                                                                                   |                                                                                                                                                    
tGGCGCTTTTTAGGGCTGAATACAACCGCATCAGGATGGCGAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGcc                                                                                                                                                     <  1:303662/149‑1 (MQ=255)
  gcgcTTTTTAGGGCTGAATACAACCGCATCAGGATGGCGAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAg                                                                                                                                                   >  1:7441/1‑149 (MQ=255)
    gcTTTTTAGGGCTGAATACAACCGCATCAGGATGGCGAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGAt                                                                                                                                                 >  2:352471/1‑149 (MQ=255)
      tttttAGGGCTGAATACAACCGCATCAGGATGGCGAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGAt                                                                                                                                                 >  1:13264/1‑147 (MQ=255)
      tttttAGGGCTGAATACAACCGCATCAGGATGGCGAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGAt                                                                                                                                                 <  2:13264/147‑1 (MQ=255)
              gCTGAATACAACCGCATCAGGATGGCGAGTAAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGt                                                                                                                                       <  2:341162/149‑1 (MQ=255)
               cTGAATACAACCGCATCAGGATGGCGAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTg                                                                                                                                      >  3:442474/1‑149 (MQ=255)
                  aaTACAACCGCATCAGGATGGCGAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGaaa                                                                                                                                   >  1:236910/1‑149 (MQ=255)
                     aCAACCGCATCAGGATGGCGAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAAcc                                                                                                                                >  4:220928/1‑149 (MQ=255)
                      cAACCGCATCAGGATGGCGAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAAccc                                                                                                                               <  1:372756/149‑1 (MQ=255)
                           gCATCAGGATGGCGAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGAttt                                                                                                                          >  3:266104/1‑149 (MQ=255)
                               cAGGATGGCGAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTg                                                                                                                      <  3:345634/149‑1 (MQ=255)
                               cAGGATGGCGAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTg                                                                                                                      <  4:98851/149‑1 (MQ=255)
                                       cGAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGAttt                                                                                                              >  2:482183/1‑149 (MQ=255)
                                        gAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGTCAGaat                                                                                                                                                >  2:111001/1‑112 (MQ=255)
                                        gAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGTCAGaat                                                                                                                                                <  1:111001/114‑3 (MQ=255)
                                        gAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGcc                                                                                                                    >  3:654191/1‑142 (MQ=255)
                                        gAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCAGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGcc                                                                                                                    <  4:654191/142‑1 (MQ=255)
                                          gTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGa                                                                                                           <  1:352471/149‑1 (MQ=255)
                                            tAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGAtttttt                                                                                                                       <  4:325092/135‑1 (MQ=255)
                                            tAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGAtttttt                                                                                                                       >  3:325092/1‑135 (MQ=255)
                                              aaaaaTTCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCa                                                                                                       >  2:338761/1‑149 (MQ=255)
                                              aaaaaTCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCa                                                                                                       >  4:347346/1‑149 (MQ=255)
                                                   tCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCaaa                                                                                                  >  1:334058/1‑149 (MQ=255)
                                                           aTCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGAtttttt                                                                                                                       <  2:228679/120‑1 (MQ=255)
                                                           aTCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGAtttttt                                                                                                                       >  1:228679/1‑120 (MQ=255)
                                                           aTCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACg                                                                                          >  4:132091/1‑149 (MQ=255)
                                                            tCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACgg                                                                                         <  3:439851/149‑1 (MQ=255)
                                                                     gtgtCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGAttt                                                                                <  4:490779/149‑1 (MQ=255)
                                                                        tCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGa                                                                             <  4:442474/149‑1 (MQ=255)
                                                                        tCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGa                                                                             <  4:69921/149‑1 (MQ=255)
                                                                        tCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGa                                                                             <  2:236910/149‑1 (MQ=255)
                                                                         ccGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGAt                                                                            <  4:236343/149‑1 (MQ=255)
                                                                               gCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCa                                                                      <  2:302352/149‑1 (MQ=255)
                                                                                 gAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGa                                                                    >  2:571230/1‑149 (MQ=255)
                                                                                 gAAGCGTCCGGACCAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGa                                                                    >  3:155742/1‑149 (MQ=255)
                                                                                     cGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGAtttt                                                                <  1:482183/149‑1 (MQ=255)
                                                                                        ccGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGAc                                                             <  1:338761/149‑1 (MQ=255)
                                                                                              gagaAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGTCAGAATCTGGCggggat                                                                                                                                     <  3:18130/71‑4 (MQ=38)
                                                                                              gagaAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGTCAGAATCTGGCggggat                                                                                                                                     >  4:18130/1‑68 (MQ=38)
                                                                                              gagaAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGa                                                              <  1:217472/142‑1 (MQ=255)
                                                                                              gagaAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGa                                                              >  2:217472/1‑142 (MQ=255)
                                                                                                gaAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGAtttt                                                                <  3:430262/138‑1 (MQ=255)
                                                                                                gaAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGAtttt                                                                >  4:430262/1‑138 (MQ=255)
                                                                                                gaAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCaga                                                     >  3:539895/1‑149 (MQ=255)
                                                                                                               gACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCAGAGTCTGACAGTTgaga                                      >  3:462749/1‑149 (MQ=255)
                                                                                                                       aTGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCAGAGTc                                                  >  2:450634/1‑129 (MQ=255)
                                                                                                                       aTGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCAGAGTc                                                  <  1:450634/129‑1 (MQ=255)
                                                                                                                          ttttGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCa                                                                      >  3:86300/1‑106 (MQ=255)
                                                                                                                          ttttGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCa                                                                      <  4:86300/106‑1 (MQ=255)
                                                                                                                                ttttACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATTCGTCAGATTTTGACGGGTCAGAGTCTGACAGTTGAGAAAATGCCGCTGCCTGaa                     >  1:118976/1‑149 (MQ=255)
                                                                                                                                ttttACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCAGAGTCTGACAGTTGAGAAAATGCCGCTGCCTGaa                     <  2:334058/149‑1 (MQ=255)
                                                                                                                                     cTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCAGAGTCTGACAGTTGAGAAAATGCCGCTGCCTGAAGCTTc                <  3:347346/149‑1 (MQ=255)
                                                                                                                                        aCGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCAGAGTCTGACAGTTGAGAAAATGCCGCTGCCTGAAGCTTCGCa             >  2:427722/1‑149 (MQ=255)
                                                                                                                                          ggATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCAGAGTCTGACAGTTGAGAAAATGCCGCTGCCTGAAGCTTCGCAAc           >  3:276312/1‑149 (MQ=255)
                                                                                                                                            aTCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCAGAGTCTGACAGTTGAGAAAATGCCGCTGCCTGAAGCTTCGCAACGt         <  2:7441/149‑1 (MQ=255)
                                                                                                                                                 cGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAAtttt                                                                                              >  2:484677/1‑59 (MQ=255)
                                                                                                                                                 cGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAAtttt                                                                                              <  1:484677/59‑1 (MQ=255)
                                                                                                                                                   ccAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCAGAGTCTGACAGTTGAGAAAATGCCGCTGCCTGAAGCTTCGCAACGTTAAGCTg  >  4:6496/1‑149 (MQ=255)
                                                                                                                                                   ccAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCAGAGTCTGACAGTTGAGAAAATGCCGCTGCCTGAAGCTTCGCAACGTTAAGCTg  >  1:694315/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGCGCTTTTTAGGGCTGAATACAACCGCATCAGGATGGCGAGTTAAAAAATCCTGTTCATCCGTCTGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCAGAGTCTGACAGTTGAGAAAATGCCGCTGCCTGAAGCTTCGCAACGTTAAGCTG  >  NC_000913/2473266‑2473561

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