Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC NC_000913 4,001,644 Δ5 bp 100% coding (220‑224/951 nt) corA → magnesium/nickel/cobalt transporter

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4001644 4001648 5 60 [0] [0] 60 corA magnesium/nickel/cobalt transporter

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 40016430 (0.000)59 (0.960) 46/492 NT 100% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001649 = 0 (0.000)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

TGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCGAACAGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATC  >  NC_000913/4001358‑4001923
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                        
tGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATct                                                                                                                                                                                                                                                                                        >  2:514668‑M1/1‑286 (MQ=255)
       cAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttc                                                                                                                                                                                                                                                                                <  2:475671‑M1/289‑11 (MQ=255)
               ttCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttga                                                                                                                                                                                                                                                                        >  1:63588‑M1/1‑271 (MQ=255)
                            ccGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaaga                                                                                                                                                                                                                                                            >  2:361582‑M1/1‑258 (MQ=255)
                                     ccGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAGCCCCTGGTAAATGCAGGATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAGCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACAGCGAAGCATCGGCACGTTTCTTTGAAGACGAAGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggt                                                                                                                                                                                                                                                  <  1:108692‑M1/289‑41 (MQ=255)
                                     ccGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggt                                                                                                                                                                                                                                                  >  2:224980‑M1/1‑249 (MQ=255)
                                     ccGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggt                                                                                                                                                                                                                                                  <  1:227437‑M1/289‑41 (MQ=255)
                                       gATATTTTAAGCATTGGGGGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaa                                                                                                                                                                                                                                                <  2:247031‑M1/289‑43 (MQ=255)
                                              tAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccact                                                                                                                                                                                                                                         <  1:514668‑M1/289‑50 (MQ=255)
                                              tAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccact                                                                                                                                                                                                                                         <  2:276945‑M1/289‑50 (MQ=255)
                                                     tGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcat                                                                                                                                                                                                                                  >  2:264588‑M1/1‑233 (MQ=255)
                                                      gggAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATccccttcttcttctttgaagatgcggaaggtcaagccggtaactccactgtggcatt                                                                                                                                                                                                                                 >  2:509992‑M1/1‑232 (MQ=255)
                                                                 tCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtg                                                                                                                                                                                                                      >  1:3109‑M1/1‑221 (MQ=255)
                                                                       tGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggt                                                                                                                                                                                                                 <  1:476554‑M1/288‑74 (MQ=255)
                                                                           cgcATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtct                                                                                                                                                                                                            >  1:361115‑M1/1‑211 (MQ=255)
                                                                             cATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccact                                                                                                                                                                                                                                         <  1:554727‑M1/258‑50 (MQ=255)
                                                                             cATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccact                                                                                                                                                                                                                                         >  2:554727‑M1/1‑209 (MQ=255)
                                                                                       ggAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcg                                                                                                                                                                                                <  2:361115‑M1/289‑91 (MQ=255)
                                                                                        gAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcg                                                                                                                                                                                                <  2:3109‑M1/288‑91 (MQ=255)
                                                                                        gAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcg                                                                                                                                                                                                <  1:224980‑M1/288‑91 (MQ=255)
                                                                                         aaaaTAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtg                                                                                                                                                                                              <  2:355197‑M1/289‑93 (MQ=255)
                                                                                          aaaTAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtga                                                                                                                                                                                             >  1:495062‑M1/1‑196 (MQ=255)
                                                                                                      gACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcc                                                                                                                                                                                 <  2:427811‑M1/289‑106 (MQ=255)
                                                                                                           ggCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgctt                                                                                                                                                                            <  1:361582‑M1/289‑111 (MQ=255)
                                                                                                            gCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgctt                                                                                                                                                                            <  1:264588‑M1/288‑111 (MQ=255)
                                                                                                                                      ttGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgta                                                                                                                                                 >  1:459549‑M1/1‑152 (MQ=255)
                                                                                                                                            aTGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagt                                                                                                                                           <  2:495062‑M1/289‑144 (MQ=255)
                                                                                                                                             tGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgag                                                                                                                                                                                            >  2:327502‑M1/1‑145 (MQ=255)
                                                                                                                                             tGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgag                                                                                                                                                                                            <  1:327502‑M1/239‑95 (MQ=255)
                                                                                                                                                agttgggattggTCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTAGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcga                                                                                                                                        <  1:217736‑M1/283‑147 (MQ=255)
                                                                                                                                                  tATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagt                                                                                                                                           >  2:299093‑M1/1‑140 (MQ=255)
                                                                                                                                                  tATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagt                                                                                                                                           <  1:299093‑M1/283‑144 (MQ=255)
                                                                                                                                                            tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgatggtagacggta                                                                                                                            <  2:63588‑M1/288‑159 (MQ=255)
                                                                                                                                                            tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgctt                                                                                                                                                                            >  1:138643‑M1/1‑130 (MQ=255)
                                                                                                                                                            tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgctt                                                                                                                                                                            <  2:138643‑M1/240‑111 (MQ=255)
                                                                                                                                                                gTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggt                                                                                                                                                                                                                 >  1:323903‑M1/1‑126 (MQ=255)
                                                                                                                                                                gTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggt                                                                                                                                                                                                                 <  2:323903‑M1/199‑74 (MQ=255)
                                                                                                                                                                                           cgcgTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtag                                                                                                                                                >  2:469598‑M1/1‑99 (MQ=255)
                                                                                                                                                                                                   aTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttctccgttgaagatgcggaagatcacgcaggtaactccactgtggcatttaccatccgtgatggtagtctgtttactctgcgtgagcgtgaactgcccgctcttcgtctgtatcgtatgcgtgcccgatgccagtcgatggtagacggcaacgcagacgagttgatgccggtacagtccgcacccaatat                                                                                    >  1:81306‑M1/1‑91 (MQ=255)
                                                                                                                                                                                                      tGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgatggtagacggtaacgcctacgagttgctgctggatctgttcgaaaccaaaatcga                                                                                 <  2:138466‑M1/289‑202 (MQ=255)
                                                                                                                                                                                                        aaCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgatggtagacggtaacgcctacgagttgctgctggatctgttcgaaaccaaaatcgaac                                                                               >  2:530525‑M1/1‑86 (MQ=255)
                                                                                                                                                                                                              gCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtggccagtcgatggtagacggtaacgccgacgagttgctgctggatctgttcgacaccaaaatcgaacagctg                                                                          >  2:210342‑M1/1‑80 (MQ=255)
                                                                                                                                                                                                              gCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgt                                                                                                                                                                       >  1:555224‑M1/1‑80 (MQ=255)
                                                                                                                                                                                                              gCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgt                                                                                                                                                                       <  2:555224‑M1/195‑116 (MQ=255)
                                                                                                                                                                                                                        ggCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgatggtagacggt                                                                                                                             <  2:176903‑M1/227‑158 (MQ=255)
                                                                                                                                                                                                                        ggCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgatggtagacggt                                                                                                                             >  1:176903‑M1/1‑70 (MQ=255)
                                                                                                                                                                                                                            aCCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgatggtagacggtaacgcctacgagttgctgctggatctgttcgaaaccaaaatcgaacagttggcagatgaaattgaa                                                           <  1:530525‑M1/289‑224 (MQ=255)
                                                                                                                                                                                                                                       aCTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgatggtagacggtaacgcctacgagttgctgctggatctgttcgaaaccaaaatcgaacagttggcagatgaaattgaaaatatctata                                                 <  1:400003‑M1/288‑234 (MQ=255)
                                                                                                                                                                                                                                         tGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcctccttcttctttgaagatgcggaagatcacgcaggtaactccactgtggcatttaccatccgtgatggtcgtctgtctactctgcgtgagcgtgaactgcccgattttcgtctgtatcgtatgcgtgcccgcagccagccgatggtagacggtaacgcccacgacttgctcctggatctgatcgaacctaaaatccaacagttggatagcacaatggaaaaccaccatagc                                               >  1:290097‑M1/1‑53 (MQ=255)
                                                                                                                                                                                                                                                cATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgatggtagacggtaacgcctacgagttgctgctggatctgttcgaaaccaaaatcgaacagttggcagatgaaattgaaaatttctatacggacctgga                                       >  1:73146‑M1/1‑46 (MQ=255)
                                                                                                                                                                                                                                                     aaGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgatggtagacggtaacgcctacgagttgctgctggatctgttcgaaaccaaaatcgaacagttggcagatgaaattgaaaatatctatagcgacctggagcagt                                  >  1:224382‑M1/1‑41 (MQ=255)
                                                                                                                                                                                                                                                         aTCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgatggtagacggtaacgcctacgagttgctgctggatctgttcgaaaccaaaatcgaacagttggcagatgaaattgaaaatatctatagcgacctgg                                         <  1:559077‑M1/278‑242 (MQ=255)
                                                                                                                                                                                                                                                         aTCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgatggtagacggtaacgcctacgagttgctgctggatctgttcgaaaccaaaatcgaacagttggcagatgaaattgaaaatatctatagcgacctgg                                         >  2:559077‑M1/1‑37 (MQ=255)
                                                                                                                                                                                                                                                           cGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagt                                                                                                                                           >  1:203593‑M1/1‑35 (MQ=255)
                                                                                                                                                                                                                                                           cGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagt                                                                                                                                           <  2:203593‑M1/178‑144 (MQ=255)
                                                                                                                                                                                                                                                                   ttCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgatggtagacggtaacgcctacgagttgctgctggatctgttcgaaaccaaaatcgaacagttggcagatgaaattgaaaatatctatagcgacctggagcagttgagccgggtgatt                    >  1:140313‑M1/1‑27 (MQ=255)
                                                                                                                                                                                                                                                                   ttCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgatggtagacggtaacgcctacgagttgctgctggatctgttcgaaaccaaaatcgaacagttggcagatgaaattgaaaatatctatagcgacctggagcagttgagccgggtgatt                    >  2:146372‑M1/1‑27 (MQ=255)
                                                                                                                                                                                                                                                                         gAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgatggtagacggtaacgcctacgagttgctgctggatctgttcgaaaccaaaatcgaacagttggcagatgaaattgaaaatatctatagcgacctggagcagttgagccgggtgattatggaa              <  2:224382‑M1/289‑269 (MQ=255)
                                                                                                                                                                                                                                                                               gacgacgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgatggtagacggtaacgcctacgagttgctgctggatctgttcgaaaccaaaatcgaacagttggcagttgaaattgaaaatatctatagcgacctggagccgtttagccgggttatttgtgaagggca         >  1:406897‑M1/1‑15 (MQ=255)
                                                                                                                                                                                                                                                                                     ggCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagccagtcgatggtagacggtaacgcctacgagttgctgctggatctgttcgaaaccaaaatcgaacagttggcagatgaaattgaaaatatctatagcgacctggagcagttgagccgggtgattatggaagggcatcagggc  <  1:146372‑M1/289‑281 (MQ=255)
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                        
TGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCGAACAGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATC  >  NC_000913/4001358‑4001923

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