New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4083704 =79 (1.020)3 (0.040) 3/252 NT 4.0% coding (2119/3051 nt) fdoG formate dehydrogenase‑O, large subunit
?NC_000913 4083730 = 70 (0.950)coding (2093/3051 nt) fdoG formate dehydrogenase‑O, large subunit
Rejected: Frequency below/above cutoff threshold.

TTATCACGGCGTGCCATCTGGTTGCCTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCA  >  NC_000913/4083556‑4083844
                                                                                                                                                    |                                                                                                                                            
ttATCACGGCGTGCCATCTGGTTGCCTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATc                                                                                                                                              >  3:299109/1‑149 (MQ=255)
ttATCACGGCGTGCCATCTGGTTGCCTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATc                                                                                                                                              <  3:243856/149‑1 (MQ=255)
acggcgagatcctctccggtatcttcctgcgcttgcgcaagatgtatgccgaacagggtggcgcgaacccggaccaggtgctgaacatgacatggaactacgccattccgcatgagccatcttcagaagaagtggcgatggagagcaac                                                                                                                                              <  3:449345‑M1/1‑1 (MQ=255)
  ggcgagatcctctccggtatcttcctgcgcttgcgcaagatgtatgccgaacagggtggcgcgaacccggaccaggtgctgaacatgacatggaactacgccattccgcatgagccatcttcagaagaagtggcgatggagagcaaCgg                                                                                                                                            <  3:581718‑M1/3‑1 (MQ=255)
   tCACGGCGTGCCATCTGGTTGCCTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGt                                                                                                                                           >  2:113352/1‑149 (MQ=255)
   tCACGGCGTGCCATCTGGTTGCCTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGt                                                                                                                                           >  1:93728/1‑149 (MQ=255)
    cACGGCGTGCCATCTGGTTGCCTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTa                                                                                                                                          <  2:73869/149‑1 (MQ=255)
    cACGGCGTGCCATCTGGTTGCCTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTa                                                                                                                                          >  2:329457/1‑149 (MQ=255)
     aCGGCGTGCCATCTGGTTGCCTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTaa                                                                                                                                         <  1:119396/149‑1 (MQ=255)
     aCGGCGTGCCATCTGGTTGCCTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTaa                                                                                                                                         >  3:414721/1‑149 (MQ=255)
     aCGGCGTGCCATCTGGTTGCCTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTaa                                                                                                                                         <  2:594778/149‑1 (MQ=255)
              cATCTGGTTGCCTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCa                                                                                                                                >  4:130268/1‑149 (MQ=255)
               aTCTGGTTGCCTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAg                                                                                                                               >  2:298066/1‑149 (MQ=255)
                tCTGGTTGCCTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAgc                                                                                                                              >  3:664200/1‑149 (MQ=255)
                    gTTGCGTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCTGCCAGCGCCt                                                                                                                          >  4:532824/1‑149 (MQ=255)
                     ttGCCTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTaa                                                                                                                                         <  1:53914/133‑1 (MQ=255)
                     ttGCCTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTaa                                                                                                                                         >  2:53914/1‑133 (MQ=255)
                         cTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCg                                                                                                                     >  1:562433/1‑149 (MQ=255)
                          ttCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGt                                                                                                                    <  3:112039/149‑1 (MQ=255)
                           tCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCTACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGtt                                                                                                                   <  4:525186/149‑1 (MQ=255)
                              ggCGTCCAGCTACCCGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGct                                                                                                                >  3:49354/1‑149 (MQ=255)
                               cttgcgcaagatgtatgccgaacagggtggcgcgaacccggaccaggtgctgaacatgacatggaactacgccattccgcatgagccatcttcagaagaagtggcgatggagagcaaCGGTAATATCGGCCAGCGCCTTACCGTTGctc                                                                                                               >  2:636124‑M1/118‑149 (MQ=255)
                                  tCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGc                                                                                                                                  <  4:344317/127‑1 (MQ=255)
                                  tCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGc                                                                                                                                  >  3:344317/1‑127 (MQ=255)
                                    cAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCAt                                                                                                          <  3:189530/149‑1 (MQ=255)
                                    cAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCAt                                                                                                          >  2:468586/1‑149 (MQ=255)
                                     aGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATc                                                                                                         <  1:641325/149‑1 (MQ=255)
                                        tACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGcc                                                                                                      <  1:263350/149‑1 (MQ=255)
                                         aCCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCa                                                                                                     >  3:26696/1‑149 (MQ=255)
                                           cGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCAct                                                                                                   >  4:529688/1‑149 (MQ=255)
                                           cGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCAct                                                                                                   >  1:89959/1‑149 (MQ=255)
                                             gCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCActtc                                                                                                 >  1:451340/1‑149 (MQ=255)
                                             gCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCActtc                                                                                                 >  2:26778/1‑149 (MQ=255)
                                             gCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCActtc                                                                                                 >  4:247746/1‑149 (MQ=255)
                                               gAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCActtctt                                                                                               <  1:222305/149‑1 (MQ=255)
                                               gAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCActtctt                                                                                               >  4:171839/1‑149 (MQ=255)
                                                 aaaTCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTtcttct                                                                                             <  1:230494/149‑1 (MQ=255)
                                                  aaTCCAGGCGCCACCGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTg                                                                                            <  1:176296/149‑1 (MQ=255)
                                                  aaTCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTg                                                                                            >  1:165678/1‑149 (MQ=255)
                                                      cagcagCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGa                                                                                        <  1:270187/149‑1 (MQ=255)
                                                       agcagcCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGAt                                                                                       >  2:152710/1‑149 (MQ=255)
                                                        gcagcCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATg                                                                                      >  3:145355/1‑149 (MQ=255)
                                                          agcCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTa                                                                                                                        >  3:535469/1‑113 (MQ=255)
                                                          agcCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTa                                                                                                                        <  4:535469/113‑1 (MQ=255)
                                                            ccACAGGAGGTTGTACCGTCATCGCGAAATTGGGCGAACGAACTAAGTTGTTGGCCATTCTTGACGATAACCGCGCCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGATGATGGCTc                                                                                  >  3:423181/1‑149 (MQ=255)
                                                              gcgaacccggaccaggtgctgaacatgacatggaactacgccattccgcatgagccatcttcagaagaagtggcgatggagagcaaCGGTAATATCGGCCAGCGCCTTACCGTTGc                                                                                                                 <  4:675990‑M1/30‑1 (MQ=255)
                                                              gcgaacccggaccaggtgctgaacatgacatggaactacgccattccgcatgagccatcttcagaagaagtggcgatggagagcaaCGGTAATATCGGCCAGCGCCTTACCGTTGc                                                                                                                 >  3:675990‑M1/87‑116 (MQ=255)
                                                              acaGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCAt                                                                                <  2:93728/149‑1 (MQ=255)
                                                              acaGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCAt                                                                                >  2:138056/1‑149 (MQ=255)
                                                              acaGGAGGTTGTACCGTCAACGCGCAATTGGGCGAACGAAATAAGTTGTTGGCCTTTCTGGGCGATAACCGCACCGGTTGCAGGACCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTGCATCACCACTTCTTCTGAAGATGGCTCAt                                                                                <  4:172822/149‑1 (MQ=255)
                                                                  gaggTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGt                                                                                                                                           <  3:195400/86‑1 (MQ=255)
                                                                  gaggTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGt                                                                                                                                           >  4:195400/1‑86 (MQ=255)
                                                                    ggTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGaa                                                                          <  4:548982/149‑1 (MQ=255)
                                                                           ccGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGTTCATGTGGAATGGCGTa                                                                   <  1:364287/149‑1 (MQ=255)
                                                                            cGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGAt                                                                                       <  4:483930/128‑1 (MQ=255)
                                                                            cGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGAt                                                                                       <  3:509426/128‑1 (MQ=255)
                                                                            cGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGAt                                                                                       >  3:483930/1‑128 (MQ=255)
                                                                            cGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGAt                                                                                       >  4:509426/1‑128 (MQ=255)
                                                                            cGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAg                                                                  >  4:615716/1‑149 (MQ=255)
                                                                                  cgcgCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCcat                                                            <  2:498862/149‑1 (MQ=255)
                                                                                     gcAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCcatgtc                                                         <  2:451340/149‑1 (MQ=255)
                                                                                       aaTTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCcatgtcat                                                       <  3:634070/149‑1 (MQ=255)
                                                                                             gCGAAGGAACTAAGTTGTTGGCCTTTCTTGACGATAACGGCAGCGGTTTCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAg                                                 >  2:234340/1‑149 (MQ=255)
                                                                                              cgaacgaaCTAAGTTGTTGGCCTTTCTTGACGATAACCTCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGc                                                <  4:616456/149‑1 (MQ=255)
                                                                                                aacgaacTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCCCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCAc                                              <  2:217246/149‑1 (MQ=255)
                                                                                                  cgaacTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCt                                            <  3:171839/149‑1 (MQ=255)
                                                                                                   gaacTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTg                                           >  4:52948/1‑149 (MQ=255)
                                                                                                     acTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGt                                         >  2:279763/1‑149 (MQ=255)
                                                                                                       tAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTcc                                       <  3:57878/149‑1 (MQ=255)
                                                                                                       tAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTcc                                       >  3:545720/1‑149 (MQ=255)
                                                                                                       tAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTcc                                       >  2:647653/1‑149 (MQ=255)
                                                                                                         aGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCgg                                     <  2:161200/149‑1 (MQ=255)
                                                                                                         aGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCgg                                     <  2:511808/149‑1 (MQ=255)
                                                                                                         aGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCgg                                     <  1:367030/149‑1 (MQ=255)
                                                                                                               tgGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTcgc                               <  2:56635/149‑1 (MQ=255)
                                                                                                               tgGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTcgc                               >  3:103457/1‑149 (MQ=255)
                                                                                                                  ccTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGcc                            <  1:152710/149‑1 (MQ=255)
                                                                                                                   cTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCa                           >  4:74062/1‑149 (MQ=255)
                                                                                                                   cTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCa                           >  3:229/1‑149 (MQ=255)
                                                                                                                   cTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCa                           >  2:150165/1‑149 (MQ=255)
                                                                                                                    tttCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCAc                          <  4:452496/149‑1 (MQ=255)
                                                                                                                    tttCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCAc                          <  4:664200/149‑1 (MQ=255)
                                                                                                                        ttGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTg                      >  2:24592/1‑149 (MQ=255)
                                                                                                                             gATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCgg                 >  3:194877/1‑149 (MQ=255)
                                                                                                                               tAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCa               <  1:150165/149‑1 (MQ=255)
                                                                                                                                 aCCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATa             <  1:26778/149‑1 (MQ=255)
                                                                                                                                  ccGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATAc            >  4:449345/1‑149 (MQ=255)
                                                                                                                                        cGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTg      <  2:197801/149‑1 (MQ=255)
                                                                                                                                         ggTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTgc     <  3:74062/149‑1 (MQ=255)
                                                                                                                                         ggTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTgc     <  2:89959/149‑1 (MQ=255)
                                                                                                                                         ggTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTgc     <  1:45300/149‑1 (MQ=255)
                                                                                                                                            tGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCGGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCa  >  1:653916/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
TTATCACGGCGTGCCATCTGGTTGCCTTCCGGCGTCCAGCTACCGGCGAAAATCCAGCAGCCACAGGAGGTTGTACCGTCATCGCGCAATTGGGCGAACGAACTAAGTTGTTGGCCTTTCTTGACGATAACCGCACCGGTTGCCGGATCGGTAATATCGGCCAGCGCCTTACCGTTGCTCTCCATCGCCACTTCTTCTGAAGATGGCTCATGCGGAATGGCGTAGTTCCATGTCATGTTCAGCACCTGGTCCGGGTTCGCGCCACCCTGTTCGGCATACATCTTGCGCA  >  NC_000913/4083556‑4083844

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