Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 3,992,705 IS5 (+) +4 bp 5.3% coding (1553‑1556/2547 nt) cyaA → adenylate cyclase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2288919 =NA (NA)11 (0.050) 11/322 NT 5.4% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 3992705 = 188 (0.940)coding (1553/2547 nt) cyaA adenylate cyclase
* ? NC_000913 = 2290113NA (NA)10 (0.050) 10/322 NT 5.0% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 = 3992708 186 (0.930)coding (1556/2547 nt) cyaA adenylate cyclase
Rejected: Frequency below/above cutoff threshold.

Only 100 of 245 total aligned reads displayed.

TGCTGCATGCAGTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACA  >  NC_000913/2289843‑2290382
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                              
tgctgcATGCAGTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATaa                                                                                                                                                                                                                                                                            >  2:674140/1‑275 (MQ=31)
tgctgcATGCAGTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATa                                                                                                                                                                                                                                                                             >  4:41890/1‑274 (MQ=31)
tgctgcATGCAGTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATAACCGGGACTTTTTCCCCCCTTCCATac                                                                                                                                                                                                                                                                            >  4:463664/1‑274 (MQ=9)
     catgactggAAGGCGTAGCATGGTGTCCAGCGGATAAGGTCGCCCGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAGTCTGCTCCATGCGGGACAAGAAAATCTCGTTTCTGGTCTGACGGCGCGTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACg                                                                                                                                                                                                                                                                          <  4:101261/264‑1 (MQ=255)
            tGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccac                                                                                                                                                                                                                                                                <  2:463664/275‑1 (MQ=35)
                      gCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTAATATCAGGGACTTGTTCGCACCTTCCATACCCCTGTACCCACCAGAACACat                                                                                                                                                                                                                                                      >  4:18164/1‑275 (MQ=255)
                                            ttgCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCGTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGaaa                                                                                                                                                                                                                                <  4:606415/273‑1 (MQ=255)
                                               ccGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCCCCAGAACAGATATTGGGGAACgaaaaagaaaaaacg                                                                                                                                                                                                                             >  4:62506/1‑272 (MQ=255)
                                                          cAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTg                                                                                                                                                                                                                   <  2:91562/274‑1 (MQ=255)
                                                                 ggggTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAACATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACACAGAGAAACAGAAACAGATTGATGcctg                                                                                                                                                                                                            >  2:120389/1‑271 (MQ=255)
                                                                 ggggTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTAGCACCTTCCATAACGCGGTAGCCACCAGAACAGATCTTGCGGAACGACAAAGCGAAACAGAACCAGCTTGCTGCATTg                                                                                                                                                                                                           >  2:716519/1‑275 (MQ=255)
                                                                     tAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTa                                                                                                                                                                                                                                                                     <  4:4427/213‑1 (MQ=34)
                                                                            gCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATc                                                                                                                                                                                                 <  2:488696/274‑1 (MQ=255)
                                                                                           ccaTGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTg                                                                                                                                                                                  <  2:613671/274‑1 (MQ=255)
                                                                                                 tttGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTaa                                                                                                                                                                           <  4:509978/275‑1 (MQ=255)
                                                                                                             gAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTa                                                                                                                                                               <  4:659950/275‑1 (MQ=255)
                                                                                                              aaTCTGCTCCATACGGGACAAGAAAATCACTTTTCTGGTCTTACGCCGCTTACTGATGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTa                                                                                                                                                               <  2:423867/274‑1 (MQ=255)
                                                                                                                    cTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCAATGAGCTTTCATCCTATGAAATTAATTGCCGTTAAAAGCATTGGGTACAGaaaa                                                                                                                                                        >  2:560241/1‑275 (MQ=255)
                                                                                                                        aTGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGAATTGGGTACAGaaactgtcttt                                                                                                                                                 >  1:26422/1‑270 (MQ=255)
                                                                                                                                   gAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTg                                                                                                                                                                                                                                                                       >  5:17393/1‑149 (MQ=35)
                                                                                                                                     aaaTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGt                                                                                                                                                                                                                                                                      <  7:437226/148‑1 (MQ=38)
                                                                                                                                     aaaTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGt                                                                                                                                                                                                                                                                      <  7:878337/148‑1 (MQ=37)
                                                                                                                                     aaaTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTa                                                                                                                                                                                                                                                                     <  5:560372/149‑1 (MQ=38)
                                                                                                                                            ttttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccacca                                                                                                                                                                                                                                                              >  6:740362/1‑149 (MQ=255)
                                                                                                                                            ttttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGAGGAACCCTGTTCTATGGCTCCAGATGAGAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccacca                                                                                                                                                                                                                                                              >  5:988479/1‑149 (MQ=255)
                                                                                                                                               tCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGaa                                                                                                                                                                                                                                                           <  6:656112/149‑1 (MQ=255)
                                                                                                                                                    tCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAaatta                                                                                                                                                                                      >  2:757500/1‑213 (MQ=255)
                                                                                                                                                      tGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACAAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTaa                                                                                                                      >  4:294133/1‑275 (MQ=255)
                                                                                                                                                        aCGGNGCTTNCTGCTGAATTNACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTg                                                                                                                                                                                  <  4:757500/213‑1 (MQ=255)
                                                                                                                                                        aCGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGata                                                                                                                                                                                                                                                     <  5:22874/146‑1 (MQ=255)
                                                                                                                                                         cGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTg                                                                                                                                                                                                                                                  <  7:385092/148‑1 (MQ=255)
                                                                                                                                                             gcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGt                                                                                                               >  2:253547/1‑275 (MQ=255)
                                                                                                                                                              cTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGaa                                                                                                                                                                                                                                             <  7:314010/148‑1 (MQ=255)
                                                                                                                                                              cTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGaa                                                                                                                                                                                                                                             >  7:870300/1‑148 (MQ=255)
                                                                                                                                                              cTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGaa                                                                                                                                                                                                                                             >  7:772683/1‑148 (MQ=255)
                                                                                                                                                               ttactttttaaTTCACTTTCGACGTAAGTAAGNTNGTGACTCCTGTTTAACCCTGTGCTATGGCTCCAGATGACAAACATNGTCTCATATCAGGGACTTTCTCGCACCTTACATAACGCTGTAGCGACCAGAACAGATATTGCGGAACGANAAAGAGAAACAGAACCAGCTTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAg                                                                                                         <  3:50204/270‑1 (MQ=255)
                                                                                                                                                               ttACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACg                                                                                                                                                                                                                                           >  5:866455/1‑149 (MQ=255)
                                                                                                                                                                   tgctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccac                                                                                                                                                                                                                                                                <  5:956063/124‑1 (MQ=255)
                                                                                                                                                                    gctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACaaa                                                                                                                                                                                                                                      <  5:467317/149‑1 (MQ=255)
                                                                                                                                                                     ctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAag                                                                                                                                                                                                                                     >  6:875274/1‑149 (MQ=255)
                                                                                                                                                                     aggAATTCACTGTCGGCGAAGGTAAGTTGATGACTCAGGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTAGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGaa                                                                                                       <  4:716519/273‑1 (MQ=255)
                                                                                                                                                                        aaTTCACTGTCGGCGAAGTTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGaa                                                                                                                                                                                                                                                           <  3:229056/124‑1 (MQ=255)
                                                                                                                                                                          ttCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGaaa                                                                                                                                                                                                                                >  5:106382/1‑149 (MQ=255)
                                                                                                                                                                                    gCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGa                                                                                                                                                                                                                      <  6:292553/149‑1 (MQ=255)
                                                                                                                                                                                         ggTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGa                                                                                                                                                                                                                  <  7:562155/148‑1 (MQ=255)
                                                                                                                                                                                          gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATg                                                                                                                                                                                                                >  6:120307/1‑149 (MQ=255)
                                                                                                                                                                                            aaGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGc                                                                                                                                                                                                               <  8:160802/148‑1 (MQ=255)
                                                                                                                                                                                            aaGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCa                                                                                                                                                                                                              <  6:439035/149‑1 (MQ=255)
                                                                                                                                                                                                tgatgaCTCATGAAGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGAACCTTCCATAACGCTGTAGCCACCAGAACAGAGATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGa                                                                                                                                                                                                          <  6:988479/149‑1 (MQ=255)
                                                                                                                                                                                                     aCTCATGATGAACCCTGTTCTATGGCTCCAGATTACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGtttt                                                                   <  1:57468/279‑1 (MQ=255)
                                                                                                                                                                                                      cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTc                                                                                                                                                                                                    <  5:354208/149‑1 (MQ=255)
                                                                                                                                                                                                      cTCATGATGAAACCTTTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGAACCTTCAATAACGCAGTAGCCACCATAACAGATATTGCGGAACGGCAAAGAGAAACAGAACCAGATTGATGCATTGAGCTATc                                                                                                                                                                                                    <  6:1006999/149‑1 (MQ=25)
                                                                                                                                                                                                       tCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTc                                                                                                                                                                                                    <  6:857513/148‑1 (MQ=255)
                                                                                                                                                                                                        cATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTGGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCa                                                                                                                                                                                                   <  8:451641/148‑1 (MQ=255)
                                                                                                                                                                                                        cATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAgaga                                                                                                                                                                                                                                  <  7:497972/117‑1 (MQ=255)
                                                                                                                                                                                                        cATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAgaga                                                                                                                                                                                                                                  >  8:497972/1‑117 (MQ=255)
                                                                                                                                                                                                        cATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGa                                                                                                                                                                                                                  <  8:618927/133‑1 (MQ=255)
                                                                                                                                                                                                        cATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCa                                                                                                                                                                                                   >  8:371940/1‑148 (MQ=255)
                                                                                                                                                                                                         atgatgAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCAt                                                                                                                                                                                                  <  7:916899/148‑1 (MQ=255)
                                                                                                                                                                                                            atgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTa                                                                                                                                                                                              <  5:680962/149‑1 (MQ=255)
                                                                                                                                                                                                            atgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTa                                                                                                                                                                                              >  6:794102/1‑149 (MQ=255)
                                                                                                                                                                                                             tgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTa                                                                                                                                                                                              >  8:237023/1‑148 (MQ=255)
                                                                                                                                                                                                              gaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAagag                                                                                                                                                                                                                                   <  8:792425/110‑1 (MQ=255)
                                                                                                                                                                                                              gaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAagag                                                                                                                                                                                                                                   >  7:792425/1‑110 (MQ=255)
                                                                                                                                                                                                              gaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCAACTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCtt                                                               <  2:381910/274‑1 (MQ=255)
                                                                                                                                                                                                                aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGa                                                                                                                                                                                                                                                       <  6:972025/88‑1 (MQ=255)
                                                                                                                                                                                                                  ccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACaaa                                                                                                                                                                                                                                      >  6:309994/1‑103 (MQ=255)
                                                                                                                                                                                                                  ccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACaaa                                                                                                                                                                                                                                      <  5:309994/103‑1 (MQ=255)
                                                                                                                                                                                                                  ccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCt                                                                                                                                                                                                       >  8:276749/1‑134 (MQ=255)
                                                                                                                                                                                                                    tGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAaatt                                                                                                                                                                                       >  8:842229/1‑148 (MQ=255)
                                                                                                                                                                                                                     gTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTaaaaa                                                                                                                   >  1:206948/1‑215 (MQ=255)
                                                                                                                                                                                                                     gTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTaaaaa                                                                                                                   <  3:206948/215‑1 (MQ=255)
                                                                                                                                                                                                                           tGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGAt                                                 >  2:444880/1‑275 (MQ=255)
                                                                                                                                                                                                                               tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCt                                                                                                                                                                                >  6:934979/1‑144 (MQ=255)
                                                                                                                                                                                                                                ccAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTg                                                                                                                                                                                                                                                  >  7:142516/1‑77 (MQ=255)
                                                                                                                                                                                                                                 cAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTaaaa                                                                                                                                                                         <  6:106382/149‑1 (MQ=255)
                                                                                                                                                                                                                                     tGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCAt                                                                                                                                                                                                  >  7:239842/1‑120 (MQ=255)
                                                                                                                                                                                                                                       aCAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTg                                                                                                                                                                   >  5:104641/1‑149 (MQ=255)
                                                                                                                                                                                                                                        cAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGCAACAGAACCAGATTGATGCATTGAGCTTTCACCCTATGCAATTAATTGCTGTTNAAAGCATTGGGTACAGAAAATACCCATAGCCACATACCCGGAGTCCGTTTTTAAAAACTGTTTAAAGAAATGCANAAGTATTGTTATTGATTTTTTAGTTGTTTTNNTTAGTTAGAATGTGAATCAAAATTcctttttnt                                >  3:84419/1‑272 (MQ=255)
                                                                                                                                                                                                                                          aaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTggg                                                                                                                                                                 >  7:924374/1‑148 (MQ=255)
                                                                                                                                                                                                                                           aCATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTa                                                                                                                                                               <  5:794102/149‑1 (MQ=255)
                                                                                                                                                                                                                                             aTGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACa                                                                                                                                                             >  6:624645/1‑149 (MQ=255)
                                                                                                                                                                                                                                             aTGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACa                                                                                                                                                             >  6:554558/1‑149 (MQ=255)
                                                                                                                                                                                                                                                aTCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTAc                                                                                                                                                              <  3:120845/145‑1 (MQ=255)
                                                                                                                                                                                                                                                  ctcATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCAt                                                                                                                                                                                                             >  5:907217/1‑96 (MQ=255)
                                                                                                                                                                                                                                                  ctcATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGaaaa                                                                                                                                                        >  5:704665/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    cATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAAt                                                                                                                                                       >  8:528982/1‑148 (MQ=255)
                                                                                                                                                                                                                                                    cATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATa                                                                                                                                                      <  5:100945/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    cATATCAGGGACTTGTTCGCACCATCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAAt                                                                                                                                                       >  7:63833/1‑148 (MQ=255)
                                                                                                                                                                                                                                                       atCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATAccc                                                                                                                                                   >  6:272772/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           gggACTTGTTCGCACCTTCCATAACGCTGTAGCCCCCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTNAAAGCATNGGGNACAGAANATACCCATAGCTCAATACACGGAGTCAGTTTTTaaaaa                                                                                                                   >  3:33022/1‑177 (MQ=255)
                                                                                                                                                                                                                                                           gggACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAg                                                                                                                                               <  5:120307/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           gggACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTaaaaa                                                                                                                   <  1:33022/177‑1 (MQ=255)
                                                                                                                                                                                                                                                            tgACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGACCAGAGCTTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCGAGGCAAGTAATTCCTATTAAAAGCATTTGGTACAAAAAATACCCATAGCTCCATGCCAGGAGTCGCTTTTTGAAAACTGTTTAAATAAATGCCCAAGTATTCTGATTGATGTTTTAGTTGTTTTTCTTCATGTGAAGCTGATGAAAAATTCCTCTTTTTAAATAAAAATGATGcc               >  4:405104/2‑275 (MQ=255)
                                                                                                                                                                                                                                                            ggACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAg           >  1:100382/1‑279 (MQ=255)
                                                                                                                                                                                                                                                               cTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGa        >  1:132084/1‑279 (MQ=255)
                                                                                                                                                                                                                                                                ttGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTc                                                                                                                                            >  5:894290/1‑147 (MQ=255)
                                                                                                                                                                                                                                                                       cACCTTCCATAACGCTGTAGCAACCAGAACAGATATTGCAGAACGAAAAAGAAAAACAGAAACAGATTGATGCATTGAGCTATCATCATATGAAATTAATTGCTGTTAAAAGCAGTGGGTACAGAAAATACCCATAGCTCCATAGCCGGAGTCAGTTTTTAAAAACTGTTTAAAGGAATGCACAAGTATTGTGATTGATTTTTTGGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACa  <  1:84419/277‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ccTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCAtt                                                                                                                                                                    >  6:898777/1‑114 (MQ=255)
                                                                                                                                                                                                                                                                            tCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTc                                                                                                                               >  8:1058559/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                              
TGCTGCATGCAGTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACA  >  NC_000913/2289843‑2290382

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