Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,466,540 9636 bp→ACCC intS[ypdJ] 15 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2466540 2476175 9636 39 [0] [0] 41 intS–[ypdJ] intS,yfdG,yfdH,yfdI,yodE,yfdK,yfdM,yfdN,yfdO,yfdP,yfdQ,yfdR,yfdS,[yfdT],[ypdJ]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 24665390 (0.000)39 (0.500)
+ACCC
38/486 0.6 100% intergenic (+156/‑6) argW/intS tRNA‑Arg/prophage CPS‑53 integrase
?NC_000913 2476176 = 0 (0.000)coding (304/306 nt)
coding (83/141 nt)
yfdT
ypdJ
protein YfdT
putative uncharacterized protein YpdJ

ACTCATGCCGCATTGTCCTCTTAGTTAAATGGATATAACGAGCCCCTCCTAAGGGCTAATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTCCCTGTGTCAGGGAAAAAGGTATGGCAGCTTCGCTACAAGATTGACGGTAAGGAGAAAATCCTGACCGTCGGAAAATATCCGCTTATGACTTTGCAGGAGGCAAGGGATAAAGCATGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAA  >  NC_000913/2466295‑2466788
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         
cctCATGCCGCATTGTCCTCTTAGTTAAATGGATATAACGAGCCCCTCCTAAGGGCTAATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATACAATAATTTTACCCACCGGATTTTacccta                                                                                                                                                                                                                                                     <  2:352101‑M1/250‑7 (MQ=255)
    aTGCCGCATTGTCCTCTTAGTTAAATGGATATAACGAGCCCCTCCTAAGGGCTAATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatag                                                                                                                                                                                                                                                 <  1:622417‑M1/251‑11 (MQ=255)
                    ttAGTTAAATGGATATAACGAGCCCCTCCTAAGGGCTAATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctaac                                                                                                                                                                                                                                 >  1:401404‑M1/1‑225 (MQ=255)
                      aGTTAAATGGATATAACGAGCCCCTCCTAAGGGCTAATTGCAGGTTCGAGTCCTTCAGGGGACTCCATGTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGTTTCTCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTGACCGACCGGATTTTgccctgatagagccatacagcctcacac                                                                                                                                                                                                                               <  2:484535‑M1/251‑29 (MQ=255)
                            aTGGATATAACGAGCCCCTCCTAAGGGCTAATTGCAGGTTCGATTCCTGCAGGGGACGCCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatg                                                                                                                                                                                                                         <  1:575991‑M1/251‑35 (MQ=255)
                                   taACGAGCCCCTCCTAAGGGCTAATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctc                                                                                                                                                                                                                                   >  2:63425‑M1/1‑210 (MQ=255)
                                   taACGAGCCCCTCCTAAGGGCTAATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctc                                                                                                                                                                                                                                   <  1:63425‑M1/234‑25 (MQ=255)
                                    aaCGAGCCCCTCCTAAGGGCTAATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgt                                                                                                                                                                                                                 >  1:403909‑M1/1‑209 (MQ=255)
                                       gAGCCCCTCCTAAGGGCTAATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtc                                                                                                                                                                                                                >  2:204959‑M1/1‑206 (MQ=255)
                                          ccccTCCTAATGGCTAATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTCCCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagttt                                                                                                                                                                                                           <  1:540065‑M1/251‑49 (MQ=255)
                                                   aGGGCTAATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaat                                                                                                                                                                                                   >  2:203981‑M1/1‑194 (MQ=255)
                                                         aaTTGCAGTTTCGATTCCTGCAGGTGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCATTGACAGCCGCACGCCATGTCGGGTAAAAAGTGGATAAAGTAAGTGTACCCACTGGATTTGaccctaatagagccatacagcctcactctcgatgaggtctgtcagtttcttaaaatagcctga                                                                                                                                                                                            <  2:189217‑M1/251‑64 (MQ=255)
                                                            tGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatc                                                                                                                                                                                          <  1:345420‑M1/250‑66 (MQ=255)
                                                                      aTTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccg                                                                                                                                                                               >  1:143995‑M1/1‑175 (MQ=255)
                                                                         ccTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtc                                                                                                                                                                             <  1:428036‑M1/250‑79 (MQ=255)
                                                                                              cAGTTCGCTCCCATCCGTCCCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttattacct                                                                                                                                                       <  2:365595‑M1/251‑101 (MQ=255)
                                                                                              cAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcggatttttattacct                                                                                                                                                       >  2:303091‑M1/1‑151 (MQ=255)
                                                                                                       cccATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttattacctgatttgc                                                                                                                                                >  1:425440‑M1/1‑142 (MQ=255)
                                                                                                                  cAGTCCGCAAAATCCCCTGCATATCAAGCATCCCGTAGATTTACAGTACGTCATGGACCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGAATTTTaccctaatagagccatacagcctcacactcgatgaggcctggcagtttcttaaaatatcctgagctaccatcgccgtcatagagcgtatttttatttccggatttgcatgttcgattcca                                                                                                                                   <  2:53356‑M1/251‑121 (MQ=255)
                                                                                                                      ccGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttattacctgatttgcaggttcgattccctatt                                                                                                                               >  1:586426‑M1/1‑127 (MQ=255)
                                                                                                                            aaTCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttattaccttcttttcaggttcgattccctattcgg                                                                                                                            >  2:594‑M1/1‑121 (MQ=255)
                                                                                                                                       tatCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttattacctgatttgcaggttcgattccctattcggagatagcactcatg                                                                                                              >  1:562444‑M1/1‑110 (MQ=255)
                                                                                                                                                          ttACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttattacctgatttgcaggttcgattccctattcggagatagcactcatgcaacacgaactacagcct                                                                                            <  1:594‑M1/250‑160 (MQ=255)
                                                                                                                                                            aCAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttatcacctgatttgcaggttcgattccctattcggagatagcactcatgcaacacgaacgacagcctgat                                                                                         >  2:374645‑M1/1‑89 (MQ=255)
                                                                                                                                                                  cGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAACAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttattacctgatttgcaggttcgattccctattcggagatagcactcatgcaacacgaactacagcctgattcactg                                                                                   >  1:549751‑M1/1‑83 (MQ=255)
                                                                                                                                                                         gTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttattacctgatttgcaggttcgattccctattcggagatagcactcatgcaacacgaactacagcctgattcactggttgatt                                                                            <  1:294645‑M1/251‑176 (MQ=255)
                                                                                                                                                                                       gTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTacccttatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttattacctgatttgcagtttcgattccctattcggagatagcactcatgcaacacgaactacagcctgattcactggttgatttgaaattcatcatg                                                              <  2:528917‑M1/251‑190 (MQ=255)
                                                                                                                                                                                         tGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCTCCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttattacctgatttgcaggttcgattccctattcggagatagcactcatgcaacacgaactacagcctgattcactggttgatttgaaattcagcatgg                                                             <  1:344608‑M1/250‑191 (MQ=255)
                                                                                                                                                                                                 gCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatcttccatcgccttcatagagcgtatttttattacctgatttgcaggttcgattccctattcggagatagcactcatgcaacacgaactacagcctgattcgctggttgatttgaaattcatcatggctgatactg                                                    <  1:455032‑M1/251‑200 (MQ=255)
                                                                                                                                                                                                      ccATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctactagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagcgcgtatttttattacctgatttgcaggttcgattccctattcggagatagcactcatgcaacacgaactacagcctgattcactggttgatttgaaattcatcatggctgatactggcttt                                               >  1:515785‑M1/1‑47 (MQ=255)
                                                                                                                                                                                                       cATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttattacctgatttgcaggttcgattccctattcggagatagcactcatgcaacacgaactacagcctgattcactggttgatttgaaattcatcatggctgatactggctttg                                              <  2:277990‑M1/251‑206 (MQ=255)
                                                                                                                                                                                                            cGGGTAAAAAGTGGATAAAATAATTTTACACACCGGCTTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcttagagcgtatttttattacctgatttgcaggttcgattccctattcggagatagcactcatgcaacacgaactacagcctgattcactggttgatttgaatttcatcatggctgatactggctttcgtaa                                          >  2:381992‑M1/1‑41 (MQ=255)
                                                                                                                                                                                                                  aaaaGTGGATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttattacctgatttgcaggttcgattccctattcggagatagcactcatgcaacacgaactacagcctgattcactggttgatttgaaattcatcatggctgatactggctttggtaaaaccttc                                   <  1:388961‑M1/251‑217 (MQ=255)
                                                                                                                                                                                                                        ggATAAAATAATTTTACCCACCGGATTTTaccctaatagagccatacagccgcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatctccgtcatagagcgtatttttattacctgatttgcaggttcgattccctatgcggagatatcactcatgcaacacgaactacagcctgattcactggttgatttgaaattcatcatggctgatactggctttggtaaaaccttcatctat                             <  2:401404‑M1/251‑223 (MQ=255)
                                                                                                                                                                                                                                   ttttACCAACCGGATTTTaccctaatagagcaatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttattacctgatttgcaggttcgattccctattcggagatagcactcatgcaacacgaactacagcctgattcactggttgatttgcaattcatcatggctgatactggcattggtaaaaccttcatctatgaccggatta                   >  2:148752‑M1/1‑18 (MQ=255)
                                                                                                                                                                                                                                        cccACCGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcgtaaaagatcctgatctaccatcgccgtcatagagcgtagttttattacctgatttgcaggttcgattccctattcggagatagcactcatgcaacacgaactacagcctgattcactggttgatttgaaatacatcatggctgatactggcgttggtaaaaccttcatctatgaccggattaagtcag             <  2:577895‑M1/251‑239 (MQ=255)
                                                                                                                                                                                                                                             cGGATTTTaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttattacctgatttgcaggttcgattccctattcggagatagcactcatgcaacacgaactacagcctgattcactggttgatttgaaattcatcatggctgatactggctttggtaaaaccttcatctatgaccggattaagtcaggcgac        <  2:143995‑M1/251‑244 (MQ=255)
                                                                                                                                                                                                                                                 ttttaccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttattacctgatttgcaggttcgattccctattcggagatagcactcatgcaacacgaactacagcctgattcactggttgatttgaaattcatcatggctgatactggctttggtaaaaccttcatctatgaccggattaagtcaggcgacctgc    >  1:341865‑M1/1‑4 (MQ=255)
                                                                                                                                                                                                                                                    taccctaatagagccatacagcctcacactcgatgaggcctgtcagtttcttaaaatatcctgatctaccatcgccgtcatagagcgtatttttattacctgatttgcaggttcgattccctattcggagatagcactcatgcaacacgaactacagcctgattcactggttgatttgaaattcatcatggctgatactggctttggtaaaaccttcatctatgaccggattaagtcaggcgacctgcca  <  1:174262‑M1/250‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         
ACTCATGCCGCATTGTCCTCTTAGTTAAATGGATATAACGAGCCCCTCCTAAGGGCTAATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTCCCTGTGTCAGGGAAAAAGGTATGGCAGCTTCGCTACAAGATTGACGGTAAGGAGAAAATCCTGACCGTCGGAAAATATCCGCTTATGACTTTGCAGGAGGCAAGGGATAAAGCATGGACTGCGAGGAAAGACATCTCGGTTGGCATCGATCCTGTAA  >  NC_000913/2466295‑2466788

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