Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,066,152 Δ12,778 bp insH6[yoeF] insH6, yoeA, insD3, insC3, yoeA, isrC, flu, yeeR, yeeS, yeeT, cbeA, cbtA, yeeW, [yoeF]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2066152 2078929 12778 39 [0] [0] 39 insH6–[yoeF] insH6,yoeA,insD3,insC3,yoeA,isrC,flu,yeeR,yeeS,yeeT,cbeA,cbtA,yeeW,[yoeF]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 20661510 (0.000)39 (0.510) 33/492 0.9 100% intergenic (‑387/+154) cobU/insH6 cobinamide‑P guanylyltransferase/cobinamide kinase/IS5 family transposase and trans‑activator
?NC_000913 2078930 = 0 (0.000)coding (356/357 nt) yoeF putative uncharacterized protein YoeF

TTTCTGGCAACAAGCAGGAATCTCAGCGCTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCTTTGGTGAAAAGTTGGTTACGGGAAGCTGGAGCGCACAACGTGTTAGTAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGC  >  NC_000913/2078682‑2079166
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            
aataagatggcattacggaatagccagactttctcgtgttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAAt                                                                                                                                                                                                                                            <  1:71684‑M2/3‑1 (MQ=255)
  taagatggcattacggaatagccagactttctcgtgttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAAttt                                                                                                                                                                                                                                          >  1:234175‑M2/247‑251 (MQ=255)
      atggcattacggaatagccagactttctcgtgttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGAt                                                                                                                                                                                                                                      <  1:518057‑M2/9‑1 (MQ=255)
                        cagactttctcgtgttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagatcttgtaaaacaacgtaatagaattatccggaatggttgcgactatAATTTCGATTGGTTCTCGgatgaggat                                                                                                                                                                                                                    <  1:544577‑M2/27‑4 (MQ=255)
                          gactttctcgtgttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCggatgaggatga                                                                                                                                                                                                                  <  1:168304‑M2/29‑1 (MQ=255)
                                     gttaacgtccgacaaaatgcgtcattggcagcggtttctgtggagtgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGAt                                                                                                                                                                                                        >  2:613110‑M2/212‑250 (MQ=255)
                                     gttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGGTTGGTTCTCGGATGAGGATGATGACGGTGAt                                                                                                                                                                                                        >  2:329525‑M2/212‑250 (MQ=255)
                                           gtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATacac                                                                                                                                                                                                    >  1:94492‑M2/206‑248 (MQ=255)
                                                      tgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaat                                                                                                                                                                                      >  2:104994‑M2/195‑251 (MQ=255)
                                                                         tctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgcaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgcaatagaattatccggaatggaggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCCATGAGCAATGAGGGTTGCCgg                                                                                                                                                                   >  2:427853‑M2/176‑251 (MQ=255)
                                                                         tctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCg                                                                                                                                                                    <  2:461380‑M2/75‑1 (MQ=255)
                                                                            gtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCa                                                                                                                                                                <  2:123354‑M2/79‑1 (MQ=255)
                                                                            gtggattgacacattgcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCa                                                                                                                                                                <  2:375484‑M2/79‑1 (MQ=255)
                                                                                   gacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacgtgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTc                                                                                                                                                          <  2:343095‑M2/85‑1 (MQ=255)
                                                                                             attcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATaaa                                                                                                                                                >  2:444803‑M2/156‑250 (MQ=255)
                                                                                                cagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACgg                                                                                                                                            <  1:104994‑M2/99‑1 (MQ=255)
                                                                                                             aattcgactgcccgaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatcttcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTActtcgc                                                                                                                                <  2:320028‑M2/111‑1 (MQ=255)
                                                                                                                gctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTgcttcgcttc                                                                                                                             <  2:504328‑M2/114‑1 (MQ=255)
                                                                                                                  tactgccggaacaaccagctcaagggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatatgcgctgtaaaagaacttgtaacacaacgtaatagaattatccggactggtggagcctagAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGATTGCCGGGCAACCCTCCTTGAATAAAACGGGAATTACTTCGCTTCGcc                                                                                                                          >  2:469128‑M2/135‑251 (MQ=255)
                                                                                                                  tactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGcc                                                                                                                          >  1:518658‑M2/135‑251 (MQ=255)
                                                                                                                                       acgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGc                                                                                                     <  1:505198‑M2/138‑1 (MQ=255)
                                                                                                                                        cgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCt                                                                                                    >  1:555803‑M2/113‑251 (MQ=255)
                                                                                                                                        cgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCt                                                                                                    >  1:46907‑M2/113‑251 (MQ=255)
                                                                                                                                         gggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttttaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCtt                                                                                                   <  2:194384‑M2/140‑1 (MQ=255)
                                                                                                                                          ggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCttt                                                                                                  <  2:423907‑M2/141‑1 (MQ=255)
                                                                                                                                                                      gattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACTGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCt                                                                      <  2:245790‑M2/169‑1 (MQ=255)
                                                                                                                                                                       attcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaccgtaatcgaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGATTTTTGGAGATATCGCGGCGTTCtt                                                                     >  2:594909‑M2/82‑251 (MQ=255)
                                                                                                                                                                          caacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTgg                                                                  >  1:401196‑M2/79‑251 (MQ=255)
                                                                                                                                                                                       agatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCAtt                                                     >  2:410265‑M2/66‑251 (MQ=255)
                                                                                                                                                                                                 tgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGTGATATCGCGGCGTTCTTTGTAGAGCTCGGCATTTTTGttgatg                                           >  2:77245‑M2/56‑251 (MQ=255)
                                                                                                                                                                                                  gtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGt                                          >  2:378831‑M2/55‑251 (MQ=255)
                                                                                                                                                                                                  gtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTGgatgat                                            <  2:565957‑M2/195‑1 (MQ=255)
                                                                                                                                                                                                         aacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCAGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTc                                   <  2:555803‑M2/204‑1 (MQ=255)
                                                                                                                                                                                                                          aatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTggagc                                                                                    <  2:26153‑M2/155‑2 (MQ=255)
                                                                                                                                                                                                                          aatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGgaga                                                                                    >  1:26153‑M2/31‑185 (MQ=255)
                                                                                                                                                                                                                             agaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCAc               <  1:595538‑M2/224‑1 (MQ=255)
                                                                                                                                                                                                                                   atccggaatggtggcgactatAATTTCGGTTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCa         <  1:329525‑M2/230‑1 (MQ=255)
                                                                                                                                                                                                                                         aatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTgg   >  1:298167‑M2/16‑251 (MQ=255)
                                                                                                                                                                                                                                          atggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTACAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAGAGCGATCTTCATAGTCACCTTTCATGCTGGc  >  2:355524‑M2/15‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            
TTTCTGGCAACAAGCAGGAATCTCAGCGCTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCTTTGGTGAAAAGTTGGTTACGGGAAGCTGGAGCGCACAACGTGTTAGTAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGC  >  NC_000913/2078682‑2079166

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