Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 364,237 Δ1,164 bp coding (906‑2069/3075 nt) lacZ ← beta‑galactosidase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 364237 365400 1164 40 [0] [0] 40 lacZ beta‑galactosidase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3642360 (0.000)40 (0.930) 37/468 0.3 100% coding (2070/3075 nt) lacZ beta‑galactosidase
?NC_000913 365401 = 0 (0.000)coding (905/3075 nt) lacZ beta‑galactosidase

TCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCGGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGTAGGTAGTCACGCAACTCGCCGCACATCTGAACTTCAGCCTCCAGTACAGCGCGGCTGAAATCATCATTAAAGCGAGTGGCAACATGGAA  >  NC_000913/365165‑365630
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                      
gctgtgaaagaaagcctgtcgggcggttaaattgccaacgcttattacccagctcgatgcaaaaatccatttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGttt                                                                                                                                                                                                                          <  1:335215‑M2/15‑1 (MQ=255)
   gtgaaagaaagcctgactggcggttaaattgccaacgcttattacccagctcgatgcaaaaatccatttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTTCg                                                                                                                                                                                                                       <  2:89596‑M2/18‑1 (MQ=255)
      aaagaaagcctgacgggcggttaaattgccaacgcttattacccagctcgatgcaaaaatccatttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTTCGACg                                                                                                                                                                                                                    <  2:228420‑M2/21‑1 (MQ=255)
           aagcctgactggcggttaaattgccaacgcttattacccagctcgatgcaaaaagccatttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTTCGACGTTCAg                                                                                                                                                                                                               <  2:107690‑M2/26‑1 (MQ=255)
             gcctgactggcggttaaattgccaacgcttattacccagctcgatgcaaaaatccatttcgctggtggtcagatgcgggatggcgtggtacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTACGACGTTCAGAc                                                                                                                                                                                                             <  2:437900‑M2/28‑1 (MQ=255)
                  actttcggttaaattgccaactcttattacccagctcgatgcaaaaatccatttcgctggtggtcagatgcgggatggcgtggtacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtacggtgagccagggttgcccggcgctcCACAGTTTCGGGTTTTCGACGTTCAGACGTAgt                                                                                                                                                                                                        <  2:290523‑M2/33‑1 (MQ=255)
                     ggcggttaaattgccaacgcttattacccagctcgatgcaaaaatccatttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccaccgctgcccggcgctgatgtgcccggcttctggccatgcgggcgcgttcggttgcactacgcgtactgtgagccagaggtgccccgcgctcCACCGTTTCGGGGTTTTggcgttcagaaggcggttg                                                                                                                                                                                                     >  2:314166‑M2/216‑242 (MQ=255)
                      gcggttaaattgccaacgcttattacccagctcgatgcaaaaatccatttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgcccggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACCGTTTCGGGTTTTCGACGTTCAGCCGTAGtgtg                                                                                                                                                                                                     >  1:100696‑M2/215‑250 (MQ=255)
                      gcggttaaattgccaacgcttattacccagctcgatgcaaaaatccatttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcgggcgcgttcggttgcattacgcgtactgtgagccagagttgcccggggctcCTCAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGa                                                                                                                                                                                                   <  1:61659‑M2/37‑1 (MQ=255)
                             aattgccaacgcttattacccagctcgatgcaaaaatccatttcgctggtggtcagatgcgggatggcgtgggacgcggcggggcgcgtcacactgaggttttccgccaggcgcccctgctgcccggcgctggtgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtcctgtgagcccgaggtgcccggcgctcCCCAGTTTCGGGGTTTCGACGTTCAGACGTAGTGTTACGCGATc                                                                                                                                                                                            >  2:453073‑M2/208‑251 (MQ=255)
                                 gccaacgcttattacccagctcgatgcaaaaatccatttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGc                                                                                                                                                                                         >  1:341837‑M2/204‑250 (MQ=255)
                                 gccaacgcttattacccagctcgatgcaaaaatccatttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGGGACGCGATCGGc                                                                                                                                                                                         >  2:145222‑M2/204‑250 (MQ=255)
                                      cgcttattaccccgctcgatgcaaaaatccatttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACCGCTTCGGGTTTtcgacgttcagccgtagtgtgcgatcggcagagcgtcg                                                                                                                                                                                    >  2:233173‑M2/199‑234 (MQ=255)
                                      cgcttattacccagctcgatgcaaaaatccatttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggctctcCACAGTTTCGGGTTTTCGACGTTCAGACGTAGtgtg                                                                                                                                                                                                     <  1:233173‑M2/36‑1 (MQ=255)
                                          gattgccaagcgcgatgcaaaaatccatttcgctggtggtcagatgggggttggcgtgggacggggcggggatcgtcacactgaggttttccgccagacgccacttctgccagtcgcttatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgcgcCACAGTTTCGTGTTTGCGACGTTCAGACGTAGTGTGACGCGATCGGCAtaaccatc                                                                                                                                                                                <  2:38348‑M2/56‑3 (MQ=255)
                                           attacccagcgcgatgcaaaaatccatgtcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTGCGACGGTCAGACGTAGTGTGACGCGATCGGCATAAccacca                                                                                                                                                                               <  1:63365‑M2/57‑1 (MQ=255)
                                             tacccagctcgatgcaaaaatcaatttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgtttggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTTCGACGTTCAGACGTAGGGGGAAGGCAATGGCATAACCACCACGc                                                                                                                                                                            >  1:106961‑M2/192‑251 (MQ=255)
                                                ccagctcgatgcaaaaatccatttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTc                                                                                                                                                                          >  1:12960‑M2/189‑250 (MQ=255)
                                                ccagctcgatgcaaaaatccatttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGGTTTCGACGTTCCGACGGAGTGGGACGCGGTcggcataaccacccccctcc                                                                                                                                                                         >  2:361581‑M2/189‑244 (MQ=255)
                                                         tgcaaaaatccatttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgagggtttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGGTTTTGCCGTTCAGACGTAGTGGGACgcgatcggcataaccaccacccacatccaaaatt                                                                                                                                                                >  1:252376‑M2/180‑237 (MQ=255)
                                                                    atttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagcgttgcccggcgctcCACAGTTTCGGGTTCTCGACCTTCCGACGTAGTGTGACGCGATCGGCCTAACCACCACGCTCATCGATAATTTCACCCCCGa                                                                                                                                                      >  2:7381‑M2/169‑250 (MQ=255)
                                                                     tttcgctggtggtcagatgcgggatggcgtgagacgcggcggggagcgtcacactgcgggtttccgccagacgccactgctaccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACCGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGGTCGGCATAACCACCACGCTCATCGATAATTTcaccgccgca                                                                                                                                                     >  1:223991‑M2/168‑248 (MQ=255)
                                                                      ttcgctggtggtcagatgcgggatggcgtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGaaa                                                                                                                                                    <  1:361581‑M2/84‑1 (MQ=255)
                                                                                              ggcgtgggacgcggcggggagcgtcacactggggttttccgcccgacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcggactgtgagccagagttgcccggcgctcCCCAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGACTAACCACCACGCCCCTCGATAATTTCAACGccgaaagacgcggtgcaccctgcgcccagcg                                                                                                                           >  2:159959‑M2/143‑236 (MQ=255)
                                                                                                 gtgggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgggcccggcttctgaccatgcggtcgcgttcgggtgcactacgcggactgtgagccagagttgccccgcgctcCACAGTTTCGGGTTTTTGACGTTCAGACGTAGTGGGGAGCGAACGGCATAACCACCACGCCCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCAGGCGACCtgcggtt                                                                                                                        >  2:10493‑M2/140‑248 (MQ=255)
                                                                                                    ggacgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGGATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCAc                                                                                                                     <  1:162027‑M2/115‑1 (MQ=255)
                                                                                                       cgcggcggggagcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCt                                                                                                                  <  1:354439‑M2/118‑1 (MQ=255)
                                                                                                                  gcgtcacactgaggttttccgccagacgccactgctgccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCCGGCGACCTGCGTTTCACCCTGCCATAAAGa                                                                                                        >  2:420552‑M2/123‑250 (MQ=255)
                                                                                                                                         ggacgccactgctgcccggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggtggcactacgcgtcctgtgagccagagtttcccggcgctcCACAGTTTCGGGTTTTCGACTTTCAGACGTGGTGTGACGCGGTCGGCATAACCACCACGCTCATCGATACTTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGTAGGTAGTCACg                                                                                <  1:59327‑M2/152‑1 (MQ=255)
                                                                                                                                               cactgttgccaggggctgatttgcccggtttctgcccaggcggtcgcgttcggttgcactccgcgtactgtgagccagagttgcccggcgctcCCCAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGAGCGGCATAACCCCCCCGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGTAGGTAGTCACGCAACTc                                                                          <  1:307355‑M2/158‑1 (MQ=255)
                                                                                                                                                       ccaggcgctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGTAGGTAGTCACGCAACTCGCCGcac                                                                   >  1:107289‑M2/86‑250 (MQ=255)
                                                                                                                                                           gcgttgatgtgcccgtcttctgaccgtgcggtcgcgttcggttgccctccgcgtactgtgagccagggttccccggcgcccCACAGTTTCGGGTTTCCGCCGTTCAGACGTAGTGTGACGCGAGCGGCATACCCCCCACGCTCAGCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGTAGGTAGTCACGCAACTCGCCGCACATCTg                                                              <  2:274658‑M2/170‑1 (MQ=255)
                                                                                                                                                              ctgatgtgcccggcttctgaccatgcggtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGAAGGTAGTCACGCAACTCGCCGCACATCTGAAc                                                           >  1:405412‑M2/79‑251 (MQ=255)
                                                                                                                                                                   gtgcccggcttctgaccatgcggtcgcgttgggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGTAGGTAGTCACGCAACTCGCCGCACATCTGAACTTCa                                                       <  1:314166‑M2/177‑1 (MQ=255)
                                                                                                                                                                          gtttctgaccctggggtggggttcggtcgcactacgcgtactgtgagccagggttgcccggcgctcCACAGTTGCGGGTTTTCGACGTTCAGACGTAGTTTGACGCGATCGGAATACCCACCCCGCACACCGATAATTTAACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGTAGGTAGTCACGCAACTCGCCGCACATCTGAACTTCAGCCTCCAg                                               <  1:219809‑M2/185‑1 (MQ=255)
                                                                                                                                                                                         gtcgcgttcggttgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGTAGGTAGTCACGCAACTCGCCGCACATCTGAACTTCAGCCTCCAGTACAGCGCGGCTGa                                 >  2:309805‑M2/52‑250 (MQ=255)
                                                                                                                                                                                              gttcggctgcactacgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATAGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGTAGGTAGTCACGCAACTCGCCGCACATCTGAACTTCAGCCTCCAGTACAGAGCGGCTGAAatcat                           >  1:377270‑M2/47‑251 (MQ=255)
                                                                                                                                                                                                           ccgcgtactgtgagccagagttgcccggcgctcCACAGTTTCGGGTTTGCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGTAGGTAGTCACGCAACTCGCCGCACATCTGAACTTCAGCCTCCAGTACAGCGCGGCTGAAATCATCATTAAAGCGAGt              <  1:145222‑M2/218‑1 (MQ=255)
                                                                                                                                                                                                                       agccagagttgcccggcgctcCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGTAGGTAGTCACGCAACTCGCCGCACATCTGAACTTCAGCCTCCAGTACAGCGCGGCTGAAATCATCATTAAAGCGAGTGGCAACATGGa   >  1:314916‑M2/22‑250 (MQ=255)
                                                                                                                                                                                                                        gccagagttgcccggcgctcCACAGTTTCGGGTTACCGCCGTTCAGTCGTAGTGTCACGCGATCGGCCTAACCACCACGCTACTCGATAATTTCACCGCCGCCAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGTAGGTAGTCACGCCACTCCCCGCACATCTGCACTTCAGCCTCCAGTACAGACCCGCTGAAATCCTCATTAAAGCGAGCGGCAACATGGaa  >  1:422605‑M2/21‑250 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                      
TCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCGGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGTAGGTAGTCACGCAACTCGCCGCACATCTGAACTTCAGCCTCCAGTACAGCGCGGCTGAAATCATCATTAAAGCGAGTGGCAACATGGAA  >  NC_000913/365165‑365630

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