Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 3,651,029 IS2 (+) +5 bp 91.4% intergenic (+367/‑258) arsC → / → yhiS arsenate reductase/putative uncharacterized protein YhiS

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3187426NA (NA)42 (1.350) 31/414 NT 91.1% noncoding (1331/1331 nt) IS2 repeat region
?NC_000913 3651029 = 4 (0.130)intergenic (+367/‑262) arsC/yhiS arsenate reductase/putative uncharacterized protein YhiS
* ? NC_000913 = 36510334 (0.130)43 (1.380) 35/414 NT 91.3% intergenic (+371/‑258) arsC/yhiS arsenate reductase/putative uncharacterized protein YhiS
?NC_000913 4498181 = NA (NA)noncoding (1/1331 nt) IS2 repeat region

TGACAGGCGTAGCGCATCAGGCGATTTTGCTATTTACACAGTACTCCCTGCGTGAGATTCCAATTATCGCGTCCAGCATGGTGTATCAGTGAGCTGCTTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATG  >  NC_000913/3650790‑3651293
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                    
tGACAGGCGTAGCGCATCAGGCGATTTTGCTATTTACACAGTACTCCCTGCGTGAGATTCCAATTATCGCGTCCAGCATGGTGTATCAGTGAGCTGCTTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctcaatacccaacaaa                                                                                                                                                                                                                                          >  1:138072‑M1/1‑244 (MQ=255)
   gagagtaacggaatagcagagagcttcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaacctttcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTc                                                                                                                                                                                                                                                                     <  1:95588‑M2/6‑1 (MQ=255)
               aTCAGGCGATTTTGCTATTTACACAGTACTCCCTGCGTGAGATTCCAATTATCGCGTCCAGCATGGTGTATCAGTGAGCTGCTTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaataacacttaaa                                                                                                                                                                                                                           >  1:28966‑M1/1‑229 (MQ=255)
                   ggCGATTTTGCTATTTACACAGTACTCCCTGCGTGAGATTCCAATTATCGCGTCCAGCATGGTGTATCAGTGAGCTGCTTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaacccccaaacaaacaaaaaaacctaaatcag                                                                                                                                                                                                                       >  2:284121‑M1/1‑225 (MQ=255)
                   cagagagcttcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCAtttt                                                                                                                                                                                                                       >  1:304131‑M2/221‑272 (MQ=255)
                           ttcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATc                                                                                                                                                                                                                                                           >  1:42182‑M2/213‑228 (MQ=255)
                           tGCTATTTACACAGTACTCCCTGCGTGAGATTCCAATTATCGCGTCCAGCATGGTGTATCAGTGAGCTGCTTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctcc                                                                                                                                                                                                                                                >  2:133281‑M1/1‑217 (MQ=255)
                            tcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCa                                                                                                                                                                                                                                                          <  2:42182‑M2/17‑1 (MQ=255)
                            gCTATTTACACAGTACTCCCTGCGTGAGATTCCAATTATCGCGTCCAGCATGGTGTATCAGTGAGCTGCTTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctcca                                                                                                                                                                                                                                               <  1:133281‑M1/239‑24 (MQ=255)
                                gaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttcagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAAc                                                                                                                                                                                                          >  1:400803‑M2/208‑272 (MQ=255)
                                gaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGAATTTTGATAGTTGTTAAc                                                                                                                                                                                                          >  1:383676‑M2/208‑272 (MQ=255)
                                gaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGAATTTTGATAGTTGTTAAc                                                                                                                                                                                                          >  1:250316‑M2/208‑272 (MQ=255)
                                gaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTAACACACTATCAGAGAACAAAATAACTTTTTAAAAAGAACGCATTTTTAAAGTTGTtaaa                                                                                                                                                                                                          >  2:149126‑M2/208‑271 (MQ=255)
                                        aaaagcgtgacgacatcgataacatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCt                                                                                                                                                                                                  <  2:304131‑M2/73‑1 (MQ=255)
                                               cTGCGTGAGATTCCAATTATCGCGTCCAGCATGGTGTATCAGTGAGCTGCTTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagact                                                                                                                                                                                           <  1:129201‑M1/272‑76 (MQ=255)
                                                   tacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACtt                                                                                                                                                                                                                                                                      >  2:217045‑M2/189‑193 (MQ=255)
                                                    acatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTc                                                                                                                                                                                                                                                                     <  1:217045‑M2/6‑1 (MQ=255)
                                                         agtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaattggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTa                                                                                                                                                                                 >  1:142257‑M2/183‑272 (MQ=255)
                                                                tgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGAATTCCTTTACGAAACTATCAATAACCcttt                                                                                                                                                                          >  2:375232‑M2/176‑272 (MQ=255)
                                                                  acagaaccagacgggttaacggcagcaaagaaccttgcagaggcggtcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCttttt                                                                                                                                                                         <  1:182735‑M2/98‑1 (MQ=255)
                                                                   ccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAAAAACTTTTTAAAATGAGCGCATTTTGATAGTTTTTAAAGATTTCCTTTATGATAACATCACTaccccttctct                                                                                                                                                                       >  2:40848‑M2/173‑268 (MQ=255)
                                                                    caaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACg                                                                                                                                                                                             <  1:40848‑M2/78‑1 (MQ=255)
                                                                     cgTCCAGCATGGTGTATCAGTGAGCTGCTTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAAAAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacgtaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgat                                                                                                                                                                      <  1:284121‑M1/271‑97 (MQ=255)
                                                                      aaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTaca                                                                                                                                                                    <  1:149126‑M2/103‑1 (MQ=255)
                                                                                    acggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATAACATCCttcacctttttttccaccacaacaacttct                                                                                                                                                      >  2:182735‑M2/156‑259 (MQ=255)
                                                                                    aTCAGTGAGCTGCTTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTCACACACTATc                                                                                                                                                                                                                                                           >  1:128230/1‑171 (MQ=255)
                                                                                     tCAGTGAGCTGCTTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTCACACACTATCa                                                                                                                                                                                                                                                          <  2:128230/171‑1 (MQ=255)
                                                                                         tGAGCTGCTTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaaaagatgaatgaattcttagggccggggaa                                                                                                                                                 >  2:117650‑M1/1‑155 (MQ=255)
                                                                                             tgggttagttTTCAGCGTTACAGCGAATTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagccattggagaacagatgattgatgtcttagggccggagaaacgc                                                                                                                                             <  2:28966‑M1/268‑122 (MQ=255)
                                                                                             aagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAaaattacaggtgggggaa                                                                                                                                             >  2:369496‑M2/147‑264 (MQ=255)
                                                                                             aagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGAt                                                                                                                                             >  1:246077‑M2/147‑272 (MQ=255)
                                                                                              tgctTAGTTATCAGCGTTACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgc                                                                                                                                             <  2:48883‑M1/271‑122 (MQ=255)
                                                                                              tgctTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagttattggaaaaaagaataatgaagtctaaaggcctggaaaacaca                                                                                                                                            >  2:175561‑M1/1‑150 (MQ=255)
                                                                                              tgctTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgc                                                                                                                                             >  1:48883‑M1/1‑150 (MQ=255)
                                                                                              tgctTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgat                                                                                                                                                                  >  2:298978‑M1/1‑150 (MQ=255)
                                                                                              tgctTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaa                                                                                                                                                                                                           >  1:382247‑M1/1‑150 (MQ=255)
                                                                                               gctTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatg                                                                                                                                                                 <  1:298978‑M1/250‑102 (MQ=255)
                                                                                               gctTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaat                                                                                                                                                                                                          <  2:382247‑M1/209‑61 (MQ=255)
                                                                                                ctTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcaga                                                                                                                                          >  2:330351‑M1/1‑148 (MQ=255)
                                                                                                  tAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttt                                                                                                                                                                                                 >  2:353955‑M1/1‑146 (MQ=255)
                                                                                                   aGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagttttt                                                                                                                                                                                                <  1:353955‑M1/215‑71 (MQ=255)
                                                                                                      tATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggaaaaacgcaaaagggct                                                                                                                                    >  2:153838‑M1/1‑142 (MQ=255)
                                                                                                          atcgaTACAGCGAAGTATTTTAAAGACGTAATGGATTATGATTTATCAGTGGTTGACATAACAAATGATTAAATAATTATAAGTTACATTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactaggcattggagaacagaggattgatgtcttagggccggagaaacgcagactgcgtacc                                                                                                                                 <  1:175561‑M1/269‑134 (MQ=255)
                                                                                                             gATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagt                                                                                                                                                                                                                      >  1:419553‑M1/1‑135 (MQ=255)
                                                                                                             gATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagt                                                                                                                                                                                                                      <  2:419553‑M1/183‑49 (MQ=255)
                                                                                                                         tagtTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagac                                                                                                                                         >  2:327603‑M1/1‑123 (MQ=255)
                                                                                                                         taaccctatggcaaccgcatagtgcgctggtttatcgctggccacgggaatatcggcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCa                                                                                                                 <  2:246077‑M2/154‑1 (MQ=255)
                                                                                                                          agtTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacg                                                                                                                                        <  1:327603‑M1/248‑127 (MQ=255)
                                                                                                                           gtTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcg                                                                                                                 <  1:153838‑M1/270‑150 (MQ=255)
                                                                                                                            tttAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaa                                                                                                                                                                                                           >  1:44452‑M1/1‑120 (MQ=255)
                                                                                                                             ttAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaat                                                                                                                                                                                                          <  2:44452‑M1/179‑61 (MQ=255)
                                                                                                                                     atccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCtttt                                                                                                                                                                          >  2:428290‑M2/107‑203 (MQ=255)
                                                                                                                                      tccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCttttt                                                                                                                                                                         <  1:428290‑M2/98‑1 (MQ=255)
                                                                                                                                      gcggctttgggcgcggggttatcgctctccacgggaatattttcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGTGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACg                                                                                                    <  2:250316‑M2/167‑1 (MQ=255)
                                                                                                                                                   tttATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattgaagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaaga                                                                                                                    >  1:248799‑M1/1‑97 (MQ=255)
                                                                                                                                                    ttATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattgaagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagat                                                                                                                   <  2:248799‑M1/243‑148 (MQ=255)
                                                                                                                                                    tgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGAtt                                                                                                                                                                                                      >  1:137733‑M2/92‑160 (MQ=255)
                                                                                                                                                     gggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGAttt                                                                                                                                                                                                     <  2:137733‑M2/70‑1 (MQ=255)
                                                                                                                                                           tGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgac                                                                               <  1:117650‑M1/272‑184 (MQ=255)
                                                                                                                                                              ctcgccacgggaatatctgcggcagcgtgcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTaga                                                                              <  1:369496‑M2/189‑1 (MQ=255)
                                                                                                                                                                   cacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTacaca                                                                                                                                                                  >  1:425424‑M2/77‑181 (MQ=255)
                                                                                                                                                                   cacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTacaca                                                                                                                                                                  <  2:425424‑M2/105‑1 (MQ=255)
                                                                                                                                                                   cacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTAc                                                                                                                                                        >  1:363648‑M2/77‑191 (MQ=255)
                                                                                                                                                                   cacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTatat                                                                                                                   >  1:332025‑M2/77‑228 (MQ=255)
                                                                                                                                                                    acgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATc                                                                                                                  <  2:332025‑M2/153‑1 (MQ=255)
                                                                                                                                                                      acaaATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcg                                                                    <  1:330351‑M1/272‑195 (MQ=255)
                                                                                                                                                                          attattAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagac                                                                                                                                                                                            >  2:427768‑M1/1‑74 (MQ=255)
                                                                                                                                                                          attattAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATc                                                                                                                  >  1:212565/1‑222 (MQ=255)
                                                                                                                                                                           ttattaAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagact                                                                                                                                                                                           <  1:427768‑M1/148‑76 (MQ=255)
                                                                                                                                                                           ttattaAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCa                                                                                                                 <  2:212565/222‑1 (MQ=255)
                                                                                                                                                                                   gcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAAt                                                       >  1:214471‑M2/61‑272 (MQ=255)
                                                                                                                                                                                    cagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAAtt                                                      <  1:375232‑M2/213‑1 (MQ=255)
                                                                                                                                                                                                       ccATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccggc                                                           >  2:117722‑M1/1‑45 (MQ=255)
                                                                                                                                                                                                        cATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccggca                                                          <  1:117722‑M1/248‑205 (MQ=255)
                                                                                                                                                                                                        aagtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTag                                                                               >  2:390138‑M2/40‑227 (MQ=255)
                                                                                                                                                                                                         agtgataacagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTaga                                                                              <  1:390138‑M2/189‑1 (MQ=255)
                                                                                                                                                                                                             tattattaCACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataag                                                                                                                                                                                                                       >  2:121175‑M1/1‑39 (MQ=255)
                                                                                                                                                                                                              attattaCACATCTAAAACATCACAAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagt                                                                                                                                                                                                                      <  1:121175‑M1/86‑49 (MQ=255)
                                                                                                                                                                                                                 cagatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTTTTCAAATTAACCAATTTAAAAATCTTAGCATAAACCa                         >  2:342153‑M2/31‑272 (MQ=255)
                                                                                                                                                                                                                  agatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATc                                                                                                                  >  1:360812‑M2/30‑182 (MQ=255)
                                                                                                                                                                                                                  agatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGa                                                                                                           >  1:222376‑M2/30‑189 (MQ=255)
                                                                                                                                                                                                                   gatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCa                                                                                                                 <  2:360812‑M2/154‑1 (MQ=255)
                                                                                                                                                                                                                   gatgtctggaaatataggggcaaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGat                                                                                                          <  2:222376‑M2/161‑1 (MQ=255)
                                                                                                                                                                                                                           gaaatataggggcaaatccaAACTTTAAACACTATCAGAGAACATAATAACTTTTTAAAATGAGCTCATTTTGATAGTTGTTAACGATTTCCTTTACGATAATATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAAt               <  1:342153‑M2/252‑1 (MQ=255)
                                                                                                                                                                                                                                     caAAAAAAATCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatacta                                                                                                                                                                                                      >  2:85245‑M1/1‑15 (MQ=255)
                                                                                                                                                                                                                                      aaaaaaaaTCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactag                                                                                                                                                                                                     <  1:85245‑M1/79‑66 (MQ=255)
                                                                                                                                                                                                                                       aaaaaaaTCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccggca                                                          <  2:375929‑M1/217‑205 (MQ=255)
                                                                                                                                                                                                                                       aaaaaaaTCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgccc                                                              >  1:401718‑M1/1‑13 (MQ=255)
                                                                                                                                                                                                                                        aaatccaCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGATCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCAAAAAAAAAAAGAATAAAAAGGAATTTTatt  >  2:311546‑M2/8‑271 (MQ=255)
                                                                                                                                                                                                                                        aaaaaaTCACTTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccg                                                             <  2:401718‑M1/213‑202 (MQ=255)
                                                                                                                                                                                                                                                  tttagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggat                                                                                  >  1:319631‑M1/1‑2 (MQ=255)
                                                                                                                                                                                                                                                  tttagactggccccctgaatctccagacaaccaatatcacttaaataagggatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagcttt                                                                                             <  2:249420‑M1/171‑170 (MQ=255)
                                                                                                                                                                                                                                                   ttagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatg                                                                                 <  2:319631‑M1/182‑182 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                    
TGACAGGCGTAGCGCATCAGGCGATTTTGCTATTTACACAGTACTCCCTGCGTGAGATTCCAATTATCGCGTCCAGCATGGTGTATCAGTGAGCTGCTTAGTTATCAGCGATACAGCGCAGTAGTTTAAAGACGTACTGGATTATGATTTATCAGTGGTTTACACAACAAATTATTAAATAATTATAAGATACCTTAACCCATTTTATTATTACACATCTAAAACATCACAAAAAAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATG  >  NC_000913/3650790‑3651293

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ 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.