Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,012,563 (G)6→7 pseudogene (62/279 nt) intG → putative defective phage integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,012,5571.G97.6% 137.6 / ‑3.3 41pseudogene (56/279 nt)intGputative defective phage integrase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (19/21);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCAAAAATCGCTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTA  >  NC_000913/2012418‑2012682
                                                                                                                                            |                                                                                                                             
tCAAAAATCGCTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCa                                                                                                                                <  1:632719/140‑1 (MQ=255)
  aaaaaTCGCTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTCATGTTCAgg                                                                                                                             <  2:511262/140‑1 (MQ=255)
    aaaTCCTTCATATCCACCGCCGGGAATATGAAAAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAgggg                                                                                                                           <  1:793226/140‑1 (MQ=255)
         gCTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAgggggggtt                                                                                                                      <  2:818213/140‑4 (MQ=255)
         gCTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAgggggggtt                                                                                                                      <  1:9971/140‑4 (MQ=255)
         gCTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAgggggggtt                                                                                                                      <  1:519827/140‑4 (MQ=255)
          cTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAgggggggtta                                                                                                                     >  1:216885/1‑136 (MQ=255)
          cTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAgggggggtta                                                                                                                     >  2:121520/1‑136 (MQ=255)
           tCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAgggggggttat                                                                                                                    <  2:128537/140‑6 (MQ=255)
           tCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAgggggggttat                                                                                                                    <  1:568712/140‑6 (MQ=255)
                    tCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATcg                                                                                                                   >  2:96549/1‑133 (MQ=255)
                    tCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATcg                                                                                                                   <  1:96549/133‑1 (MQ=255)
                     cGCCCGGAATATGAATAAAAACATCAACAAAGAATAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGTTATCGCGTATACgg                                                                                                          <  1:807905/140‑1 (MQ=255)
                       cccGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGa                                                                                                         <  1:821387/140‑1 (MQ=255)
                       cccGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGa                                                                                                         >  2:662435/1‑140 (MQ=255)
                       cccGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGa                                                                                                         >  1:677657/1‑140 (MQ=255)
                            aaTATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTg                                                                                                    <  1:761295/140‑1 (MQ=255)
                                 gAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATgg                                                                                               >  1:453122/1‑140 (MQ=255)
                                            tCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGtt                                                                                    >  2:677961/1‑140 (MQ=255)
                                               acaaAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTAAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAg                                                                                 >  2:158778/1‑140 (MQ=255)
                                                caaAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAgg                                                                                <  2:804065/140‑1 (MQ=255)
                                                  aaGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTCCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGAt                                                                              >  2:648850/1‑140 (MQ=255)
                                                   aGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGt                                                                                     >  1:77564/1‑132 (MQ=255)
                                                   aGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGt                                                                                     <  2:77564/132‑1 (MQ=255)
                                                              tCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCAtt                                                                  >  1:259515/1‑140 (MQ=255)
                                                                cAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCAtttt                                                                <  1:6422/140‑1 (MQ=255)
                                                                  gagGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTaa                                                              <  1:662435/140‑1 (MQ=255)
                                                                         cTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGaa                                                       <  2:216885/140‑1 (MQ=255)
                                                                          tGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGaaa                                                      <  2:85828/140‑1 (MQ=255)
                                                                          tGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGaaa                                                      <  1:445802/140‑1 (MQ=255)
                                                                            ttGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAAtt                                                    >  2:332234/1‑140 (MQ=255)
                                                                               cGTCTGTGGTCGTAACAACATCCGACGTGGCGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGAATACGGATCGTGGATGGCAGAAAACGATCAGGATCAGGTATCCATTTTAAACCAGAAATtgag                                                 >  2:404126/1‑140 (MQ=255)
                                                                                 tCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATtgagtg                                               >  2:549133/1‑140 (MQ=255)
                                                                                          gTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGccc                                      <  1:121520/140‑1 (MQ=255)
                                                                                           tAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGccc                                      <  1:4516/139‑1 (MQ=255)
                                                                                           tAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGccc                                      >  2:4516/1‑139 (MQ=255)
                                                                                                       aCGTGGTGATGTGTCAAATGAGACGCTCTGATGTCCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCa                         <  1:648850/140‑1 (MQ=255)
                                                                                                        cGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCACGCCCCAc                        >  1:435166/1‑140 (MQ=255)
                                                                                                             tgATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGt                   <  1:332234/140‑1 (MQ=255)
                                                                                                                  tgtCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGAt              <  1:549133/140‑1 (MQ=255)
                                                                                                                     cAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCgg           <  2:529659/140‑1 (MQ=255)
                                                                                                                        aTGAGTCGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATg        <  1:404126/140‑1 (MQ=255)
                                                                                                                        aTGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATg        <  2:453122/140‑1 (MQ=255)
                                                                                                                        aTGAGACGCTCTGATGTTAAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATg        <  1:158778/140‑1 (MQ=255)
                                                                                                                            gaCGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAAt    >  2:416068/1‑140 (MQ=255)
                                                                                                                             aCGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAAtt   >  1:643673/1‑140 (MQ=255)
                                                                                                                              cGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTa  >  2:675331/1‑140 (MQ=255)
                                                                                                                              cGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTa  >  1:583093/1‑140 (MQ=255)
                                                                                                                                     aTGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAAt    <  2:601828/131‑1 (MQ=255)
                                                                                                                                     aTGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAAt    >  1:601828/1‑131 (MQ=255)
                                                                                                                                      tGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATg        <  1:504481/126‑1 (MQ=255)
                                                                                                                                      tGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATg        >  2:504481/1‑126 (MQ=255)
                                                                                                                                            |                                                                                                                             
TCAAAAATCGCTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTA  >  NC_000913/2012418‑2012682

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