Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,112,162 A→C pseudogene (780/849 nt) yghF ← pseudogene, secretion pathway protein, C‑type protein homology;putative transport; Not classified; putative general secretion pathway for protein export (GSP)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,112,1620AC100.0% 208.0 / NA 57pseudogene (780/849 nt)yghFpseudogene, secretion pathway protein, C‑type protein homology;putative transport; Not classified; putative general secretion pathway for protein export (GSP)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (27/30);  total (27/30)

ATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAACCGACAAT  >  NC_000913/3112016‑3112305
                                                                                                                                                  |                                                                                                                                               
atTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGt                                                                                                                                              >  1:886278/1‑150 (MQ=255)
 ttCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTa                                                                                                                                             <  1:932248/150‑1 (MQ=255)
    aaCGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTg                                                                                                                                         <  1:182400/151‑1 (MQ=255)
         ggCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCa                                                                                                                                    <  2:590056/151‑1 (MQ=255)
              aaTGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATc                                                                                                                                  >  1:678255/1‑148 (MQ=255)
              aaTGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCaga                                                                                                                               <  2:123533/151‑1 (MQ=255)
               aTGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCagag                                                                                                                              <  2:349299/151‑1 (MQ=255)
                 gCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCa                                                                                                                            <  2:742596/151‑1 (MQ=255)
                          cGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGCAACTGCAGCATCAGAGCAATCATTGc                                                                                                                    >  1:498225/1‑150 (MQ=255)
                           gATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTc                                                                                                                  >  2:887176/1‑151 (MQ=255)
                            atatCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCg                                                                                                                 <  2:404065/151‑1 (MQ=255)
                              atCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTg                                                                                                               >  2:300813/1‑151 (MQ=255)
                                   ggTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGAt                                                                                                            <  2:635848/149‑1 (MQ=255)
                                   ggTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCa                                                                                                          <  2:305677/151‑1 (MQ=255)
                                      tCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTg                                                                                                       >  1:257634/1‑151 (MQ=255)
                                       cAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTg                                                                                                       <  2:445420/150‑1 (MQ=255)
                                       cAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGa                                                                                                      <  2:513936/151‑1 (MQ=255)
                                        aGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGaa                                                                                                     <  1:895174/151‑1 (MQ=255)
                                           tcttctTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAAt                                                                                                   >  1:465710/1‑150 (MQ=255)
                                            cttcttCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATcc                                                                                                 >  1:667578/1‑151 (MQ=255)
                                               cttcCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCctgc                                                                                              >  2:227080/1‑151 (MQ=255)
                                                ttcCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCctgct                                                                                             >  1:407443/1‑151 (MQ=255)
                                                 tcCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCctgctg                                                                                            >  1:291638/1‑151 (MQ=255)
                                                 tcCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCctgct                                                                                             <  1:467604/150‑1 (MQ=255)
                                                       aaGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGAttt                                                                                        >  2:579112/1‑149 (MQ=255)
                                                         ggCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAgcg                                                                                    >  2:612939/1‑151 (MQ=255)
                                                            gggCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCa                                                                                  >  2:518911/1‑150 (MQ=255)
                                                                  aaTAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGa                                                                           <  1:303295/151‑1 (MQ=255)
                                                                      gTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATc                                                                       >  2:398242/1‑151 (MQ=255)
                                                                      gTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATc                                                                       >  2:912626/1‑151 (MQ=255)
                                                                         ggTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGcc                                                                    >  1:255007/1‑151 (MQ=255)
                                                                            agcgcagcgcGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGcctt                                                                  <  1:675731/150‑1 (MQ=255)
                                                                             gcgcagcgcGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGccttcc                                                                >  1:707443/1‑151 (MQ=255)
                                                                             gcgcagcgcGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGcctt                                                                  <  1:523880/149‑1 (MQ=255)
                                                                              cgcagcgcGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGccttcct                                                               >  2:427/1‑151 (MQ=255)
                                                                              cgcagcgcGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGccttcct                                                               >  2:460665/1‑151 (MQ=255)
                                                                                 agcgcGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTg                                                             >  1:373978/1‑150 (MQ=255)
                                                                                         ggAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCt                                                     <  2:598935/150‑1 (MQ=255)
                                                                                           aaaTGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGc                                                  <  2:835100/151‑1 (MQ=255)
                                                                                                  gTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAAc                                           >  2:387378/1‑151 (MQ=255)
                                                                                                      cgcgcACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACagag                                       <  1:286737/151‑1 (MQ=255)
                                                                                                        cgcACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAgaga                                      >  2:66036/1‑150 (MQ=255)
                                                                                                         gcACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAAc                                    <  2:625405/151‑1 (MQ=255)
                                                                                                          cACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGaa                                     <  1:925432/149‑1 (MQ=255)
                                                                                                            cccTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACg                                   <  2:713594/149‑1 (MQ=255)
                                                                                                            cccTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGa                                  <  2:669578/150‑1 (MQ=255)
                                                                                                                tGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTg                              <  2:464682/150‑1 (MQ=255)
                                                                                                                tGCGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGc                             >  1:250241/1‑151 (MQ=255)
                                                                                                                  cGTAAAACCGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGccc                           <  1:225867/151‑1 (MQ=255)
                                                                                                                          cGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTc                   <  2:587819/151‑1 (MQ=255)
                                                                                                                          cGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTc                   <  1:417822/151‑1 (MQ=255)
                                                                                                                          cGTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTc                   <  2:153013/151‑1 (MQ=255)
                                                                                                                           gTAAGTTGAATGGAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCg                  <  2:727376/151‑1 (MQ=255)
                                                                                                                                       gAATCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAACCg       >  2:200539/1‑150 (MQ=255)
                                                                                                                                        aaTCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAACCGAc     >  1:454196/1‑151 (MQ=255)
                                                                                                                                        aaTCCATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAACCGAc     >  1:769802/1‑151 (MQ=255)
                                                                                                                                           ccATTGACGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAACCGACAAt  <  2:50970/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAACCGACAAT  >  NC_000913/3112016‑3112305

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