Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,041,159 G→A 8.6% intergenic (+41/‑216) yedZ → / → zinT inner membrane heme subunit for periplasmic YedYZ reductase/zinc and cadmium binding protein, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,041,1590GA 8.6% 219.2 / 9.9 70intergenic (+41/‑216)yedZ/zinTinner membrane heme subunit for periplasmic YedYZ reductase/zinc and cadmium binding protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base A (3/3);  total (35/35)
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.80e-01

TCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCT  >  NC_000913/2041012‑2041299
                                                                                                                                                   |                                                                                                                                            
tCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGAtgt                                                                                                                                             <  2:287148/149‑1 (MQ=255)
  aTCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGACGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTgg                                                                                                                                           >  2:359408/1‑149 (MQ=255)
            ggctggctGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGcc                                                                                                                                 <  2:1331710/149‑1 (MQ=255)
                           ttttAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGcc                                                                                                                                 <  2:668999/134‑1 (MQ=255)
                           ttttAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGcc                                                                                                                                 >  1:668999/1‑134 (MQ=255)
                                  cTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGc                                                                                                                        >  1:1086251/1‑136 (MQ=255)
                                  cTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGc                                                                                                                        <  2:1086251/136‑1 (MQ=255)
                                      cGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCa                                                                                                       <  1:359408/149‑1 (MQ=255)
                                           taAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAgg                                                                                                     <  2:145647/146‑1 (MQ=255)
                                           taAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAgg                                                                                                     >  1:145647/1‑146 (MQ=255)
                                             agaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGcc                                                                                                                                 >  3:356008/1‑116 (MQ=255)
                                             agaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGcc                                                                                                                                 <  4:356008/116‑1 (MQ=255)
                                              gaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAg                                                                                                      <  1:135104/142‑1 (MQ=255)
                                              gaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAg                                                                                                      >  2:135104/1‑142 (MQ=255)
                                                 gTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGAt                                                                                                              <  1:868165/131‑1 (MQ=255)
                                                 gTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGAt                                                                                                              >  2:1350959/1‑131 (MQ=255)
                                                 gTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGAt                                                                                                              >  2:868165/1‑131 (MQ=255)
                                                 gTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGAt                                                                                                              <  1:1350959/131‑1 (MQ=255)
                                                  ttGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGgcgc                                                                                                                                        >  3:757098/1‑104 (MQ=255)
                                                  ttGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGgcgc                                                                                                                                        <  4:757098/104‑1 (MQ=255)
                                                  ttGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGgc                                                                                                                                          >  3:6238/1‑102 (MQ=255)
                                                  ttGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGgc                                                                                                                                          <  4:481525/102‑1 (MQ=255)
                                                  ttGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGgc                                                                                                                                          <  4:6238/102‑1 (MQ=255)
                                                  ttGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGgc                                                                                                                                          >  3:481525/1‑102 (MQ=255)
                                                      gTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTg                                                                                       >  4:141086/1‑149 (MQ=255)
                                                       ttCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGt                                                                                      >  2:1084437/1‑149 (MQ=255)
                                                          ccTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGcc                                                                                                                                 >  4:1209776/1‑103 (MQ=255)
                                                          ccTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGcc                                                                                                                                 <  3:1209776/103‑1 (MQ=255)
                                                          ccTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGtcat                                                                                   <  3:141086/149‑1 (MQ=255)
                                                           cTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGAt                                                                                                              >  3:568128/1‑121 (MQ=255)
                                                           cTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGAt                                                                                                              <  4:568128/121‑1 (MQ=255)
                                                           cTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGccacat                                                                                                                             <  1:455273/106‑5 (MQ=255)
                                                           cTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGccacat                                                                                                                             >  2:455273/1‑102 (MQ=255)
                                                                     ggTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAAt                                                                        >  2:126949/1‑149 (MQ=255)
                                                                      gTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATcc                                                                                                                           <  3:1333234/97‑1 (MQ=255)
                                                                      gTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATcc                                                                                                                           >  4:1333234/1‑97 (MQ=255)
                                                                               gCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAaata                                                              <  1:126949/149‑1 (MQ=255)
                                                                                cAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAaataa                                                             <  1:1084437/149‑1 (MQ=255)
                                                                                         aacaaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAg                                                                                                                <  2:1269576/89‑1 (MQ=255)
                                                                                         aacaaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAg                                                                                                                >  1:1269576/1‑89 (MQ=255)
                                                                                         aacaaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCag                                                                                                                      >  1:1061961/1‑82 (MQ=21)
                                                                                         aacaaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCag                                                                                                                      <  2:1061961/83‑2 (MQ=21)
                                                                                           caaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCaaa                                                                         >  2:1325/1‑126 (MQ=255)
                                                                                           caaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCaaa                                                                         <  1:1325/126‑1 (MQ=255)
                                                                                            aaaTTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGt                                                                                               >  4:1319635/1‑103 (MQ=255)
                                                                                            aaaTTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGt                                                                                               <  3:1319635/103‑1 (MQ=255)
                                                                                                      tataACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGtt                                                                                                   <  2:67263/89‑1 (MQ=255)
                                                                                                      tataACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGtt                                                                                                   >  1:67263/1‑89 (MQ=255)
                                                                                                      tataACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGtcatc                                                                                  >  2:768828/1‑106 (MQ=255)
                                                                                                      tataACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGtcatc                                                                                  <  1:768828/106‑1 (MQ=255)
                                                                                                      tataACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAAc                                                                               >  3:377834/1‑109 (MQ=255)
                                                                                                      tataACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAAc                                                                               <  4:377834/109‑1 (MQ=255)
                                                                                                      tataACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAaataat                                                            >  1:1284812/1‑128 (MQ=255)
                                                                                                      tataACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAaataat                                                            <  2:1284812/128‑1 (MQ=255)
                                                                                                      tataACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCag                                                                                                                      <  4:783851/70‑2 (MQ=17)
                                                                                                      tataACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCag                                                                                                                      >  3:783851/1‑69 (MQ=18)
                                                                                                           cTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAAtttatttta                                                    <  1:14556/131‑1 (MQ=255)
                                                                                                           cTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAAtttatttta                                                    >  2:14556/1‑131 (MQ=255)
                                                                                                            ttCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCaaa                                                                         >  3:213509/1‑109 (MQ=255)
                                                                                                            ttCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCaaa                                                                         <  4:213509/109‑1 (MQ=255)
                                                                                                                   tGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATc                                                                       <  1:830801/104‑1 (MQ=255)
                                                                                                                   tGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATc                                                                       >  2:830801/1‑104 (MQ=255)
                                                                                                                       gACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAatttatt                                                       <  4:307367/116‑1 (MQ=255)
                                                                                                                       gACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAatttatt                                                       >  3:307367/1‑116 (MQ=255)
                                                                                                                                cgcgTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAatttat                                                        <  3:908591/106‑1 (MQ=255)
                                                                                                                                cgcgTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAatttat                                                        >  4:908591/1‑106 (MQ=255)
                                                                                                                                        cccTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGa                                                                                                               >  3:737637/1‑43 (MQ=255)
                                                                                                                                        cccTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGa                                                                                                               <  4:737637/43‑1 (MQ=255)
                                                                                                                                          cTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAAc                                                                               >  2:724569/1‑73 (MQ=255)
                                                                                                                                          cTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAAc                                                                               <  1:724569/73‑1 (MQ=255)
                                                                                                                                           tGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCt  >  3:110039/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
TCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCT  >  NC_000913/2041012‑2041299

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