Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,012,563 (G)6→8 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.G100.0% 134.7 / NA 39pseudogene (56/279 nt)intGputative defective phage integrase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (18/21);  total (18/21)
*NC_0009132,012,5572.G94.9% 126.6 / 0.2 39pseudogene (56/279 nt)intGputative defective phage integrase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (17/20);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

CTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCA‑‑GGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAA  >  NC_000913/2012428‑2012686
                                                                                                                                  ||                                                                                                                                 
cTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCA‑‑ggggggggtt                                                                                                                         >  2:366885/1‑136 (MQ=255)
cTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCA‑‑ggggggggtt                                                                                                                         >  1:392093/1‑136 (MQ=255)
cTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCA‑‑g                                                                                                                                  <  2:159012/131‑1 (MQ=255)
cTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCA‑‑g                                                                                                                                  >  1:159012/1‑131 (MQ=255)
  cATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCA‑‑ggggggggttat                                                                                                                       <  1:580138/140‑7 (MQ=255)
  cATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCA‑‑ggggggggttat                                                                                                                       <  1:254321/140‑7 (MQ=255)
      tCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATcgcg                                                                                                                     <  2:121030/140‑1 (MQ=255)
          tCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGtata                                                                                                                 <  1:205812/140‑1 (MQ=255)
                            aaaaCATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGa                                                                                               >  1:528209/1‑140 (MQ=255)
                            aaaaCATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGaa                                                                                              >  1:469379/1‑140 (MQ=255)
                            aaaaCATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGaa                                                                                              <  2:469379/140‑1 (MQ=255)
                            aaaaCATCAACAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGaa                                                                                              >  1:788041/1‑140 (MQ=255)
                             aaaCATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGaa                                                                                              <  1:235190/140‑1 (MQ=255)
                              aaCATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGaaa                                                                                             <  1:93456/140‑1 (MQ=255)
                                cATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAAc                                                                                           >  1:760752/1‑140 (MQ=255)
                                    aacaaAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTc                                                                                       <  1:245753/140‑1 (MQ=255)
                                          gATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATc                                                                                 <  2:528209/140‑1 (MQ=255)
                                            ttAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAg                                                                               <  1:602867/140‑1 (MQ=255)
                                                 ttATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATc                                                                          <  2:392093/140‑1 (MQ=255)
                                                   aTCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGtt                                                                                        <  2:685469/124‑1 (MQ=255)
                                                   aTCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGtt                                                                                        >  1:685469/1‑124 (MQ=255)
                                                         agGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTaa                                                                  >  1:357779/1‑140 (MQ=255)
                                                         agGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTaa                                                                  >  2:716379/1‑140 (MQ=255)
                                                                   tGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAAtt                                                        >  1:394315/1‑140 (MQ=255)
                                                                      gTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAAc                                                                <  2:106349/129‑1 (MQ=255)
                                                                      gTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAAc                                                                >  1:106349/1‑129 (MQ=255)
                                                                      gTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATtgag                                                     >  1:536835/1‑140 (MQ=255)
                                                                            ggTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGtt                                               >  2:2001/1‑140 (MQ=255)
                                                                            ggTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGtt                                               >  2:246138/1‑140 (MQ=255)
                                                                                 tAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGccc                                          <  2:788041/140‑1 (MQ=255)
                                                                                      acaTCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCccatcc                                     <  2:394315/140‑1 (MQ=255)
                                                                                      acaTCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCccatcc                                     <  1:366885/140‑1 (MQ=255)
                                                                                            gACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGccc                                          >  2:102171/1‑129 (MQ=255)
                                                                                            gACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGccc                                          <  1:102171/129‑1 (MQ=255)
                                                                                              cGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCa                             <  1:8260/140‑1 (MQ=255)
                                                                                                  gtgATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCAcgcg                         <  1:246138/140‑1 (MQ=255)
                                                                                                         gtCAAATGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGGGGTGAGAt                  <  2:262814/140‑1 (MQ=255)
                                                                                                             aaTGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGa              >  1:662642/1‑140 (MQ=255)
                                                                                                               tGAGACGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATg            <  1:716379/140‑1 (MQ=255)
                                                                                                                  gaCGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTaa         >  1:414680/1‑140 (MQ=255)
                                                                                                                  gaCGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCC‑‑GCCCCACGCGGTGAGATCGGATGTAAtt       >  1:641580/1‑140 (MQ=255)
                                                                                                                   aCGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGAt                  <  1:737458/130‑1 (MQ=255)
                                                                                                                   aCGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGAt                  >  2:737458/1‑130 (MQ=255)
                                                                                                                    cGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCAcgcg                         >  1:695927/1‑122 (MQ=255)
                                                                                                                    cGCTCTGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCAcgcg                         <  2:695927/122‑1 (MQ=255)
                                                                                                                         tGATGTTCAGGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTaacaa  <  1:127759/140‑1 (MQ=255)
                                                                                                                                  ||                                                                                                                                 
CTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCA‑‑GGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAA  >  NC_000913/2012428‑2012686

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: