Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,976,679 Δ1,824 bp IS1‑mediated [motA]flhD [motA], flhC, flhD

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1976679 1979129–1978504 1826–2451 73 [0] [31] 32 [motA]–[insA‑5] [motA],flhC,flhD,insB‑5,[insA‑5]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 19766780 (0.000)73 (0.750) 51/422 0.8 100% coding (462/888 nt) motA motility protein A
?NC_000913 = 3584195 NA (NA)noncoding (768/768 nt) IS1 repeat region

TTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCA  >  NC_000913/1976441‑1976913
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
ttCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGgg                                                                                                                                                                                                                                           >  1:958399‑M1/1‑238 (MQ=255)
  cGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggta                                                                                                                                                                                                                                         <  1:923207‑M1/240‑5 (MQ=255)
      gcTTTTCTGACGTAAAACAGTCGCTAATGGAGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaat                                                                                                                                                                                                                                       <  1:192135‑M1/238‑7 (MQ=255)
        ttttCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgact                                                                                                                                                                                                                                   <  2:1049424‑M1/240‑11 (MQ=255)
        ttttCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgact                                                                                                                                                                                                                                   <  2:212608‑M1/240‑11 (MQ=255)
          ttCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactcc                                                                                                                                                                                                                                 <  1:184264‑M1/240‑13 (MQ=255)
                    aaaCAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCAAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGAACGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACTCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattg                                                                                                                                                                                                                        >  2:476130‑M1/1‑218 (MQ=255)
                     aaCAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgat                                                                                                                                                                                                                      <  2:86019‑M1/240‑24 (MQ=255)
                        aGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgact                                                                                                                                                                                                                                   <  2:705343‑M1/224‑11 (MQ=255)
                                                gTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccga                                                                                                                                                                                           <  2:1019767‑M1/240‑51 (MQ=255)
                                                       aaTAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgtttt                                                                                                                                                                                                              <  1:186538‑M1/214‑32 (MQ=255)
                                                        aTAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttg                                                                                                                                                                                   <  2:1057869‑M1/240‑59 (MQ=255)
                                                         tAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgt                                                                                                                                                                                  >  1:848376‑M1/1‑181 (MQ=255)
                                                                                    cGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacg                                                                                                                                                       >  1:166036‑M1/1‑154 (MQ=255)
                                                                                      cATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCCAACGCCGGAAGTGAGTCCCCTACCAGCGCCATACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggta                                                                                                                                                                                                                                         >  2:16421‑M1/1‑152 (MQ=255)
                                                                                          tgCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcg                                                                                                                                                 >  2:501527‑M1/1‑148 (MQ=255)
                                                                                           gCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcg                                                                                                                                                 >  2:767800‑M1/1‑147 (MQ=255)
                                                                                                     cACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcc                                                                                                                                       <  1:1130243‑M1/239‑103 (MQ=255)
                                                                                                        ccAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagc                                                                                                                                   <  1:953528‑M1/240‑107 (MQ=255)
                                                                                                            cTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgt                                                                                                                                >  2:413487‑M1/1‑130 (MQ=255)
                                                                                                            cTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagcc                                                                                                                                  >  2:744531‑M1/1‑130 (MQ=255)
                                                                                                                      ggACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgct                                                                                                                      >  2:74934‑M1/1‑120 (MQ=255)
                                                                                                                      ggACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcag                                                                                                                                                                                                      >  1:898779‑M1/1‑120 (MQ=255)
                                                                                                                          gATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcct                                                                                                                  <  1:378761‑M1/239‑124 (MQ=255)
                                                                                                                                         aCGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcct                                                                                                                  <  2:931187‑M1/224‑124 (MQ=255)
                                                                                                                                          cgcgTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgcc                                                                                                 >  1:326221‑M1/1‑100 (MQ=255)
                                                                                                                                               gAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgag                                                                                                                                                           >  2:819161‑M1/1‑95 (MQ=255)
                                                                                                                                               gAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgag                                                                                                                                                           <  1:819163‑M1/177‑83 (MQ=255)
                                                                                                                                                aaCGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgt                                                                                                                                >  2:923180‑M1/1‑94 (MQ=255)
                                                                                                                                                    aCCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattc                                                                                        <  1:507863‑M1/239‑150 (MQ=255)
                                                                                                                                                     cccccATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctcca                                                                                                                                                                      <  1:853072‑M1/160‑72 (MQ=255)
                                                                                                                                                      ccccATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgct                                                                                     <  1:656040‑M1/240‑153 (MQ=255)
                                                                                                                                                           ttACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgt                                                                                 <  1:684913‑M1/239‑157 (MQ=255)
                                                                                                                                                               aGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttttgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagat                                                                                                             >  1:225709‑M1/1‑79 (MQ=255)
                                                                                                                                                               aGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctca                                                                                            >  1:896712‑M1/1‑79 (MQ=255)
                                                                                                                                                                 ccGCAACAATACCAAACGCCGGAAGTTAGTCCCCGACCAGCGCCCGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgc                                                                          >  2:64455‑M1/1‑77 (MQ=255)
                                                                                                                                                                 ccGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttg                                                                                                                                                                                   <  1:695697‑M1/135‑59 (MQ=255)
                                                                                                                                                                 ccGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttg                                                                                                                                                                                   >  2:695695‑M1/1‑77 (MQ=255)
                                                                                                                                                                   gCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgctt                                                                        >  2:966124‑M1/1‑75 (MQ=255)
                                                                                                                                                                    caacaaTACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgctt                                                                        >  1:917527‑M1/1‑74 (MQ=255)
                                                                                                                                                                           aCCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgatt                                                                 <  2:229250‑M1/239‑173 (MQ=255)
                                                                                                                                                                               aaCGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgt                                                             >  1:858051‑M1/1‑63 (MQ=255)
                                                                                                                                                                                    cGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagct                                                       >  1:938217‑M1/1‑58 (MQ=255)
                                                                                                                                                                                         gTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgca                                                          <  2:105464‑M1/232‑180 (MQ=255)
                                                                                                                                                                                            aGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttccctt                                                >  1:960470‑M1/1‑50 (MQ=255)
                                                                                                                                                                                              tCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcag                                             >  2:534986‑M1/1‑48 (MQ=255)
                                                                                                                                                                                                cccGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggc                                           <  2:63138‑M1/240‑195 (MQ=255)
                                                                                                                                                                                                cccGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggc                                           <  1:1134158‑M1/240‑195 (MQ=255)
                                                                                                                                                                                                cccGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggc                                           <  1:311098‑M1/240‑195 (MQ=255)
                                                                                                                                                                                                cccGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggc                                           <  1:371799‑M1/240‑195 (MQ=255)
                                                                                                                                                                                                cccGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggc                                           <  1:1064625‑M1/240‑195 (MQ=255)
                                                                                                                                                                                                  cGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgg                                         <  2:166035‑M1/240‑197 (MQ=255)
                                                                                                                                                                                                   gACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgg                                         <  2:811120‑M1/239‑197 (MQ=255)
                                                                                                                                                                                                        gcgcCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattc                                    >  1:442449‑M1/1‑38 (MQ=255)
                                                                                                                                                                                                        gcgcCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccag                                                                                                                                    <  1:122025‑M1/143‑106 (MQ=255)
                                                                                                                                                                                                        gcgcCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccag                                                                                                                                    >  2:122025‑M1/1‑38 (MQ=255)
                                                                                                                                                                                                        gcgcCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgatttcgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattc                                    >  2:322390‑M1/1‑38 (MQ=255)
                                                                                                                                                                                                         cgcCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcat                                  >  1:652643‑M1/1‑37 (MQ=255)
                                                                                                                                                                                                               aCTGTTCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgtgttatgttcagataatgcccgatgactttgtcatgcagctccgccgagtttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcatcagattcaggttatgccgctcaattcgct                                                                                     >  1:96224‑M1/1‑31 (MQ=255)
                                                                                                                                                                                                                   ttCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagnttcaggttatgccgctcaattcgctgcgtatat                                                                             >  2:18314‑M1/1‑27 (MQ=255)
                                                                                                                                                                                                                   ttCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatacagcggcc                         >  2:864801‑M1/1‑27 (MQ=255)
                                                                                                                                                                                                                   ttCGCCGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatat                                                                             >  1:425337‑M1/1‑27 (MQ=255)
                                                                                                                                                                                                                       ccGGGACTTCTGCCTCGCTTTCGggtaatgactccaacttattggtagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagcttccccttcaggcgg                                         <  2:407994‑M1/219‑197 (MQ=255)
                                                                                                                                                                                                                          ggACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattc                                                                                        <  1:1118045‑M1/169‑150 (MQ=255)
                                                                                                                                                                                                                           gACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatacagcggccagccatccg                <  2:640228‑M1/240‑222 (MQ=255)
                                                                                                                                                                                                                           gACTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatacagcggccagccatccg                <  1:615441‑M1/240‑222 (MQ=255)
                                                                                                                                                                                                                             cTTCTGCCTCGCTTTCGggtaatgactccaacttattgatagggttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcagac                                >  2:595834‑M1/1‑17 (MQ=255)
                                                                                                                                                                                                                                cTGCCTCGCTTTCGtgtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatacagcggccagccatccgtca             <  2:1030162‑M1/238‑225 (MQ=255)
                                                                                                                                                                                                                                cTGCCTCGCTTTCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatacagcggccagccatccgtcatc           >  1:154843‑M1/1‑14 (MQ=255)
                                                                                                                                                                                                                                     tCGCTTTCGggtaatgactccaacttattgata                                                                                                                                                                                                                     <  2:577743‑M1/33‑25 (MQ=255)
                                                                                                                                                                                                                                         tttCGggtaatgactccaacttattgttagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatacagcggccagccatccgtcatccatatcacc  <  2:363720‑M1/240‑236 (MQ=255)
                                                                                                                                                                                                                                         tttCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatacagcggccagccatccgtcatccatatcacc  <  2:183784‑M1/240‑236 (MQ=255)
                                                                                                                                                                                                                                         tttCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatacagcggccagccatccgtcatccatatcacc  <  1:501529‑M1/240‑236 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
TTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTTCGATTTCGAAGGTGTTCATGTGACCGCTGATAATCAGGCGCAGATAATCGACGATAAAATCAAGCATGACGCTATCCGCGAGGATGCGTGGGTAGCTGGCGAAGATCTCGCTCTCACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCA  >  NC_000913/1976441‑1976913

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: