Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 4,542,238 IS1 (–) +8 bp 5.0% coding (202‑209/597 nt) fimE → regulator for fimA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1978503 =NA (NA)9 (0.080) 9/176 NT 6.2% noncoding (768/768 nt) IS1 repeat region
?NC_000913 = 4542245 123 (1.170)coding (209/597 nt) fimE regulator for fimA
* ? NC_000913 3583428 =NA (NA)4 (0.040) 4/174 NT 2.9% noncoding (1/768 nt) IS1 repeat region
?NC_000913 4542238 = 122 (1.180)coding (202/597 nt) fimE regulator for fimA
Rejected: Frequency below/above cutoff threshold.

Only 100 of 156 total aligned reads displayed.

AATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACT  >  NC_000913/3583330‑3583524
                                                                                                  |                                                                                                
aaTTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGAGATTATAGCATTTCAATATTTACTTATGACATTAAAg                                                                                                  <  1:607209/99‑1 (MQ=255)
 aTTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAgg                                                                                                 <  1:2096982/99‑1 (MQ=255)
 aTTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAgg                                                                                                 >  1:1371044/1‑99 (MQ=255)
 aTTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAgg                                                                                                 >  2:2595166/1‑99 (MQ=255)
  ttttATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGt                                                                                                >  1:1687934/1‑99 (MQ=255)
      ataGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGAtg                                                                                            >  1:2030376/1‑99 (MQ=255)
      ataGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGAtg                                                                                            >  2:2820286/1‑99 (MQ=255)
       taGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGAtgc                                                                                           >  2:1338038/1‑99 (MQ=255)
         gAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGAtgctg                                                                                         >  2:44948/1‑99 (MQ=255)
         gAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGAtgctg                                                                                         >  1:1900713/1‑99 (MQ=255)
         gAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGAtgctg                                                                                         >  2:2749694/1‑99 (MQ=255)
          aaaTAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGAtgctgc                                                                                        >  1:651481/1‑99 (MQ=255)
          aaaTAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGAtgctgc                                                                                        <  1:325733/99‑1 (MQ=255)
           aataatGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGcc                                                                                       >  2:1990593/1‑99 (MQ=255)
           aataatGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGcc                                                                                       <  2:2671124/99‑1 (MQ=255)
               atgatgATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACt                                                                                   <  2:2480024/99‑1 (MQ=255)
               atgatgATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACt                                                                                   <  2:3037785/99‑1 (MQ=255)
                tgatgaTTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACtt                                                                                  >  1:1483881/1‑99 (MQ=255)
                  atgatTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTAc                                                                                >  1:3045575/1‑99 (MQ=255)
                  atgatTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTAc                                                                                <  1:2652273/99‑1 (MQ=255)
                   tgatTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACt                                                                               <  1:1880606/99‑1 (MQ=255)
                     atTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGa                                                                             >  1:1175364/1‑99 (MQ=255)
                      tttCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAt                                                                            >  2:2569651/1‑99 (MQ=255)
                      tttCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAt                                                                            <  1:1660346/99‑1 (MQ=255)
                       ttCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAtt                                                                           <  2:972236/99‑1 (MQ=255)
                       ttCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAtt                                                                           >  1:1253898/1‑99 (MQ=255)
                       ttCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAtt                                                                           <  2:2177545/99‑1 (MQ=255)
                       ttCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAtt                                                                           <  2:1692532/99‑1 (MQ=255)
                       ttCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAtt                                                                           <  2:1544049/99‑1 (MQ=255)
                        tCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAttt                                                                          >  2:1707173/1‑99 (MQ=255)
                         cATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTa                                                                         >  2:1227946/1‑99 (MQ=255)
                          aTAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAg                                                                        <  2:2098217/99‑1 (MQ=255)
                           tAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAgt                                                                       >  1:2448507/1‑99 (MQ=255)
                           tAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAgt                                                                       >  2:2985929/1‑99 (MQ=255)
                            aaaCCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAgtg                                                                      <  1:2476231/99‑1 (MQ=255)
                             aaCCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAgtgt                                                                     >  1:2967749/1‑99 (MQ=255)
                             aaCCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAgtgt                                                                     >  2:2561714/1‑99 (MQ=255)
                                ccTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatg                                                                  >  1:1932164/1‑99 (MQ=255)
                                ccTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatg                                                                  >  1:911853/1‑99 (MQ=255)
                                ccTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatg                                                                  >  2:2729054/1‑99 (MQ=255)
                                 cTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatga                                                                 <  1:3012716/99‑1 (MQ=255)
                                 cTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatga                                                                 >  2:2327119/1‑99 (MQ=255)
                                 cTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatga                                                                 >  1:720771/1‑99 (MQ=255)
                                    aTCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATgg                                                              >  2:896857/1‑99 (MQ=255)
                                    aTCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATgg                                                              >  2:3061722/1‑99 (MQ=255)
                                       tACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGgtgt                                                           <  2:1976639/99‑1 (MQ=255)
                                       tACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGgtgt                                                           <  2:1810640/99‑1 (MQ=255)
                                       tACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGgtgt                                                           <  2:1433113/99‑1 (MQ=255)
                                        aCAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGtt                                                          <  2:1170252/99‑1 (MQ=255)
                                         cAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGttt                                                         >  1:2750458/1‑99 (MQ=255)
                                            gAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTg                                                      >  1:2939529/1‑99 (MQ=255)
                                            gAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTg                                                      >  1:2679343/1‑99 (MQ=255)
                                             aaCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGa                                                     >  2:1442756/1‑99 (MQ=255)
                                             aaCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGa                                                     <  2:679650/99‑1 (MQ=255)
                                             aaCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGa                                                     <  2:2410753/99‑1 (MQ=255)
                                              aCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAg                                                    >  2:876454/1‑99 (MQ=255)
                                               cGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAgg                                                   >  2:407897/1‑99 (MQ=255)
                                                  ttAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGc                                                >  2:2134459/1‑99 (MQ=255)
                                                      tGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCa                                            <  1:260115/99‑1 (MQ=255)
                                                       gATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAg                                           <  1:1263258/99‑1 (MQ=255)
                                                       gATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAg                                           <  1:338141/99‑1 (MQ=255)
                                                        aTCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGt                                          <  1:98461/99‑1 (MQ=255)
                                                        aTCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGt                                          >  1:2122406/1‑99 (MQ=255)
                                                         tCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTg                                         >  1:241565/1‑99 (MQ=255)
                                                           ggggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGc                                       <  1:1812588/99‑1 (MQ=255)
                                                              gATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTc                                    >  2:260124/1‑99 (MQ=255)
                                                              gATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTc                                    >  2:2452297/1‑99 (MQ=255)
                                                               aTTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCt                                   >  1:2532184/1‑99 (MQ=255)
                                                                 tataGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGt                                 >  2:787723/1‑99 (MQ=255)
                                                                   taGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGttt                               <  2:1917976/99‑1 (MQ=255)
                                                                   taGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGttt                               >  2:1520163/1‑99 (MQ=255)
                                                                   taGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGttt                               >  2:2316357/1‑99 (MQ=255)
                                                                    aGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTc                              <  2:1871157/99‑1 (MQ=255)
                                                                       aTTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTAt                           <  1:1583834/99‑1 (MQ=255)
                                                                         ttCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCa                         <  1:1655676/99‑1 (MQ=255)
                                                                          tCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAg                        <  1:97902/99‑1 (MQ=255)
                                                                            aaTATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCt                      >  2:2199623/1‑99 (MQ=255)
                                                                             atatTTACTTATGACATCAAAGGTGATGCTGCCAACTTACTGAGATTGTGTATGATGGTTTATTTGAGGTGATCCAGTGCCTTCTGTGTCTCGCAGCTg                     >  2:2925812/1‑99 (MQ=17)
                                                                              tatTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGt                    >  2:794907/1‑99 (MQ=255)
                                                                              tatTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGt                    >  1:3050722/1‑99 (MQ=255)
                                                                              tatTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGt                    >  1:1372585/1‑99 (MQ=255)
                                                                              tatTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGt                    >  2:178537/1‑99 (MQ=255)
                                                                                tttACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTcc                  >  2:2858384/1‑99 (MQ=255)
                                                                                 ttACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTccc                 >  1:2512662/1‑99 (MQ=255)
                                                                                    cTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctcc              <  2:436805/99‑1 (MQ=255)
                                                                                     ttATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctcct             >  2:2904716/1‑99 (MQ=255)
                                                                                     ttATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctcct             <  2:1360510/99‑1 (MQ=255)
                                                                                        tGACATTAAAGGTGATGCTGCCAACTTACTGATTTGGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGtt          >  1:2812091/1‑99 (MQ=255)
                                                                                        tGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGtt          >  2:584541/1‑99 (MQ=40)
                                                                                          aCATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCa        >  1:2879047/1‑99 (MQ=38)
                                                                                           cATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCat       >  2:1613577/1‑98 (MQ=37)
                                                                                           cATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAg       >  1:1590747/1‑99 (MQ=37)
                                                                                            aTTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGc      <  2:2953741/99‑1 (MQ=35)
                                                                                             ttAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCt     <  1:1476794/99‑1 (MQ=34)
                                                                                              tAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTa    <  2:216209/99‑1 (MQ=33)
                                                                                               aaaGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTAc   <  1:1767048/99‑1 (MQ=32)
                                                                                               aaaGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTAc   <  1:1068459/99‑1 (MQ=32)
                                                                                               aaaGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTAc   <  1:859958/99‑1 (MQ=32)
                                                                                                aaGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACt  <  1:1775920/99‑1 (MQ=31)
                                                                                                aaGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACt  >  2:1160196/1‑99 (MQ=31)
                                                                                                  |                                                                                                
AATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACT  >  NC_000913/3583330‑3583524

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