Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,857,096 Δ2,538 bp coding (4‑2541/2562 nt) mutS → DNA mismatch repair protein MutS

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2857096 2859633 2538 74 [0] [0] 75 mutS DNA mismatch repair protein MutS

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 28570950 (0.000)74 (1.140) 56/276 0.1 100% coding (3/2562 nt) mutS DNA mismatch repair protein MutS
?NC_000913 2859634 = 0 (0.000)coding (2542/2562 nt) mutS DNA mismatch repair protein MutS

GCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACG  >  NC_000913/2859502‑2859773
                                                                                                                                    |                                                                                                                                           
catcaagaactcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGaa                                                                                                                                      <  2:783746‑M2/8‑1 (MQ=255)
  tcaagaactcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAaga                                                                                                                                    <  1:1093354‑M2/10‑1 (MQ=255)
    aagaactcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGc                                                                                                                                  <  2:616253‑M2/12‑1 (MQ=255)
     agaactcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGcc                                                                                                                                 <  2:1070962‑M2/13‑1 (MQ=255)
      gaactcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCt                                                                                                                                >  2:1038178‑M2/127‑140 (MQ=255)
      gaactcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCt                                                                                                                                <  2:1063142‑M2/14‑1 (MQ=255)
      gaactcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCt                                                                                                                                <  2:891025‑M2/14‑1 (MQ=255)
        actcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCTTgg                                                                                                                              <  2:759263‑M2/16‑1 (MQ=255)
          tcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGgtg                                                                                                                            >  2:918237‑M2/123‑140 (MQ=255)
          tcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGgtg                                                                                                                            <  1:331668‑M2/18‑1 (MQ=255)
           cgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGgtgt                                                                                                                           >  1:421509‑M2/122‑140 (MQ=255)
            gccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGta                                                                                                                          >  2:319747‑M2/121‑140 (MQ=255)
            gccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGta                                                                                                                          >  1:578566‑M2/121‑140 (MQ=255)
                ttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGtaataa                                                                                                                      <  2:578566‑M2/24‑1 (MQ=255)
                    cttcttcccctgaaatgattaactccggtatcatgtgagccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggtcataacctcatgCGCTTGAAGAGCCTGGTGTAATGACAAt                                                                                                                  >  2:973392‑M2/113‑140 (MQ=255)
                         tcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCg                                                                                                             <  2:800516‑M2/33‑1 (MQ=255)
                          cccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa                                                                                                            >  1:462394‑M2/107‑140 (MQ=255)
                             ctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAg                                                                                                         >  1:86941‑M2/104‑140 (MQ=255)
                               gaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTc                                                                                                       >  1:50748‑M2/102‑140 (MQ=255)
                                 aatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCtt                                                                                                     <  2:1060803‑M2/41‑1 (MQ=255)
                                    gattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTg                                                                                                  <  1:1038178‑M2/44‑1 (MQ=255)
                                      ttaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCt                                                                                                <  1:455029‑M2/46‑1 (MQ=255)
                                       taactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTa                                                                                               >  2:926419‑M2/94‑140 (MQ=255)
                                        aactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTAt                                                                                              >  1:487899‑M2/93‑140 (MQ=255)
                                           tccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCgg                                                                                           >  2:497219‑M2/90‑140 (MQ=255)
                                                tatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAata                                                                                      <  2:648568‑M2/56‑1 (MQ=255)
                                                     tgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAAc                                                                                 >  2:350766‑M2/80‑140 (MQ=255)
                                                         cgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACa                                                                             >  2:398053‑M2/76‑140 (MQ=255)
                                                            cttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACt                                                                          >  2:1087993‑M2/73‑140 (MQ=255)
                                                                    attacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaata                                                                  >  2:107923‑M2/65‑140 (MQ=255)
                                                                        caacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaataaaat                                                              >  1:623559‑M2/61‑140 (MQ=255)
                                                                         aacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaataaaata                                                             >  2:648208‑M2/60‑140 (MQ=255)
                                                                            gaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAAtaaaataaaa                                                          <  1:350766‑M2/84‑1 (MQ=255)
                                                                             aaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATaaaaa                                                         >  1:140133‑M2/56‑140 (MQ=255)
                                                                               aataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAaca                                                       <  1:973392‑M2/87‑1 (MQ=255)
                                                                               aataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAaca                                                       <  2:50748‑M2/87‑1 (MQ=255)
                                                                               aataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAaca                                                       <  2:853815‑M2/87‑1 (MQ=255)
                                                                               aataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAaca                                                       <  1:766428‑M2/87‑1 (MQ=255)
                                                                                 taaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAcc                                                     <  2:86941‑M2/89‑1 (MQ=255)
                                                                                 taaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAcc                                                     <  1:648208‑M2/89‑1 (MQ=255)
                                                                                 taaaaaccatcacaccccatttaatatcagggaaccggacataacaccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAcc                                                     <  1:918237‑M2/89‑1 (MQ=255)
                                                                                  aaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAccc                                                    <  1:319747‑M2/90‑1 (MQ=255)
                                                                                    aaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTg                                                  <  1:107923‑M2/92‑1 (MQ=255)
                                                                                              accccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGa                                        <  2:487899‑M2/102‑1 (MQ=255)
                                                                                                   atttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCt                                   >  1:13715‑M2/34‑140 (MQ=255)
                                                                                                        atatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTAttt                              >  1:113961‑M2/29‑140 (MQ=255)
                                                                                                         tatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTAtttt                             >  2:899577‑M2/28‑140 (MQ=255)
                                                                                                          atcagggaaccggacataaccccatgCGCTTGAAGATCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTa                            <  1:926419‑M2/114‑1 (MQ=255)
                                                                                                          atcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGaaa                                      >  1:745461‑M2/27‑130 (MQ=255)
                                                                                                          atcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGaaa                                      <  2:745461‑M2/104‑1 (MQ=255)
                                                                                                              gggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAgg                        >  2:38655‑M2/23‑140 (MQ=255)
                                                                                                              gggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAgg                        >  1:919586‑M2/23‑140 (MQ=255)
                                                                                                              gggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAgg                        >  2:51427‑M2/23‑140 (MQ=255)
                                                                                                               ggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGa                                        <  2:296964‑M2/102‑1 (MQ=255)
                                                                                                               ggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGa                                        >  1:296964‑M2/22‑123 (MQ=255)
                                                                                                               ggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAggg                       >  2:953263‑M2/22‑140 (MQ=255)
                                                                                                                   ccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTaaa                   >  1:134474‑M2/18‑140 (MQ=255)
                                                                                                                   ccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTaaa                   >  2:1014062‑M2/18‑140 (MQ=255)
                                                                                                                      gacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAAcc                >  2:503730‑M2/15‑140 (MQ=255)
                                                                                                                         ataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATg             >  2:74681‑M2/12‑140 (MQ=255)
                                                                                                                         ataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATg             <  2:747320‑M2/129‑1 (MQ=255)
                                                                                                                         ataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATg             <  2:266357‑M2/129‑1 (MQ=255)
                                                                                                                         ataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATg             <  2:919586‑M2/129‑1 (MQ=255)
                                                                                                                          taaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGc            <  1:1014062‑M2/130‑1 (MQ=255)
                                                                                                                           aaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGcc           >  1:488887‑M2/10‑140 (MQ=255)
                                                                                                                           aaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGcc           <  2:13715‑M2/131‑1 (MQ=255)
                                                                                                                           aaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGcc           <  1:74681‑M2/131‑1 (MQ=255)
                                                                                                                              cccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATc        >  2:439473‑M2/7‑140 (MQ=255)
                                                                                                                                catgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTa      >  2:825418‑M2/5‑140 (MQ=255)
                                                                                                                                  tgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACTCCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGc            <  1:348432‑M2/130‑1 (MQ=255)
                                                                                                                                  tgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGc            >  2:348432‑M2/3‑132 (MQ=255)
                                                                                                                                  tgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTaca    <  2:488887‑M2/138‑1 (MQ=255)
                                                                                                                                  tgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTaca    <  2:134474‑M2/138‑1 (MQ=255)
                                                                                                                                  tgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTaca    >  1:1022947‑M2/3‑140 (MQ=255)
                                                                                                                                    cGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACg  <  2:113961/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                           
GCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACG  >  NC_000913/2859502‑2859773

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑