Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 3,922,629 IS5 (–) +4 bp 15.1% intergenic (‑189/+425) atpI ← / ← rsmG ATP synthase, membrane‑bound accessory factor/16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3652036 =NA (NA)24 (0.130) 24/298 NT 14.2% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 3922632 141 (0.760)intergenic (‑192/+425) atpI/rsmG ATP synthase, membrane‑bound accessory factor/16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein
* ? NC_000913 = 3653230NA (NA)26 (0.140) 23/298 NT 15.2% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 3922629 = 141 (0.760)intergenic (‑189/+428) atpI/rsmG ATP synthase, membrane‑bound accessory factor/16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein

Only 100 of 205 total aligned reads displayed.

AACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGAC  >  NC_000913/3922355‑3922889
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                    
caatgctgcatgcagtgaatgcgtagcatggtttccagcggataaggtcgccggccattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccc                                                                                                                                                                                                                                                                      <  4:390929‑M2/1‑1 (MQ=255)
                                                                                       aCGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTcg                                                                                                                                                                                >  2:446730/1‑274 (MQ=255)
                                                                                          tccaccatgttttgccagggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgcggaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGca                                                                                                                                                                            <  2:469047‑M2/91‑1 (MQ=255)
                                                                                                 cTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCa                                                                                                                                                                       <  4:118432/273‑1 (MQ=255)
                                                                                                      gCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACCCTGGTTATAAAGCGTGCTTCAGATCACATATTGCCCATGTTCCCCCACAGCAtatttat                                                                                                                                                                >  4:53821/1‑275 (MQ=255)
                                                                                                                           atgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTc                                                                                                                                                                                                       <  2:394808‑M2/64‑1 (MQ=255)
                                                                                                                                 tataAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTaa                                                                                                                                                                                                                                                                   <  6:250123/149‑1 (MQ=255)
                                                                                                                                   taAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAAgg                                                                                                                                                                                                                                                                 <  5:487991‑M2/149‑3 (MQ=255)
                                                                                                                                   taAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAAcc                                                                                                                                                                                                                                                                 >  5:844795/1‑149 (MQ=255)
                                                                                                                                       aCTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAAggaagg                                                                                                                                                                                                                                                             <  7:32767‑M2/149‑7 (MQ=255)
                                                                                                                                       aCTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGt                                                                                                                                                                                                                                                             >  5:253036/1‑149 (MQ=255)
                                                                                                                                            ctttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATCAAATTGCGACTTTTCTGCCAACCCTGTTCGAt                                                                                                                                                                                                                                                        <  6:654792/149‑1 (MQ=255)
                                                                                                                                             tttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATc                                                                                                                                                                                                                                                       >  5:249894/1‑149 (MQ=255)
                                                                                                                                               tAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCAt                                                                                                                                                                                                                                                     <  7:1051299/149‑1 (MQ=255)
                                                                                                                                                    cTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAAggaaggtgcgaataagcgg                                                                                                                                                                                                                                                >  6:487991‑M2/1‑130 (MQ=255)
                                                                                                                                                     ttAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGaaaaa                                                                                                                                                                                                                                               >  7:240068/1‑149 (MQ=255)
                                                                                                                                                             aaGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGGTTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAAggaaggtgcgaataagcggggaaattcttc                                                                                                                                                                                                                                     <  2:716511‑M2/151‑31 (MQ=255)
                                                                                                                                                               gCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATaa                                                                                                                                                                                                                                     <  7:157307/149‑1 (MQ=255)
                                                                                                                                                                 ttCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAacac                                                                                                                                                                                                                          <  4:469646/157‑1 (MQ=255)
                                                                                                                                                                 aTCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGGTTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAAggaaggtgcgaataagcggggaaattcttctcgg                                                                                                                                                                                                                                 >  4:716511‑M2/1‑117 (MQ=255)
                                                                                                                                                                  ttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAAcgc                                                                                                                                                                                                                                  <  7:880151‑M2/37‑1 (MQ=255)
                                                                                                                                                                  tCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAAcgc                                                                                                                                                                                                                                  >  5:810485/1‑149 (MQ=255)
                                                                                                                                                                    aaaaaGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTACTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGt                                                                                                                                                                                                                                >  6:894092/1‑149 (MQ=255)
                                                                                                                                                                          cGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCTAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAacac                                                                                                                                                                                                                          <  6:1092686/149‑1 (MQ=255)
                                                                                                                                                                            gTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTg                                                                                                                                                                                                                        <  7:371362/149‑1 (MQ=255)
                                                                                                                                                                               aaaTTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTg                                                                                                                                                                                                                     >  7:841515/1‑149 (MQ=255)
                                                                                                                                                                                   tataCGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGGTCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATaa                                                                                                                                                                                                                 <  8:625176/149‑1 (MQ=255)
                                                                                                                                                                                      aCGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGc                                                                                                                                                                                                              <  8:483506/149‑1 (MQ=255)
                                                                                                                                                                                       cGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACa                                                                                                            >  4:160725/1‑246 (MQ=255)
                                                                                                                                                                                          tGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgacatttaaggagttatccccagtttttagtgagatctctcccactgacgtctcatttggtcc                                                                            <  2:276547‑M2/275‑184 (MQ=255)
                                                                                                                                                                                            cgcCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTaa                                                                                                                                                                                                                              >  2:50389/1‑127 (MQ=255)
                                                                                                                                                                                             gcCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTc                                                                                                                                                                                                       <  8:838244/149‑1 (MQ=255)
                                                                                                                                                                                              cccccGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCa                                                                                                                                                                                                      <  7:45577/149‑1 (MQ=255)
                                                                                                                                                                                                agttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGa                                                                                                                                                                                                    <  7:308829‑M2/67‑1 (MQ=255)
                                                                                                                                                                                                 gttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGAt                                                                                                                                                                                                   <  6:343504‑M2/68‑1 (MQ=255)
                                                                                                                                                                                                 ccGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgttt                                                                                                                                                                                                   >  8:579189‑M2/1‑85 (MQ=255)
                                                                                                                                                                                                 ccGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGAt                                                                                                                                                                                                   <  6:810485/149‑1 (MQ=255)
                                                                                                                                                                                                 ccGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGAt                                                                                                                                                                                                   >  5:1062400/1‑149 (MQ=255)
                                                                                                                                                                                                    tGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATcac                                                                                                                                                                                                <  5:777474/149‑1 (MQ=255)
                                                                                                                                                                                                     gATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagc                                                                                                                                                                                               <  5:621715‑M2/149‑69 (MQ=255)
                                                                                                                                                                                                        ttCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAg                                                                                                                                                                                                     <  4:39983/140‑1 (MQ=255)
                                                                                                                                                                                                          tcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTgg                                                                                                                                                                                                                       <  8:387497‑M2/48‑1 (MQ=255)
                                                                                                                                                                                                          tcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTgg                                                                                                                                                                                                                       >  7:387497‑M2/73‑120 (MQ=255)
                                                                                                                                                                                                          tcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTg                                                                                                                                                                                          >  8:218214‑M2/73‑149 (MQ=255)
                                                                                                                                                                                                            aaCAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTgcg                                                                                                                                                                                        <  6:454901/149‑1 (MQ=255)
                                                                                                                                                                                                               aaTAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttc                                                                                                                                                                                     >  6:754663‑M2/1‑71 (MQ=255)
                                                                                                                                                                                                               aaTAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAAggaaggtgcgaataagcggggaaattcttcggctgactcagtcatttcatttcttcatgtttgagccgattttttctc                                                                                                                                                                                     <  8:593827‑M2/149‑79 (MQ=255)
                                                                                                                                                                                                               aaTAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCAt                                                                                                                                                                                     >  6:868258/1‑149 (MQ=255)
                                                                                                                                                                                                               aaTAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCAt                                                                                                                                                                                     >  6:974166/1‑149 (MQ=255)
                                                                                                                                                                                                               aaTAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCAt                                                                                                                                                                                     <  7:1008302/149‑1 (MQ=255)
                                                                                                                                                                                                               aaTAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCAt                                                                                                                                                                                     >  8:117397/1‑149 (MQ=255)
                                                                                                                                                                                                                   accctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGca                                                                                                                                                                            <  4:500252‑M2/91‑1 (MQ=255)
                                                                                                                                                                                                                    gTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTcg                                                                                                                                                                                >  8:885710/1‑149 (MQ=255)
                                                                                                                                                                                                                     tAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTcgc                                                                                                                                                                               <  6:632560/149‑1 (MQ=255)
                                                                                                                                                                                                                     cctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCtt                                                                                                                                                                                                        >  6:535671‑M2/62‑124 (MQ=255)
                                                                                                                                                                                                                     cctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCtt                                                                                                                                                                                                        <  5:535671‑M2/63‑1 (MQ=255)
                                                                                                                                                                                                                         cAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCTCATATTGCGCATGTTCGCGcac                                                                                                                                                                           >  5:327370/1‑149 (MQ=255)
                                                                                                                                                                                                                          aaaaGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGcacg                                                                                                                                                                          >  8:162791/1‑148 (MQ=255)
                                                                                                                                                                                                                          aaaaGGTGAATAAATGTTTAAATATCTTTCCAGTGCATACAATTGCGACTTTTCTGCTAAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatg                                                                                                                                                                          <  7:576879‑M2/149‑90 (MQ=255)
                                                                                                                                                                                                                             atggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCa                                                                                                                                                                       >  5:186217‑M2/54‑149 (MQ=255)
                                                                                                                                                                                                                                tGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAAggaaggtgcgaataagcggggaaatt                                                                                                                                                                                                                                         >  6:995420‑M2/1‑54 (MQ=255)
                                                                                                                                                                                                                                tGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCatat                                                                                                                                                                    <  5:451777/149‑1 (MQ=255)
                                                                                                                                                                                                                                     aaaTGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCAtatttatt                                                                                                                                                               >  6:285124/1‑149 (MQ=255)
                                                                                                                                                                                                                                      aaTGTTTAAATA‑TTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCa                                                                                                                                                                       <  1:140900/139‑1 (MQ=255)
                                                                                                                                                                                                                                        tGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTAc                                                                                                                                                            >  6:143763/1‑149 (MQ=255)
                                                                                                                                                                                                                                         gTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACt                                                                                                                                                           <  6:101860/149‑1 (MQ=255)
                                                                                                                                                                                                                                           ttAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATGTTTATTTACTTg                                                                                                                                                         >  5:871948/1‑149 (MQ=255)
                                                                                                                                                                                                                                           ttAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTg                                                                                                                                                         <  7:999940/149‑1 (MQ=255)
                                                                                                                                                                                                                                              aaTATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCa                                                                                                                                                      >  6:195307/1‑149 (MQ=255)
                                                                                                                                                                                                                                              aaTATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGc                                                                                                                                        <  2:165152/163‑1 (MQ=255)
                                                                                                                                                                                                                                                   ttttCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTaa                                                                                                                                                                                                                                         >  4:670771/1‑61 (MQ=255)
                                                                                                                                                                                                                                                     ttCCAGTGCATACAATTGCGACTTTTCTGCTAAggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccatttaaggcgttatccccagtttttagtgagatctctcccactgacgtatcatttggtccgcccgaaacaggttggccagcgtgaataacatcgccagttggttatcgtttttcagcaa                 >  2:612372‑M2/1‑33 (MQ=255)
                                                                                                                                                                                                                                                       ccAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATaa                                                                                                                                                                                                                                     <  2:670771/61‑1 (MQ=255)
                                                                                                                                                                                                                                                          gTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACa                                                                                                            >  2:65322/1‑179 (MQ=255)
                                                                                                                                                                                                                                                           tGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACt                                                                                                      >  4:590355/1‑184 (MQ=255)
                                                                                                                                                                                                                                                            gCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTcg                                                                                                                                                                                >  4:169055/1‑109 (MQ=255)
                                                                                                                                                                                                                                                             cATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCa                                                                                                                                                                                                                                                      >  3:44691/1‑38 (MQ=255)
                                                                                                                                                                                                                                                             agggacttgttcgcaccttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCa                                                                                                                                       >  5:329210‑M2/22‑149 (MQ=255)
                                                                                                                                                                                                                                                              aTACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATa                                                                                                        <  4:65322/179‑1 (MQ=255)
                                                                                                                                                                                                                                                               tACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGttt                                                                                                  <  2:590355/184‑1 (MQ=255)
                                                                                                                                                                                                                                                               ggacttgttcgcaccttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAat                                                                                                                                                   >  5:672722‑M2/20‑135 (MQ=255)
                                                                                                                                                                                                                                                               ggacttgttcgcaccttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAat                                                                                                                                                   <  6:672722‑M2/116‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aTTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTa                                                                                                                                 <  5:933998/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aTTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCATTGCAAGTTTa                                                                                                                                 <  7:611993/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ttGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTAt                                                                                                                                >  5:426883/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    ttGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTAt                                                                                                                                <  8:36371/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ttGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTACAAGTTTAt                                                                                                                                >  7:467040/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      ttcgcaccttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATAAATGa                                                                                        <  1:91755‑M2/175‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTaa                                                                                                                                                                                                                              >  2:711239/1‑52 (MQ=255)
                                                                                                                                                                                                                                                                       cGACTTTTATGCTAAggcaggtgcgaataagcggggaaattcttctcggctgactcagtaatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagactatgtagaccgtttcttcgccatttaaggcgttatccccagtttttagggagatctctcccactgacgtatcatatggtccgcgcgaaacaggttggccagcgtgaataacatctacagttggttatcgtttttcagcaaccacttgtatctggc  <  2:104403‑M2/272‑258 (MQ=255)
                                                                                                                                                                                                                                                                        gtcTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTGTATGata                                                                                                                            >  7:929431/3‑149 (MQ=255)
                                                                                                                                                                                                                                                                         aCTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGatat                                                                                                                           >  8:78843/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          cTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATAtt                                                                                                                          <  7:261807/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ttttCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAacac                                                                                                                                                                                                                          <  4:711239/52‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ttttCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATAttt                                                                                                                         <  5:967898/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                            tttCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATAAATGACAt                                                                                     >  2:500213/1‑184 (MQ=255)
                                                                                                                                                                                                                                                                             ttCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCa                                                                                                                       <  8:240068/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCa                                                                                                                       >  6:404960‑M2/6‑149 (MQ=255)
                                                                                                                                                                                                                                                                              ttccCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAg                                                                                                                      >  8:310086‑M2/5‑149 (MQ=255)
                                                                                                                                                                                                                                                                                  cTAACCCTTTTAGAGCATGAAAAACTGTAAATAACGCGTAAACACTTGTGATAAAGCGTGCTTCAGTTCACATATTTCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGACATTTCAGTCTAAAAACAGATACTGTTTGAATAAATGACATTTACACAACAAACACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCACCGAt             <  4:447828/250‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                    
AACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAGCGTGCTTCAGATCACATATTGCGCATGTTCGCGCACAGCATATTTATTTACTTGGCAAATGATGCCTTTGCAAGTTTATGATATTTCAGTCTAAAAACAGATACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGAC  >  NC_000913/3922355‑3922889

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