New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1789000 =50 (0.840)3 (0.050) 3/266 NT 5.6% coding (566/1047 nt) aroH 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible
?NC_000913 1789053 = 53 (0.930)coding (619/1047 nt) aroH 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible

CGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACC  >  NC_000913/1788853‑1789147
                                                                                                                                                   |                                                                                                                                                   
aaactcgtgcagcgtatcgcaggctgcggcgatatcatcggcatgataattccgttatttgccgccacgcaaaataatgtggccatacgggttgccgctggtctgatagatggtcctctgaccatttttgtatggcgagaggaccatat                                                                                                                                                    <  1:245592‑M1/2‑1 (MQ=255)
  actcgtgcagcgtatcgcaggcttcggcgatatcatcggcatgataattcggttttttgccgcaacgcataataatgtggccatacgggttgccgctggtctgatagatcgtcatctgaccatttttgtctggcgagaggaacatATgg                                                                                                                                                  >  2:498374‑M1/146‑149 (MQ=255)
    cGAGGTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCa                                                                                                                                                >  2:923909/1‑149 (MQ=255)
         tCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATAcg                                                                                                                                           <  1:884425/149‑1 (MQ=255)
              gATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGAt                                                                                                                                      >  1:731947/1‑149 (MQ=255)
                    gTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCtgt                                                                                                                                <  2:306048/149‑1 (MQ=255)
                    gTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATACACCGCGAAATGGCTTCGGCACTCTCCTTTCCGGTAGGTTTTAAAAATGGTACAGATGGCAATACGCGGATTGCagt                                                                                                                                >  2:887987/1‑149 (MQ=255)
                                tttATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCcg                                                                                                                    >  1:389265/1‑149 (MQ=255)
                                   catcggcatgataattcggttttttgccgccacgcataataatgtggccatacgggttgccgctggtctgatagatggtcatctgaccatttttgtctgtcgagaggaacatATGGCTGGCGCGGATTGCTGTGGATGCTATCCGCGCa                                                                                                                >  1:27673‑M1/113‑149 (MQ=255)
                                    ttGCTTATTTAATCAGTTGGTGCGCGATTGGCGAACGTACTACCGAAATTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCTGTAGGTTTTAAAAATGGTAGCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCa                                                                                                                <  1:403290/149‑1 (MQ=255)
                                       cTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGcc                                                                                                             >  2:469152/1‑149 (MQ=255)
                                                 cAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCgccagcc                                                                                                   <  1:393402/149‑1 (MQ=255)
                                                   gTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAGATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGAGGATTGCTGTGGATGCTATACGCGCAGCCCGCGCCAGCCat                                                                                                 >  2:214710/1‑149 (MQ=255)
                                                   gTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATAAGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCat                                                                                                 >  2:634704/1‑149 (MQ=255)
                                                    ttGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCAACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCata                                                                                                >  2:264245/1‑149 (MQ=255)
                                                     gttttttgccgcaacgcataataatgtggccatacgggttgccgctggtctgttagatggtcatctgaccatttttgtctggcgagaggaacatATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCatat                                                                                               >  2:913765‑M1/95‑149 (MQ=255)
                                                       tttttgccgccacgcataatactgtggccatacgggttgccgctggtctgatagatggtcatctgaccagttttgtctggcgagaggaacatATGGCAATACGCGGATTGCGGGGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGt                                                                                             >  1:630063‑M1/93‑149 (MQ=255)
                                                       gggCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGCATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGt                                                                                             >  1:68000/1‑149 (MQ=255)
                                                         gcgcGATTGGCGCACGTACTACCGAAAGGCAGATCCCCCGCGAAATGGCTTCGGCACTCTCCTGGCCGGTAGGTTTTAAAACTGGTACCGATGGCACTACGCGGATGGCTGTGGAGGCTATCCGCGCAGCCCGCGCCAGCCAGATGtgc                                                                                           >  1:233872/1‑147 (MQ=255)
                                                              aTTGGCGCACGTACTACCGAAAGGCAGCACCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCAGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCtctc                                                                                      <  2:268945/149‑1 (MQ=255)
                                                                  gcgcACGTACTACCGAAAGTCAGATCCACCGCGAAATTGCTTAGGCACTCTCCTGTCCGGTATGTTTTAAAAATGGTACAGATGGCAATAAGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCa                                                                                  <  2:458639/149‑1 (MQ=255)
                                                                   cgcACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCAATCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGAAGCTATCCGCGCAGCCCGCGCCAGCCATATGGTCCTCTCGCCAg                                                                                 <  2:677671/149‑1 (MQ=255)
                                                                    gcccGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCAGATGTTACTGTCGccgcc                                                                                >  1:706362/4‑146 (MQ=255)
                                                                           ctacCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAAt                                                                         <  1:469152/149‑1 (MQ=255)
                                                                              ccGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGCTGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGt                                                                      >  1:297019/1‑149 (MQ=255)
                                                                              ccGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACTCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATgga                                                                      >  2:773806/1‑148 (MQ=255)
                                                                              ccGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATggg                                                                      >  1:33238/1‑148 (MQ=255)
                                                                              ccGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGt                                                                      >  1:301740/1‑149 (MQ=255)
                                                                              ccGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGt                                                                      >  2:131041/1‑149 (MQ=255)
                                                                               cGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTc                                                                     <  2:390025/149‑1 (MQ=255)
                                                                                 aaaGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAg                                                                   <  1:126938/149‑1 (MQ=255)
                                                                                  aaGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                                  <  1:923909/149‑1 (MQ=255)
                                                                                    gTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGATTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATg                                                                >  2:847876/1‑149 (MQ=255)
                                                                                      cacaTCCACCGCGAAAGGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGAc                                                              <  2:431747/146‑1 (MQ=255)
                                                                                        gATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATAAGTTCCTCTCGCCAGACAAAAATGGTCAGATGCCCa                                                            >  1:507341/1‑149 (MQ=255)
                                                                                         aTCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCAt                                                           >  1:392743/1‑149 (MQ=255)
                                                                                             aCCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCCTCTAt                                                       >  1:72087/1‑149 (MQ=255)
                                                                                              ccGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATc                                                      <  2:587078/149‑1 (MQ=255)
                                                                                               cgcgAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGAGGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCa                                                     <  1:981697/149‑1 (MQ=255)
                                                                                                gcgAAATGGCTTCGGCACTCTCCTGTCCGGTCGGTTTTAAAAATGGTACCGATGGCAATTCGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAATAATGGTCAGATGACCATCTATCAg                                                    <  1:971290/149‑1 (MQ=255)
                                                                                                gcgAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAg                                                    >  1:254736/1‑149 (MQ=255)
                                                                                                 cgAAATGGCTTCGGCACTCTCCTGTCCGGAAGTTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATCCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGa                                                   <  1:738547/149‑1 (MQ=255)
                                                                                                      tGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGc                                              <  1:770744/149‑1 (MQ=255)
                                                                                                              gCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAccc                                      >  2:749401/1‑149 (MQ=255)
                                                                                                               cACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGCAATGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCACGCCAGACAAAAATGGTCACATGACCAACTAACAGACCAGCGGCCACCCg                                     >  2:245592/1‑149 (MQ=255)
                                                                                                                  tctcCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                                  <  2:743839/149‑1 (MQ=255)
                                                                                                                       tGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCca                             >  2:284194/1‑149 (MQ=255)
                                                                                                                          ccGGTAGGTTTTAAAAATGGTACCGATGTCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACat                          <  1:887987/149‑1 (MQ=255)
                                                                                                                           cGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCACCCCGTATGGCCACatt                         >  1:99728/1‑149 (MQ=255)
                                                                                                                           cGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCAAGCCATATCTTCCTCTCGCCAGACAAAAATCTTCAGATGCCCATCTATCAGCCCAGCGGCAACCCGTATGGCCACatt                         >  2:801937/1‑149 (MQ=255)
                                                                                                                            ggTAGGTTTTAAAAATGGTCCCGATGGCAATACTCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCTCCAGCCATATGTTCCTCTCGCCAGTCACAAATGGTCAGATGACAATCTATCAGACCAGAGGCACCCCGTATGGCCACatta                        >  1:847052/1‑149 (MQ=255)
                                                                                                                                  ttttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCACACCAGCGGCAACCCGTATGGCCACAATATTATGc                  >  2:617669/1‑149 (MQ=255)
                                                                                                                                   tttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGt                >  1:129986/1‑149 (MQ=255)
                                                                                                                                    tttaaaaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCCACCCGTATGGCCCCAGTATTATGCGt                >  2:325776/4‑149 (MQ=255)
                                                                                                                                        aaaTGGTACCGATGGCAATTCGAGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAATATGGGCAGATGACCATCTATCAGACCAGCTGCAACCCGTATGGCCACATTATTATGCGTggcg            <  1:219413/149‑1 (MQ=255)
                                                                                                                                                ccGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCAACTATCAGACCAGCGGCAACCCGTATGGCCACATTAGTATGCGTGGCGGCaaaaaa    >  1:319279/1‑149 (MQ=255)
                                                                                                                                                  gATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAAcc  <  2:731947/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACC  >  NC_000913/1788853‑1789147

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