New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 4625972 =NA (NA)2 (0.060)
+CGACCGCTGTACTCCACCAACCCCATTGCATTTGCCGC
4/174 NT 4.8% intergenic (+4519/+45) PP_5654/PP_4091 phosphotriesterase/transposase
?NC_002947 5401729 = 57 (1.160)intergenic (+1057/+111) PP_4743/PP_4745 hypothetical protein/transposase
Rejected: Frequency below/above cutoff threshold.

CCAGCGACAACCTGATTGCCGAGGTACTCACCCAAACACAGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGCGTGCGGCAGTCGCTCCAGGACGTGGCGCAGCCACGT  >  NC_002947/4625823‑4626118
                                                                                                                                                     |                                                                                                                                                  
ccAGCGACAACCTGATTGCCGAGGTACTCACCCAAACACAGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACg                                                                                                                                                    >  2:971494/1‑150 (MQ=255)
           cTGATTGCCGAGGTACTCACCCAAACACAGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCa                                                                                                                                         >  2:900155/1‑150 (MQ=255)
           cTGATTGCCGAGGTACTCACCCAAACACAGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCa                                                                                                                                         >  2:735783/1‑150 (MQ=255)
              aTTGCCGAGGTACTCACCCAAACACAGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGa                                                                                                                                      <  1:60699/150‑1 (MQ=255)
                 gCCGAGGTACTCACCCAGACGCAGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAACCCCTCAAAAGCAGAGCATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAAcc                                                                                                                                   <  1:1186069/150‑1 (MQ=14)
                  ccGAGGTACTCACCCAAACACAGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAAccc                                                                                                                                  >  1:367518/1‑150 (MQ=255)
                   cGAGGTACTCACCCAAACACAGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAAcccc                                                                                                                                 >  1:955306/1‑150 (MQ=255)
                   cGAGGTACTCACCCAAACACAGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAAcccc                                                                                                                                 >  1:118930/1‑150 (MQ=255)
                    gAGGTACTCACCCAAACACAGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCa                                                                                                                                <  1:900155/150‑1 (MQ=255)
                               ccAAACACAGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTaa                                                                                                                                                 >  2:455000/1‑122 (MQ=255)
                               ccAAACACAGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTaa                                                                                                                                                 <  1:455000/122‑1 (MQ=255)
                               ccAAACACAGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGCTCCATGAACCCCACCTACCACAAg                                                                                                                     >  2:204864/1‑150 (MQ=255)
                               ccAAACACAGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACGCCACCTGCCACAAg                                                                                                                     <  2:955306/150‑1 (MQ=255)
                                cAAACACAGCGCCGAGGGATGTACAACGCCGCAACGGGCGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCAACCCCACGCGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTCCCACAAgg                                                                                                                    <  2:942216/150‑1 (MQ=18)
                                     acaGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGaa                                                                                                               <  2:1262810/150‑1 (MQ=255)
                                      caGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAg                                                                                                              <  2:1082034/150‑1 (MQ=255)
                                        gcgcCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCACCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAg                                                                                                                     >  1:445557/1‑141 (MQ=255)
                                        gcgcCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCAACGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAg                                                                                                                     <  2:445557/141‑1 (MQ=255)
                                                       cAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGt                                                                                                                   <  1:332836/128‑1 (MQ=255)
                                                       cAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGt                                                                                                                   >  2:332836/1‑128 (MQ=255)
                                                            cctccAATGACGAAAACGCGAACCCCGAGAAGATCATCGGTCGCTGCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCCCAAGGTGGAAGCCAAAGTTACCGTGGCATctct                                                                                        <  1:547844/147‑1 (MQ=25)
                                                              gccAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTaa                                                                                                                                                 >  1:101480/1‑91 (MQ=39)
                                                              gccAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTaa                                                                                                                                                 <  2:101480/91‑1 (MQ=39)
                                                                      cGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGt                                                                              >  1:174880/1‑150 (MQ=255)
                                                                       gAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCt                                                                                                                             >  1:662258/1‑102 (MQ=255)
                                                                       gAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCt                                                                                                                             <  2:662258/102‑1 (MQ=255)
                                                                       gAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTg                                                                                                <  1:630538/131‑1 (MQ=255)
                                                                       gAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTg                                                                                                >  2:630538/1‑131 (MQ=255)
                                                                         aaaCGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTAc                                                                                                                           >  2:242309/1‑102 (MQ=255)
                                                                         aaaCGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTAc                                                                                                                           <  1:242309/102‑1 (MQ=255)
                                                                             gCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCCTTCGCTGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATgg                                                                       >  1:1043178/1‑150 (MQ=255)
                                                                                    tCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCagcagc                                                                <  2:78192/150‑1 (MQ=255)
                                                                                    tCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCagcagc                                                                <  2:174880/150‑1 (MQ=255)
                                                                                      gagaaGATCATCGGTGGCTTCCTTGAACTACAACACAAGGGCAAGACCTATTTGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAgcgc                                                              >  2:137508/4‑150 (MQ=18)
                                                                                              catcGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGtt                                                                             <  1:42753/127‑1 (MQ=255)
                                                                                              catcGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGtt                                                                             >  2:42753/1‑127 (MQ=255)
                                                                                                 cGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAg                                                                                                        >  2:517525/1‑97 (MQ=255)
                                                                                                 cGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAg                                                                                                        <  1:517525/97‑1 (MQ=255)
                                                                                                     cGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGCTCCATGAACCCCACCTCCCACCAGGTGGAAGCCAAAGTTACCGTCGCATCTCTGGCGAACAGTTCCATGGCCGCAGCGCTTCGTAGTCTGCTAc                                               <  2:367518/150‑1 (MQ=21)
                                                                                                     cGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCATCGTAGTCTTCTAc                                               <  2:118930/150‑1 (MQ=255)
                                                                                                        ttCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGa                                            <  1:971494/150‑1 (MQ=255)
                                                                                                           cTTGACCTACCCCACACGGGCAAGACCTATTCGGGACGCAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGa                                         >  1:4526/1‑150 (MQ=255)
                                                                                                           cTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGa                                         >  1:372488/1‑150 (MQ=255)
                                                                                                              gACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGc                                      >  2:399821/1‑150 (MQ=255)
                                                                                                                ccTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATg                                                                        >  2:1178589/1‑114 (MQ=255)
                                                                                                                ccTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATg                                                                        <  1:1178589/114‑1 (MQ=255)
                                                                                                                 cTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGCGTg                                   >  2:1052229/1‑150 (MQ=255)
                                                                                                                 cTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGCGTg                                   >  2:367676/1‑150 (MQ=255)
                                                                                                                   aCCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGCGTGCg                                 >  2:826163/1‑150 (MQ=255)
                                                                                                                     cccACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGaa                                                                                                               >  1:743335/1‑70 (MQ=255)
                                                                                                                     cccACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGaa                                                                                                               <  2:743335/70‑1 (MQ=255)
                                                                                                                      ccACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAg                                        >  1:507473/1‑140 (MQ=255)
                                                                                                                      ccACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAg                                        <  2:507473/140‑1 (MQ=255)
                                                                                                                      ccACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGCGTGCGGCa                              <  2:1043178/150‑1 (MQ=255)
                                                                                                                      ccACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGCGTGCGGCa                              <  1:826163/150‑1 (MQ=255)
                                                                                                                       cacaCGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGa                                                                                                                <  1:563428/67‑1 (MQ=255)
                                                                                                                       cacaCGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGa                                                                                                                >  2:563428/1‑67 (MQ=255)
                                                                                                                        acacGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGCGTGCGGCAGt                            >  2:879223/1‑150 (MQ=255)
                                                                                                                        acacGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGCGTGCGGCAGt                            <  1:216120/150‑1 (MQ=255)
                                                                                                                            gggCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTAACACAAGGTGGAAGCCAAAGTTACCGTGACATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGCGTGCGGCAGTCGCt                        <  1:1215268/150‑1 (MQ=255)
                                                                                                                                 aGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGCGTGCGGCAGTCGCTCCAgg                   >  2:895753/1‑150 (MQ=255)
                                                                                                                                 aGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGCGTGCGGCAGTCGCTCCAgg                   >  1:1201769/1‑150 (MQ=255)
                                                                                                                                       aTTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGCGTGCGGCAGTCGCTCCAGGACGTgg             <  1:367676/150‑1 (MQ=40)
                                                                                                                                            ggACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGCGTGCGGCAGTCGCTCCAGGACGTGGCGCAg        >  2:1012151/1‑150 (MQ=37)
                                                                                                                                                acaACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGCGTGCGGCAGTCGCTCCAGGACGTGGCGCAGCCAc    >  2:127464/1‑150 (MQ=33)
                                                                                                                                                 caACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCaaa                                                                                                         >  2:1084756/1‑48 (MQ=40)
                                                                                                                                                 caACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCaaa                                                                                                         <  1:1084756/48‑1 (MQ=40)
                                                                                                                                                  aaCGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGCGTGCGGCAGTCGCTCCAGGACGTGGCGCAGCCACGt  >  1:637152/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                                                                  
CCAGCGACAACCTGATTGCCGAGGTACTCACCCAAACACAGCGCCGAGGGATGTCCAACGCCGCCAATGACGAAAACGCTAACCTCGAAAAGATCATCGGTCGCTTCCTTGACCTACCCCACACGGGCAAGACCTATTCGGGACACAACGTAACCGATCCATGAACCCCACCTACCACAAGGTGGAAGCCAAAGTTACCGTGGCATCTCTGGCGAACAGTTCCATGGCAGCAGCGCTTCGTAGTCTGCTACCGATGAGGCGTGCGGCAGTCGCTCCAGGACGTGGCGCAGCCACGT  >  NC_002947/4625823‑4626118

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.