Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,372,569 G→A 6.6% intergenic (‑1154/‑418) PP_5636 ← / → PP_3849 hypothetical protein/hemolysin‑type calcium‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,372,5690GA 6.6% 173.8 / 2.3 60intergenic (‑1154/‑418)PP_5636/PP_3849hypothetical protein/hemolysin‑type calcium‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (30/26);  new base A (2/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

CATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGAACATCCATTGTAAATGGCGGGGTTAATGACTCTTGGA  >  NC_002947/4372422‑4372711
                                                                                                                                                   |                                                                                                                                              
cATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACgg                                                                                                                                               >  7:142755/1‑149 (MQ=255)
                 cGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGAtt                                                                                                                              >  5:106379/1‑149 (MQ=255)
                         aTTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGAt                                                                                                                               >  8:202346/1‑140 (MQ=255)
                         aTTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGAt                                                                                                                               <  7:202346/140‑1 (MQ=255)
                               ttATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCGCCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATc                                                                                                                   >  1:220730/1‑146 (MQ=255)
                                    gATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACAt                                                                                                           >  7:190361/1‑149 (MQ=255)
                                      tAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCACTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATc                                                                                                          <  8:142755/148‑1 (MQ=255)
                                        aaTCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTa                                                                                                       <  7:223826/149‑1 (MQ=255)
                                        aaTCATCAATCAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTa                                                                                                       <  4:73780/149‑1 (MQ=255)
                                          tcatcaATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAgt                                                                                                     <  6:250427/149‑1 (MQ=255)
                                           catcaATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAgtg                                                                                                    >  4:248712/1‑149 (MQ=255)
                                             tcaATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGcc                                                                                                                                        >  1:24548/1‑111 (MQ=255)
                                             tcaATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGcc                                                                                                                                        <  2:24548/111‑1 (MQ=255)
                                             tcaATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAgtgt                                                                                                   >  3:196375/1‑148 (MQ=255)
                                               aaagAGCCCCGCGCCAGACACTACCGCACCCAAACTGTCCACACTAACAATACCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCAt                                                                                                <  6:370816/146‑1 (MQ=255)
                                               aaTGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCAt                                                                                                >  1:184849/1‑149 (MQ=255)
                                                 tGAGACACGCTCCAGACACTACCGCACCCAAATTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCAttt                                                                                              <  8:149092/149‑1 (MQ=255)
                                                     ccACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCAGGTCTTACAAAACACGGCCAAGCCACGGACGATTTATATGGCATCAAGAACATCGTAGTGTCATTTGa                                                                                            >  7:145767/1‑147 (MQ=255)
                                                     ccACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGAAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAACCCAAGGACGATTTAGCTGGCATCAAGAACATCGATGTGTCATTTGa                                                                                            >  6:291046/1‑147 (MQ=255)
                                                        cGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACAGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGa                                                                                        >  1:2622/1‑148 (MQ=255)
                                                        cGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACAGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGa                                                                                        >  1:8670/1‑148 (MQ=255)
                                                             cAGACACTACCGCACCCAAACTGTCATCCCTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCCCGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTg                                                                                  >  2:228491/1‑149 (MQ=255)
                                                               gACACTACCGCCCCCAAACTGGCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGc                                                                                <  8:352719/149‑1 (MQ=255)
                                                               gACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGGCATTTGACAGAAATCTGGc                                                                                >  3:44445/1‑149 (MQ=255)
                                                                  acTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCtt                                                                              >  8:213821/1‑148 (MQ=255)
                                                                      ccGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACATTAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTc                                                                          <  6:142456/148‑1 (MQ=255)
                                                                         cACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCata                                                                       >  3:100630/1‑148 (MQ=255)
                                                                            ccAAACTGTCATCACTAACAATTACCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATac                                                                    >  7:51589/1‑148 (MQ=255)
                                                                              aaaCTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACt                                                                 <  8:1805/149‑1 (MQ=255)
                                                                              aaaCTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACt                                                                 <  2:220730/149‑1 (MQ=255)
                                                                               aaCTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTc                                                                          >  4:39059/1‑139 (MQ=255)
                                                                               aaCTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTc                                                                          <  3:39059/139‑1 (MQ=255)
                                                                                    ttatgaCTAACAATTCCCTCGATAGTCAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCTAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATg                                                            <  4:44445/143‑1 (MQ=255)
                                                                                      atcaCTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACa                                                         >  6:241906/1‑149 (MQ=255)
                                                                                          cTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTa                                                     >  7:124674/1‑149 (MQ=255)
                                                                                            aaCAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACg                                                   >  3:37848/1‑149 (MQ=255)
                                                                                             aCAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAAt                                                             <  3:12260/138‑1 (MQ=255)
                                                                                             aCAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAAt                                                             >  4:12260/1‑138 (MQ=255)
                                                                                             aCAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGc                                                  >  1:59216/1‑149 (MQ=255)
                                                                                                   cccTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCgaga                                             <  8:190361/148‑1 (MQ=255)
                                                                                                   cccTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCgag                                              >  2:244533/1‑147 (MQ=255)
                                                                                                     cTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGagaaga                                          >  8:198518/1‑149 (MQ=255)
                                                                                                        gATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACg                                        <  7:198518/148‑1 (MQ=255)
                                                                                                         aTAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGaa                                      <  4:196375/149‑1 (MQ=255)
                                                                                                          tAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGAAc                                     <  7:213821/149‑1 (MQ=255)
                                                                                                             tAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGAACATc                                  >  6:332094/1‑149 (MQ=255)
                                                                                                               aaTCTTTCGTCAGGGCGGGCTTGTCTTACAATAGACGGCCAAGCGACGGAGGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGAACATcc                                 <  8:145767/148‑1 (MQ=255)
                                                                                                                     tCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGAACATCCATTGt                            <  5:332094/147‑1 (MQ=255)
                                                                                                                           gggcgggcTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGAACATCCATTGTAAATGGCg                    <  4:37848/149‑1 (MQ=255)
                                                                                                                           gggcgggcTGGTCTTACAATAGACAGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGAACATCCATTGTAAATGGCg                    <  2:8670/149‑1 (MQ=255)
                                                                                                                           gggcgggcTGGTCTTACAATAGACAGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGAACATCCATTGTAAATGGCg                    <  2:2622/149‑1 (MQ=255)
                                                                                                                            ggcgggcTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGAACATCCATTGTAAATGGCgg                   <  3:248712/149‑1 (MQ=255)
                                                                                                                              cgggcTGGTCTTACAATAGACGGCCACGCCCCGGACGATTTATCTGGCATCAAGAACATCGTAGTGGCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGAACATCAATTGGAAAGGGCgggg                 >  2:117363/1‑149 (MQ=255)
                                                                                                                                 gcTGGACTTACAATAGACGGACAAGCCACGGACGATTTATCTTGCATCAAGAACATCGTAGTCTCATTTGACAGAAATCTGGCTAACTAATATACACTGAAAGACACTTACGCCGAGAAGACGAACATCCATAGTAAATGGAGGGGtt               >  8:281757/1‑148 (MQ=255)
                                                                                                                                   tGGTCTTACAATAGACGGCCAAGCCACTGACGATTTAGATGGCATCAAGAACATCGTTGTGTCATTTGACAGTTTTCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGCCGAACCTCCATTGTAAATGGCGGGGTTAAt            <  1:117363/149‑1 (MQ=255)
                                                                                                                                   tGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGATGAACATCCATTGTAAATGGCGGGGTTAAt            >  7:166633/1‑149 (MQ=255)
                                                                                                                                           cAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGAACATCCATTGTAAATGGCGGGGTTAATGACTCtt     <  4:100630/148‑1 (MQ=255)
                                                                                                                                           cAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGAACATCCATTGTAAATGGCGGGGTTAATGACTCTTg    <  2:184849/149‑1 (MQ=255)
                                                                                                                                            aaTAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGAACATCCATTGTAAATGGCGGGGTTAATGACTCTTg    >  6:28554/1‑148 (MQ=255)
                                                                                                                                              tAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGAACATCCATTGTAAATGGCGGGGTTAATGACTCTTGGa  >  8:212519/1‑148 (MQ=255)
                                                                                                                                                gACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGAACATCCATTGTAAATGGCGGGGTTAATGACTCTTGGa  >  5:304142/1‑146 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCTGGCATCAAGAACATCGTAGTGTCATTTGACAGAAATCTGGCTTACTCATATACACTGAATGACAGTTACGCCGAGAAGACGAACATCCATTGTAAATGGCGGGGTTAATGACTCTTGGA  >  NC_002947/4372422‑4372711

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