Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 1,426,494 IS2 (+) +5 bp 15.3% pseudogene (41‑45/933 nt) tmpR → Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1426494 =105 (0.800)22 (0.160) 18/176 NT 16.5% pseudogene (41/933 nt) tmpR Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein
?NC_000913 = 1651548 NA (NA)noncoding (706/706 nt) IS2 repeat region
* ? NC_000913 = 1426498105 (0.790)16 (0.110) 15/178 NT 12.6% pseudogene (45/933 nt) tmpR Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein
?NC_000913 = 2070271 NA (NA)noncoding (1/1331 nt) IS2 repeat region
Rejected: Frequency below cutoff threshold.

AGTCACTTTCCCTGCGCTATTTGCAGCAGTTTCCCCAGCCTGCGTCATTTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCt‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1426591‑1426494
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTAT  <  NC_000913/1651548‑1651449
                                                                                                                                                                                                      
AGTCACTTTCCCTGCGCTATTTGCAGCAGTTTCCCCAGCCTGCGTCATTTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTG                                                                                                    <  2:2967176/100‑1
   CACTTTCCCTGCGCTATTTGCAGCAGTTTCCCCAGCCTGCGTCATTTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGAT                                                                                                 >  1:3495769/1‑100
           CTGCGCTATTTGCAGCAGTTTCCCCAGCCTGCGTCATTTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCT                                                                                         >  2:85431/1‑100
           CTGCGCTATTTGCAGCAGTTTCCCCAGCCTGCGTCATTTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCT                                                                                         <  2:3866295/100‑1
             GCGCTATTTGCAGCAGTTTCCCCAGCCTGCGTCATTTTGACTAATCAGGAGGTCAGACTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCC                                                                                          >  1:1908034/1‑97
                 TATTTGCAGCAGTTTCCCCAGCCTGCGTCATTTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTT                                                                                   >  2:2299254/1‑100
                  ATTTGCAGCAGTTTCCCCAGCCTGCGTCATTTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTC                                                                                  <  2:2399949/100‑1
                  ATTTGCAGCAGTTTCCCCAGCCTGCGTCATCTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTC                                                                                  <  2:2391977/100‑1
                       CAGCAGTTTCCCCAGCCTGCGTCATTTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTCCAGAC                                                                             <  1:1320893/100‑1
                             TTTCCCCAGCCTGCGTCATTTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTCCAGACATCTGT                                                                       <  2:3723438/100‑1
                                   CAGCCTGCGTCATTTTGACTAATGAGGAGGTCAGATTGTTAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCAC                                                                 <  2:1143831/100‑1
                                   CAGCCTGCGTCATTTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCAC                                                                 <  1:2892241/100‑1
                                     GCCTGCGTCATTTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTT                                                               <  2:1995466/100‑1
                                     GCCTGCGTCATTTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTT                                                               <  2:1534575/100‑1
                                     GCCTGCGTCATTTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTT                                                               <  2:61613/100‑1
                                        TGCGTCATTTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAAC                                                            <  2:1753122/100‑1
                                              ATTTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTA                                                      >  2:2803713/1‑100
                                                    ACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCGGTGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCC                                                <  1:2605335/100‑1
                                                       AATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGC                                             >  2:1377716/1‑100
                                                               GGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACCTAACCCATTACAAGCCCGCTGCCGCAG                                     <  1:1666468/100‑1
                                                                   AGATTGTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATAT                                 <  2:2509025/100‑1
                                                                        GTCAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCG                            <  2:1726798/100‑1
                                                                          CAGCATTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGT                           >  1:1610865/1‑99
                                                                               TTATTTTTCGCTCCAGTGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAG                     >  2:465991/1‑100
                                                                                             AGTGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGC       <  1:3531414/100‑1
                                                                                               TGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCAC     <  2:3165155/100‑1
                                                                                                GCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCACAGATATTCCCGTGGCGAGCGATAACCCAGCGCACT    >  2:3292215/1‑100
                                                                                                 CTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTA   >  2:3240013/1‑100
                                                                                                 CTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTA   >  1:1451413/1‑100
                                                                                                 CTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACT    >  2:2069542/1‑99
                                                                                                 CTGGATTTGCCCCTATATTTCCAGACATCTGTTAGCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTA   <  1:2042504/100‑1
                                                                                                  TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTAT  >  1:1091355/1‑100
                                                                                                  TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTAT  <  2:73355/100‑1
                                                                                                  TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTAT  >  2:2922072/1‑100
                                                                                                  TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTAT  >  2:2896896/1‑100
                                                                                                  TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTAT  <  1:2551708/100‑1
                                                                                                  TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTAT  <  2:1992622/100‑1
                                                                                                  TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTAT  >  2:123321/1‑100
                                                                                                  TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTAT  <  2:1225198/100‑1
                                                                                                  TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTAT  <  1:3245657/100‑1
                                                                                                                                                                                                      
AGTCACTTTCCCTGCGCTATTTGCAGCAGTTTCCCCAGCCTGCGTCATTTTGACTAATGAGGAGGTCAGATTGTCAGCATTATTTTTCGCTCCAGTGCt‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1426591‑1426494
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTAT  <  NC_000913/1651548‑1651449

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.