Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,034,536 G→A pseudogene (486/663 nt) yedS → pseudogene, outer membrane protein homology, putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,034,5360GA100.0% 241.7 / NA 74pseudogene (486/663 nt)yedSpseudogene, outer membrane protein homology, putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/36);  total (38/36)

GTCGTATATGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAGGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTATGTCAACA  >  NC_000913/2034402‑2034671
                                                                                                                                      |                                                                                                                                       
gtCGTATATGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                                     <  1:1328086/139‑1 (MQ=255)
         gACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGa                                                                                                                            >  2:448502/1‑139 (MQ=255)
           cGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAAc                                                                                                                          >  1:1043082/1‑139 (MQ=255)
           cGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAAc                                                                                                                          >  1:837012/1‑139 (MQ=255)
           cGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAAc                                                                                                                          >  1:440453/1‑139 (MQ=255)
           cGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAAc                                                                                                                          >  1:1218115/1‑139 (MQ=255)
           cGTAGAAGCATGGACCCATATACTACCAGAATTTGTTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACTATGGCTTCTTTGGTTAAGTCGATGGTTTGAAc                                                                                                                          >  2:565051/1‑139 (MQ=255)
                aaGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGa                                                                                                                                    <  2:1131568/124‑1 (MQ=255)
                aaGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGa                                                                                                                                    >  1:1131568/1‑124 (MQ=255)
                aaGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgc                                                                                                                     <  2:1218115/139‑1 (MQ=255)
                aaGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgc                                                                                                                     <  1:887328/139‑1 (MQ=255)
                aaGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActt                                                                                                                        >  2:266523/1‑136 (MQ=255)
                aaGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActt                                                                                                                        <  1:266523/136‑1 (MQ=255)
                     tGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCCGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgc                                                                                                                     >  2:1534684/1‑134 (MQ=255)
                     tGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCCGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgc                                                                                                                     <  1:1534684/134‑1 (MQ=255)
                                gctgcCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGt                                                                                                     <  1:1228074/139‑1 (MQ=255)
                                    ccAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGtaata                                                                                                 >  1:487180/1‑139 (MQ=255)
                                        aaTTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGa                                                                                             <  2:837012/139‑1 (MQ=255)
                                               tggCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAgtggt                                                                                      <  2:1043082/139‑1 (MQ=255)
                                                  cGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGa                                                                                   >  2:1178352/1‑139 (MQ=255)
                                                     tACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGAt                                                                                                                                   <  2:538851/88‑1 (MQ=255)
                                                     tACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGAt                                                                                                                                   >  1:538851/1‑88 (MQ=255)
                                                           tGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgc                                                                                                                     <  2:711126/96‑1 (MQ=255)
                                                           tGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgc                                                                                                                     >  1:711126/1‑96 (MQ=255)
                                                            gCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGc                                                                         <  1:448502/139‑1 (MQ=255)
                                                                gcgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGaa                                                                                                                           >  1:208137/1‑85 (MQ=255)
                                                                gcgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGaa                                                                                                                           <  2:208137/85‑1 (MQ=255)
                                                                gcgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGa                                                                     >  2:1179361/1‑139 (MQ=255)
                                                                 cgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttg                                                                                                                      <  2:203072/89‑1 (MQ=255)
                                                                 cgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttg                                                                                                                      >  1:203072/1‑89 (MQ=255)
                                                                 cgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGaa                                                                    >  2:1078001/1‑139 (MQ=255)
                                                                           ttCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTa                                                          <  1:1178352/139‑1 (MQ=255)
                                                                            tCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAg                                                         <  1:1179361/139‑1 (MQ=255)
                                                                                  aCGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                                     <  2:252371/57‑1 (MQ=255)
                                                                                  aCGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                                     >  1:252371/1‑57 (MQ=255)
                                                                                     gTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAatat                                                                                                           <  2:815793/80‑1 (MQ=255)
                                                                                     gTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAatat                                                                                                           >  1:815793/1‑80 (MQ=255)
                                                                                     gTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATgg                                                >  2:811565/1‑139 (MQ=255)
                                                                                           cTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAg                                          >  1:1424286/1‑139 (MQ=255)
                                                                                           cTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAg                                          >  2:177702/1‑139 (MQ=255)
                                                                                              aTGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTTAACTTTGCTTTGCAATCTCa                                                                                                         <  1:1042080/73‑1 (MQ=38)
                                                                                              aTGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCa                                                                                                         >  2:1042080/1‑73 (MQ=255)
                                                                                                ggTGTTGCTACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGtttcattat                                                                                             >  1:996505/1‑75 (MQ=255)
                                                                                                ggTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGtttcattat                                                                                             <  2:996505/83‑9 (MQ=255)
                                                                                                ggTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCttt                                     >  2:1336595/1‑139 (MQ=255)
                                                                                                    ttGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAgtgg                                                                                       >  2:616734/1‑85 (MQ=255)
                                                                                                    ttGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAgtgg                                                                                       <  1:616734/85‑1 (MQ=255)
                                                                                                    ttGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCAtt                                                                              <  1:99577/94‑1 (MQ=255)
                                                                                                    ttGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCAtt                                                                              >  2:99577/1‑94 (MQ=255)
                                                                                                      gCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAg                               >  1:602149/1‑139 (MQ=255)
                                                                                                            tACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGa                                                                     >  2:1061719/1‑95 (MQ=255)
                                                                                                            tACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGa                                                                     <  1:1061719/95‑1 (MQ=255)
                                                                                                                 tAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAatat                                                                                                           <  1:974042/52‑1 (MQ=255)
                                                                                                                 tAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAatat                                                                                                           >  2:974042/1‑52 (MQ=255)
                                                                                                                    cAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGa                                                   <  2:1407353/105‑1 (MQ=255)
                                                                                                                    cAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGa                                                   >  1:1407353/1‑105 (MQ=255)
                                                                                                                     aaTGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCaaa                                <  1:511845/123‑1 (MQ=255)
                                                                                                                     aaTGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCaaa                                >  2:511845/1‑123 (MQ=255)
                                                                                                                        ggCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAatat                                                                                                           >  1:755807/1‑45 (MQ=255)
                                                                                                                        ggCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAatat                                                                                                           <  2:755807/45‑1 (MQ=255)
                                                                                                                         gCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAatat                                                                                                           >  1:483870/1‑44 (MQ=255)
                                                                                                                         gCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAatat                                                                                                           <  2:483870/44‑1 (MQ=255)
                                                                                                                         gCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCgg            <  1:1336595/139‑1 (MQ=255)
                                                                                                                         gCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCgg            >  2:1447029/1‑139 (MQ=255)
                                                                                                                          cttcttTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGt           <  2:487180/139‑1 (MQ=255)
                                                                                                                          cttcttTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGt           <  2:440453/139‑1 (MQ=255)
                                                                                                                           ttcttTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAatat                                                                                                           <  1:242306/42‑1 (MQ=255)
                                                                                                                           ttcttTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAatat                                                                                                           >  2:242306/1‑42 (MQ=255)
                                                                                                                              tttGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTATGt       <  2:1424286/139‑1 (MQ=255)
                                                                                                                                 ggTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGa                                                                                             >  1:196507/1‑50 (MQ=255)
                                                                                                                                 ggTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGa                                                                                             <  2:196507/50‑1 (MQ=255)
                                                                                                                                 ggTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTATGTCaa    <  1:122297/139‑1 (MQ=255)
                                                                                                                                 ggTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTATGTCaa    <  1:811565/139‑1 (MQ=255)
                                                                                                                                   ttAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTATGTCAACa  >  1:1059654/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                       
GTCGTATATGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAGGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTATGTCAACA  >  NC_000913/2034402‑2034671

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