Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC NC_000913 3,909,311 Δ774 bp 100% [pstA][pstC] [pstA], [pstC]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3909311 3910084 774 34 [0] [0] 34 [pstA]–[pstC] [pstA], [pstC]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 39093100 (0.000)30 (0.910) 28/484 NT 100% coding (129/891 nt) pstA phosphate ABC transporter permease
?NC_000913 3910085 = 0 (0.000)coding (313/960 nt) pstC phosphate ABC transporter permease

ACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCG  >  NC_000913/3909024‑3909589
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                       
acaGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgc                                                                                                                                                                                                                                                                                       >  1:27100‑M1/1‑287 (MQ=255)
   gACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgca                                                                                                                                                                                                                                                                                    >  2:266206‑M1/1‑284 (MQ=255)
     caaacaccaatggccgGCGCAGAGAGCGGAAGGTCGTTAATGAAGCGAATCACTCCTGCCAGCCAGGATTTCCGACCATATTCCGCCAGATAAACCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaag                                                                                                                                                                                                                                                                                  <  2:61463‑M1/282‑8 (MQ=255)
              aTCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtaag                                                                                                                                                                                                                                                                         >  1:317866‑M1/1‑273 (MQ=255)
                       caagaaggcagaaggTCGTAAGTGAAGCGAATCCCTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacaggg                                                                                                                                                                                                                                                                 <  2:317866‑M1/281‑25 (MQ=255)
                                     gTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactc                                                                                                                                                                                                                                                   >  2:304820‑M1/1‑250 (MQ=255)
                                                                                          ccAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatc                                                                                                                                                                                             >  1:165986‑M1/1‑197 (MQ=255)
                                                                                                     cccGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgc                                                                                                                                                                                  <  1:266206‑M1/289‑104 (MQ=255)
                                                                                                               tGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcag                                                                                                                                                                                                                               >  2:93433‑M1/1‑176 (MQ=255)
                                                                                                               tGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcag                                                                                                                                                                                                                               <  1:93433‑M1/234‑59 (MQ=255)
                                                                                                                      cggcgTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcg                                                                                                                                                                 >  2:43365‑M1/1‑169 (MQ=255)
                                                                                                                          gcaccAAAAACCGTGGCCCACAAAATTAACAGCCCCCGCCCCGCCAAAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGTCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACAGTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtg                                                                                                                                                              <  2:165986‑M1/286‑124 (MQ=255)
                                                                                                                          gTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtg                                                                                                                                                              <  1:275053‑M1/288‑124 (MQ=255)
                                                                                                                            aCCGAATACCGTCGCCCACAAAATAACCAGCCCCCTACACGCCAGGGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgca                                                                                                                                                            <  2:27100‑M1/288‑126 (MQ=255)
                                                                                                                                       tGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggt                                                                                                                                                               >  2:21935‑M1/1‑152 (MQ=255)
                                                                                                                                       tGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggt                                                                                                                                                               <  1:21935‑M1/274‑123 (MQ=255)
                                                                                                                                       tGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagacatcgttcggtgcatcccactctttg                                                                                                                                                >  2:318121‑M1/1‑152 (MQ=255)
                                                                                                                                        ggCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGCCATTAAAATCCAAATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtccccaaagtaccgtagctcggcaacagcgcccagtagatatcgttcggtgcatcccactctttgg                                                                                                                                               >  2:213695‑M1/1‑151 (MQ=255)
                                                                                                                                          cccACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGACCAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatgagcagcgcgataaactaagtcaccaacgcaccgtagatcgccaccagcgccacgtagatatcgatcggtgcatccacctctttggt                                                                                                                                              >  1:298530‑M1/1‑149 (MQ=255)
                                                                                                                                                         gCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaccgaagt                                                                                                                                                                                                                >  1:276534‑M1/1‑134 (MQ=255)
                                                                                                                                                         gCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagt                                                                                                                                                                                                                <  2:276534‑M1/207‑74 (MQ=255)
                                                                                                                                                            cGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggaaccagcgccccgtagatatcgttcggttcatcccactctttggtccataggaaagcacgacca                                                                                                                           >  2:89616‑M1/1‑131 (MQ=255)
                                                                                                                                                                         agCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtccataggaaagccagaccaaatttctgaatg                                                                                                               >  2:182459‑M1/1‑118 (MQ=255)
                                                                                                                                                                                    taccaccacctTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtccataggaaagccagaccaaatttctgaatgctcggccaggag                                                                                                   <  1:304820‑M1/288‑183 (MQ=255)
                                                                                                                                                                                      ccaccaccTGCCGTATTGGGCGGCGGCGTCATGTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtccataggaaagccagaccaaatttctgaatgctcggccaggagga                                                                                                 <  1:89616‑M1/289‑185 (MQ=255)
                                                                                                                                                                                                        ggcggcggcGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtccataggaaagccagaccaaatttctgaatgctcggccaggaggagatgatcagagagacaat                                                                               >  2:83483‑M1/1‑87 (MQ=255)
                                                                                                                                                                                                             cggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtccataggaaagccagaccaaatttctgaatgctcggccaggaggagatgatcag                                                                                        <  1:71859‑M1/275‑194 (MQ=255)
                                                                                                                                                                                                             cggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtccataggaaagccagaccaaatttctgaatgctcggccaggaggagatgatcag                                                                                        >  2:71859‑M1/1‑82 (MQ=255)
                                                                                                                                                                                                                                       gCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtccataggaaagccagaccaaatttctgaatgctcggccaggaggagatgatcagagagacaataatgccacccaacatcaatagcacaatcagc                                                >  1:6826‑M1/1‑56 (MQ=255)
                                                                                                                                                                                                                                              aCCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtccataggaaagccagaccaaatttctgaatgctcggccaggaggagatgatcagagagacaataatgccacccaacatcaatagcacaatcagcgccgcca                                         <  1:43365‑M1/289‑241 (MQ=255)
                                                                                                                                                                                                                                                                gATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtccataggaaagccagaccaaatttctgaatgctcggccaggaggagatgatcagagagacaataatgccacccaacatcaatagcacaatcagcgccgccagttttaccagcacgct                         <  1:182459‑M1/287‑257 (MQ=255)
                                                                                                                                                                                                                                                                                 tCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtccataggaaagccagaccaaatttctgaatgctcggccaggaggagatgatcagagagacactaatgccacccaacatcaatagcacaatcagccccgccagttttaccagccagctgaacattctgtcgaccttt      >  2:26743‑M1/1‑14 (MQ=255)
                                                                                                                                                                                                                                                                                      aTCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgcccagtagatatcgttcggtgcatcccactctttggtccataggaaagccagaccaaatttctgaatgctcggcccggaggagatgatcagagagacaataatgccacccaacatcaatagcacaatcagccccgccagttttaccagcacgctgaaaattaagtccccatttttac  >  2:109736‑M1/1‑9 (MQ=255)
                                                                                                                                                                                                                                                                                      aTCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactaaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtccacaggaaagccagaccaaatttctgaatgctcggccaggaggagatgatcagagagacaataatgccacacaacatcaatagcccaatcagcgccgcaagttttaccagaacgctgaaaataatgtagccctttttac  >  1:244222‑M1/1‑9 (MQ=255)
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                       
ACAGACCAACCACAATCGACGGCGCAGAGAGCAGAATGTCGTTAATGAAGCGAATCACTTCTGCCAGCCAGGATTTACGACCATATTCCGCCAGATAAATCCCCGCCATAATGCCCAGCGGCGTACCGAATACCGTGGCCCACAAAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCG  >  NC_000913/3909024‑3909589

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