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.G93.9% 163.3 / 3.2 50pseudogene (56/279 nt)intGputative defective phage integrase
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base G (22/24);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

AAAAATCGCTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAA  >  NC_000913/2012420‑2012686
                                                                                                                                          |                                                                                                                                    
aaaaaTCGCTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAgg                                                                                                                                    <  1:87137/140‑1 (MQ=255)
  aaaTCGCTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAgggg                                                                                                                                 <  1:564970/140‑1 (MQ=255)
      cGCTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAgggggggt                                                                                                                           >  2:409248/1‑138 (MQ=255)
      cGCTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAgggggggt                                                                                                                           >  1:70901/1‑138 (MQ=255)
      cGCTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAgggg                                                                                                                                 >  2:366622/1‑136 (MQ=255)
      cGCTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAgggg                                                                                                                                 <  1:366622/136‑1 (MQ=255)
      cGCTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGATGAAGTGGCAAATGAGGCGCTGTGATGTTCAgggggggt                                                                                                                           >  2:198144/1‑138 (MQ=255)
         tCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAgggggggttat                                                                                                                        <  1:418125/140‑6 (MQ=255)
         tCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAgggggggttat                                                                                                                        <  2:436313/140‑6 (MQ=255)
                 aTCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGtata                                                                                                                 <  2:733465/140‑1 (MQ=255)
                  tCGCCCGGAATATGAATAAAAACATCAACAATGCTGAAATTATCCAGATGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATAc                                                                                                                <  2:330036/140‑1 (MQ=255)
                   cGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACg                                                                                                               >  1:646184/1‑140 (MQ=255)
                     cccGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGTCGCTCTGATGTGCAGGGGGGGTTATCGCGTATACGGa                                                                                                             >  1:330036/1‑140 (MQ=255)
                            tatGAATAAAAACATCAACAAACATTAAATTGTCCCGAGCTTACTGTTTCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGa                                                                                                      <  2:218037/140‑1 (MQ=255)
                                aaTAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGCTCAGGGGAGGTGATTGCGTATACGGATCGTGGATGGc                                                                                                  >  1:292245/1‑140 (MQ=255)
                                    aaaaCATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGaa                                                                                              <  2:70901/140‑1 (MQ=255)
                                    aaaaCATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGaa                                                                                              >  1:142565/1‑140 (MQ=255)
                                    aaaaCATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGaa                                                                                              >  1:623598/1‑140 (MQ=255)
                                         aTCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGt                                                                                         <  1:449845/140‑1 (MQ=255)
                                         aTCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGt                                                                                         <  1:522722/140‑1 (MQ=255)
                                         aTCAACAAAGATTAAATTATCCAGACGTTGCTGTTGCCTCTGTGGTCGTAACAACATCCGACGTGGTGATGCGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGt                                                                                         <  1:358123/140‑1 (MQ=255)
                                          tCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGtt                                                                                        <  1:409248/140‑1 (MQ=255)
                                          tCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGt                                                                                         >  2:608789/1‑139 (MQ=255)
                                          tCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGt                                                                                         <  1:608789/139‑1 (MQ=255)
                                                   aTTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGTGGTTATCGCGTAAACGGATCGTGGAGGGCAAAAAACGTTCAGGATCAg                                                                               >  2:226628/1‑140 (MQ=255)
                                                                gagGTTACTGTTGCGTACGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTaa                                                                  <  2:337501/140‑1 (MQ=255)
                                                                 agGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTaaa                                                                 >  1:658304/1‑140 (MQ=255)
                                                                 agGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTaaa                                                                 >  2:129005/1‑140 (MQ=255)
                                                                      aCTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGa                                                            >  1:406473/1‑140 (MQ=255)
                                                                             cGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATtgag                                                     >  1:599706/1‑140 (MQ=255)
                                                                               tCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATtgagtg                                                   >  1:403247/1‑140 (MQ=255)
                                                                                  gtgGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAATCGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTgagtgagt                                                <  1:129005/140‑1 (MQ=255)
                                                                                  gtgGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTgagtgagt                                                >  1:140965/1‑140 (MQ=255)
                                                                                    ggTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGttt                                              >  2:659144/1‑140 (MQ=255)
                                                                                        gTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGccc                                          <  2:658304/140‑1 (MQ=255)
                                                                                         tAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGcccc                                         >  2:585205/1‑140 (MQ=255)
                                                                                         tAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGATATCCATTTTATACGAGAAATCGAGTGAGTTTGcccc                                         >  1:304825/1‑140 (MQ=255)
                                                                                               caTCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCccatccat                                   >  1:79538/1‑140 (MQ=255)
                                                                                               caTCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCccatccat                                   >  1:58216/1‑140 (MQ=255)
                                                                                                 tCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCTTTTAAAAACAGAAATTGAGTGAGTTTGCCCCATCCATGc                                 >  2:307776/1‑140 (MQ=255)
                                                                                                 tCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGc                                 <  2:577014/140‑1 (MQ=255)
                                                                                                 tCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATAAGGTATCCGTTTGAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGc                                 >  2:217106/1‑140 (MQ=255)
                                                                                                     aCGTGGTGATGTGTCAAATGAGTCGCTCTTATGTTCATGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCa                             <  1:226628/140‑1 (MQ=255)
                                                                                                     aCGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCa                             >  1:387258/1‑140 (MQ=255)
                                                                                                           tgATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGt                       <  2:623598/140‑1 (MQ=255)
                                                                                                                tgtCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGAt                  <  2:140965/140‑1 (MQ=255)
                                                                                                                tgtCAAATGAGACGCTCTGATGTTCAGGGGGAGTTATCGAGTATACGGCTCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGAt                  <  2:304825/140‑1 (MQ=255)
                                                                                                                  tCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCAcgcg                         >  1:714880/1‑131 (MQ=255)
                                                                                                                  tCAAATGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCAcgcg                         <  2:714880/131‑1 (MQ=255)
                                                                                                                   cAAATGAGGCGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATTCATGCCCCACGCGGTGAGATCgg               <  2:403247/140‑1 (MQ=255)
                                                                                                                   cAAATGAGACGCTCTCATGTTCAGGAGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCgg               <  1:307776/140‑1 (MQ=255)
                                                                                                                     aaTGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGAt             >  2:180033/1‑140 (MQ=255)
                                                                                                                      aTGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATg            <  1:53127/140‑1 (MQ=255)
                                                                                                                      aTGAGACGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGAATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATg            <  2:292245/140‑1 (MQ=255)
                                                                                                                           aCGCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAAtt       >  2:9964/1‑140 (MQ=255)
                                                                                                                             gCTCTGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTaa     <  2:406473/140‑1 (MQ=255)
                                                                                                                               tctGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTaaca   <  2:58216/140‑1 (MQ=255)
                                                                                                                                ctGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTaacaa  <  1:553017/140‑1 (MQ=255)
                                                                                                                                ctGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTaacaa  <  1:659144/140‑1 (MQ=255)
                                                                                                                                ctGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTaacaa  <  2:387258/140‑1 (MQ=255)
                                                                                                                                ctGATGTTCAGGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTaacaa  <  1:180033/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                    
AAAAATCGCTCATATCCATCGCCCGGAATATGAATAAAAACATCAACAAAGATTAAATTATCCAGAGGTTACTGTTGCGTCTGTGGTCGTAACAACATCCGACGTGGTGATGTGTCAAATGAGACGCTCTGATGTTCAGGGGGGTTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAA  >  NC_000913/2012420‑2012686

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