Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,018,356 (GTC)4→3 coding (2297‑2299/3099 nt) ygfK → putative Fe‑S subunit oxidoreductase subunit

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3018347 3018349 3 88 [0] [0] 87 ygfK putative Fe‑S subunit oxidoreductase subunit

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3018350 =0 (0.000)80 (1.140) 68/280 0.1 100% coding (2291/3099 nt) ygfK putative Fe‑S subunit oxidoreductase subunit
?NC_000913 = 3018346 0 (0.000)coding (2287/3099 nt) ygfK putative Fe‑S subunit oxidoreductase subunit

TACAAGAGATGCCCGCATGGCGCGAAGAGTATGAAGAAGCGTTGCACGACGGCGTAGAGTTCCGTTTCCTGAATAATCCGGAACGTTTCGATGCTGATGGCACCTTAACCTTGCGCGTTATGTCGCTTGGCGAACCGGATGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGG  >  NC_000913/3018200‑3018497
                                                                                                                                                      |                                                                                                                                                   
cactacaagagatgcccgcatggcgcgaagagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagg                                                                                                                                                     <  2:1003674‑M1/1‑1 (MQ=255)
 actacaagagatgcccgcatggcgcgaagagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaaggt                                                                                                                                                    >  2:1116168‑M1/150‑151 (MQ=255)
 actacaagagatgcccgcatggcgcgaagagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaaggt                                                                                                                                                    >  2:691996‑M1/150‑151 (MQ=255)
    acaagagatgcccgcatggcgcgaagagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaaggtcgt                                                                                                                                                 <  1:772876‑M1/5‑1 (MQ=255)
     caagagatgcccgcatggcgcgaagagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaaggtcg                                                                                                                                                  >  2:635841‑M1/146‑149 (MQ=255)
      aagagatgcccgcatggcgcgaagagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaaggtcgtcg                                                                                                                                               <  2:739445‑M1/7‑1 (MQ=255)
       agagatgcccgcatggcgcgaagagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaaggtcgtcgt                                                                                                                                              <  2:278835‑M1/8‑1 (MQ=255)
          gatgcccgcatggcgcgaagagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCg                                                                                                                                           >  1:1163815‑M1/141‑151 (MQ=255)
             gcccgcatggcgcgaagagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGtt                                                                                                                                        >  1:797018‑M1/138‑151 (MQ=255)
               ccgcatggcgcgaagagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGa                                                                                                                                      >  1:384772‑M1/136‑151 (MQ=255)
                cgcatggcgcgaagagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGa                                                                                                                                      <  2:995325‑M1/16‑1 (MQ=255)
                 gcatggcgcgaagagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGaaa                                                                                                                                    >  1:765043‑M1/134‑151 (MQ=255)
                    tggcgcgaagagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAAcc                                                                                                                                  <  1:256238‑M1/20‑1 (MQ=255)
                          gaagagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGaa                                                                                                                            <  2:998918‑M1/26‑1 (MQ=255)
                            agagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACa                                                                                                                         >  2:1001456‑M1/123‑151 (MQ=255)
                             gagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAAc                                                                                                                          <  2:390691‑M1/28‑1 (MQ=255)
                             gagtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAg                                                                                                                        >  1:667868‑M1/122‑151 (MQ=255)
                              agtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGt                                                                                                                       >  1:964566‑M1/121‑151 (MQ=255)
                               gtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGt                                                                                                                       <  2:729055‑M1/31‑1 (MQ=255)
                               gtatgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGt                                                                                                                       <  2:370361‑M1/31‑1 (MQ=255)
                                 atgaagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAac                                                                                                                    <  1:848081‑M1/34‑1 (MQ=255)
                                    aagaagcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACt                                                                                                                 >  2:747107‑M1/115‑151 (MQ=255)
                                        agcgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCtt                                                                                                             <  1:494375‑M1/41‑1 (MQ=255)
                                          cgttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTg                                                                                                            <  2:638369‑M1/42‑1 (MQ=255)
                                            ttgcacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGt                                                                                                           <  1:188907‑M1/43‑1 (MQ=255)
                                               cacgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACa                                                                                                      >  2:989426‑M1/104‑151 (MQ=255)
                                                 cgacggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGc                                                                                                    <  2:801980‑M1/50‑1 (MQ=255)
                                                    cggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTg                                                                                                 <  1:1096032‑M1/53‑1 (MQ=255)
                                                     ggcgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTg                                                                                                 <  2:820184‑M1/53‑1 (MQ=255)
                                                       cgtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATc                                                                                              <  2:1152566‑M1/56‑1 (MQ=255)
                                                        gtagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATc                                                                                              <  2:126140‑M1/56‑1 (MQ=255)
                                                         tagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCAc                                                                                            <  2:543028‑M1/58‑1 (MQ=255)
                                                         tagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCAc                                                                                            <  1:301676‑M1/58‑1 (MQ=255)
                                                         tagagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCAc                                                                                            <  2:749441‑M1/58‑1 (MQ=255)
                                                          agagttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCAcc                                                                                           <  2:350790‑M1/59‑1 (MQ=255)
                                                             gttccgtttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGcc                                                                                        <  1:966879‑M1/62‑1 (MQ=255)
                                                                    ttcctgaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGt                                                                                  >  2:448316‑M1/83‑150 (MQ=255)
                                                                         gaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAAca                                                                             >  2:852703‑M1/78‑150 (MQ=255)
                                                                         gaataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAAc                                                                              >  2:608227‑M1/78‑149 (MQ=255)
                                                                          aataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAAcagc                                                                           >  2:321638‑M1/77‑151 (MQ=255)
                                                                          aataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAAcagc                                                                           <  1:480704‑M1/75‑1 (MQ=255)
                                                                           ataatccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAAcagc                                                                           <  2:661264‑M1/75‑1 (MQ=255)
                                                                           ataatccggaacgtgtcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGACGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAAcagc                                                                           <  2:82629‑M1/75‑1 (MQ=255)
                                                                                ccggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAATAGCAGGAt                                                                      <  2:288844‑M1/80‑1 (MQ=255)
                                                                                  ggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACt                                                                   <  1:631587‑M1/83‑1 (MQ=255)
                                                                                  ggaacgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACt                                                                   >  2:455461‑M1/69‑151 (MQ=255)
                                                                                      cgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGa                                                                 <  1:494442‑M1/85‑1 (MQ=255)
                                                                                      cgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGa                                                                 <  2:322106‑M1/85‑1 (MQ=255)
                                                                                      cgtttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGa                                                                 <  2:942987‑M1/85‑1 (MQ=255)
                                                                                         ttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGcc                                                             >  2:1129690‑M1/62‑150 (MQ=255)
                                                                                         ttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGcc                                                             >  1:1004181‑M1/62‑150 (MQ=255)
                                                                                         ttcgatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGcc                                                             >  1:216996‑M1/62‑150 (MQ=255)
                                                                                            gatgctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACCCTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGa                                                         >  1:1131270‑M1/59‑151 (MQ=255)
                                                                                               gctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATg                                                      <  2:151752‑M1/96‑1 (MQ=255)
                                                                                               gctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATg                                                      <  2:1154447‑M1/96‑1 (MQ=255)
                                                                                                ctgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATg                                                      <  1:125045‑M1/96‑1 (MQ=255)
                                                                                                 tgatggcaccttaaccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCg                                                    <  2:758265‑M1/98‑1 (MQ=255)
                                                                                                               ccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCt                                      <  1:959381‑M1/112‑1 (MQ=255)
                                                                                                               ccttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCt                                      >  1:350705‑M1/40‑151 (MQ=255)
                                                                                                                 ttgcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTgg                                    <  1:476858‑M1/114‑1 (MQ=255)
                                                                                                                   gcgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGAc                                  >  2:76723‑M1/36‑151 (MQ=255)
                                                                                                                    cgcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGAc                                  <  1:1065402‑M1/116‑1 (MQ=255)
                                                                                                                     gcgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACa                                 <  2:913496‑M1/117‑1 (MQ=255)
                                                                                                                      cgttatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACaa                                <  2:1144025‑M1/118‑1 (MQ=255)
                                                                                                                         tatgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACaaaaa                             >  2:687307‑M1/30‑150 (MQ=255)
                                                                                                                          atgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAGCAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAAc                            >  2:256294‑M1/29‑150 (MQ=255)
                                                                                                                           tgtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACgg                          <  2:869271‑M1/124‑1 (MQ=255)
                                                                                                                            gtcgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACgg                          >  1:278766‑M1/27‑150 (MQ=255)
                                                                                                                              cgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCt                        >  2:617512‑M1/25‑150 (MQ=255)
                                                                                                                              cgcttggcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTg                       >  2:494499‑M1/25‑151 (MQ=255)
                                                                                                                                    gcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGcc                    <  1:339239‑M1/130‑1 (MQ=255)
                                                                                                                                    gcgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGa                 >  2:311092‑M1/19‑151 (MQ=255)
                                                                                                                                     cgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTg                       >  2:126456‑M1/18‑144 (MQ=255)
                                                                                                                                     cgaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTg                       <  1:126423‑M1/127‑1 (MQ=255)
                                                                                                                                      gaaccggatgataaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTTCCGCTGGACAAAAACGGCTGGCcagtc                >  1:342545‑M1/17‑148 (MQ=255)
                                                                                                                                      gaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGAc                >  2:592220‑M1/17‑150 (MQ=255)
                                                                                                                                      gaaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACg               <  1:333051‑M1/135‑1 (MQ=255)
                                                                                                                                       aaccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGact               >  2:84374‑M1/16‑149 (MQ=255)
                                                                                                                                         ccggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCg            <  2:849763‑M1/138‑1 (MQ=255)
                                                                                                                                          cggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCg            <  1:802186‑M1/138‑1 (MQ=255)
                                                                                                                                           ggatgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGAc          <  1:692587‑M1/140‑1 (MQ=255)
                                                                                                                                             atgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCa        >  1:370272‑M1/10‑151 (MQ=255)
                                                                                                                                              tgagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCAt       <  1:8452‑M1/143‑1 (MQ=255)
                                                                                                                                               gagaaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATa      >  1:739279‑M1/8‑151 (MQ=255)
                                                                                                                                                  aaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAAt    >  1:638219‑M1/5‑150 (MQ=255)
                                                                                                                                                  aaagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATg   <  1:689330‑M1/147‑1 (MQ=255)
                                                                                                                                                   aagGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATgg  <  1:90802‑M1/148‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TACAAGAGATGCCCGCATGGCGCGAAGAGTATGAAGAAGCGTTGCACGACGGCGTAGAGTTCCGTTTCCTGAATAATCCGGAACGTTTCGATGCTGATGGCACCTTAACCTTGCGCGTTATGTCGCTTGGCGAACCGGATGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGG  >  NC_000913/3018200‑3018497

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