Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 477,430 T→C E396E (GAA→GAG pdeB ← c‑di‑GMP phosphodiesterase PdeB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913477,4300TC98.3% 208.6 / NA 60E396E (GAA→GAGpdeBc‑di‑GMP phosphodiesterase PdeB
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (23/35);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGTCGACATCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATGCTGCTGTGGATG  >  NC_000913/477295‑477561
                                                                                                                                       |                                                                                                                                   
aaTGTCGACATCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAg                                                                                                                                 >  2:52173/1‑140 (MQ=255)
  tGTCGACATCCAGATCCCGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTg                                                                                                                               >  1:138190/1‑140 (MQ=255)
   gTCGACATCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGa                                                                                                                              >  2:145519/1‑140 (MQ=255)
       aCATCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAAATATCGGGGAGCTAGTTGTCGGATCGGCAAAATCGCGCTCAGTGAGTTc                                                                                                                          >  2:192998/1‑140 (MQ=255)
          tCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAg                                                                                                                       <  1:541493/140‑1 (MQ=255)
          tCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTACCGGGGCGCAGGTCTCCGGATCGGCATACTCCTGCTCAGTGAGTTCAAg                                                                                                                       >  2:353014/1‑140 (MQ=255)
            cAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAgcg                                                                                                                     <  2:293522/140‑1 (MQ=255)
              gATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGcgcg                                                                                                                   >  2:132786/1‑140 (MQ=255)
               aTCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGa                                                                                                                  >  1:101902/1‑140 (MQ=255)
                tCCTGTAAATAACTTGAACTTGCATGCCCAGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGAt                                                                                                                 <  1:340455/140‑1 (MQ=255)
                tCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGAt                                                                                                                 <  2:680058/140‑1 (MQ=255)
                tCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGAt                                                                                                                 <  1:397433/140‑1 (MQ=255)
                tCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGAt                                                                                                                 <  1:428461/140‑1 (MQ=255)
                tCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGAt                                                                                                                 <  1:52173/140‑1 (MQ=255)
                tCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGAt                                                                                                                 <  1:523972/140‑1 (MQ=255)
                tCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGAt                                                                                                                 <  1:158040/140‑1 (MQ=255)
                tCCCGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGAt                                                                                                                 <  2:138190/140‑1 (MQ=255)
                   tGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGAtctg                                                                                                              >  1:588459/1‑140 (MQ=255)
                     tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGAtctgtc                                                                                                            >  1:426142/1‑140 (MQ=255)
                     tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGAt                                                                                                                 <  2:533839/135‑1 (MQ=255)
                     tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGAt                                                                                                                 >  1:533839/1‑135 (MQ=255)
                     tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTAGCGGTCTGTGAGTGCAAGCGCGAtctgtc                                                                                                            >  1:211990/1‑140 (MQ=255)
                           aCTTAAACTTGCATAGCCCGTAGCACAATGATCAAGACAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                                      <  1:345626/140‑1 (MQ=255)
                               aaaCTTGACTACCCCGTCCCAATATCATCAGTAGACATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAAc                                                                                                  <  2:355883/140‑1 (MQ=255)
                               aaaCTTGAATACCCCGTCCCAAAATCATCAATATACATTTCATGGCCCACCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAAc                                                                                                  <  2:195085/140‑1 (MQ=255)
                               aaaCTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAAc                                                                                                  <  2:41327/140‑1 (MQ=255)
                                 aCTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCt                                                                                                <  2:577406/140‑1 (MQ=255)
                                    tGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGAt                                                                                             <  2:144687/140‑1 (MQ=255)
                                           cccGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGAt                                                                                             <  1:689491/133‑1 (MQ=255)
                                           cccGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGAt                                                                                             >  2:689491/1‑133 (MQ=255)
                                           cccGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtgatt                                                                                      >  2:128487/1‑140 (MQ=255)
                                           cccGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtgatt                                                                                      >  1:185983/1‑140 (MQ=255)
                                            ccGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGTGAGGTCAAGCGGGATCTGTCAGGGAGTAACCTGATAGtgattg                                                                                     >  2:296806/1‑140 (MQ=255)
                                             cGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCGAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtg                                                                                         <  1:33766/135‑1 (MQ=255)
                                             cGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCGAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtg                                                                                         >  2:33766/1‑135 (MQ=255)
                                                aCCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCa                                                                                 >  2:180280/1‑140 (MQ=255)
                                                     aaTCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCa                                                                            <  1:145519/140‑1 (MQ=255)
                                                                        aTGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGAttt                                                         >  1:52789/1‑140 (MQ=255)
                                                                         tGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGAtttt                                                        >  2:648499/1‑140 (MQ=255)
                                                                            cccGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTc                                                     >  1:514089/1‑140 (MQ=255)
                                                                             ccGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTcc                                                    >  2:619353/1‑140 (MQ=255)
                                                                             ccGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTcc                                                    <  2:170676/140‑1 (MQ=255)
                                                                              cGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGCTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCg                                                   >  2:550623/1‑140 (MQ=255)
                                                                                ccTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAg                                                 <  1:619353/140‑1 (MQ=255)
                                                                                ccTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAg                                                 <  2:426142/140‑1 (MQ=255)
                                                                                   cccGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTg                                              <  2:514089/140‑1 (MQ=255)
                                                                                     cGGTAGCGAGAAATTCTCGGGGCGCTGGTCTTCGGATTGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAg                                            <  1:192998/140‑1 (MQ=255)
                                                                                         aGCGAGAAATTATCGGGGCGCTAGTTTTCGGTTCGCCAATCTCGCGCTAAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACg                                        <  2:293092/140‑1 (MQ=255)
                                                                                          gCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACgg                                       <  1:132786/140‑1 (MQ=255)
                                                                                          gCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACgg                                       >  1:709541/1‑140 (MQ=255)
                                                                                          gCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACgg                                       >  1:151577/1‑140 (MQ=255)
                                                                                           cGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACggg                                      <  1:180280/140‑1 (MQ=255)
                                                                                             agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACggg                                      <  2:624832/138‑1 (MQ=255)
                                                                                             agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACggg                                      >  1:624832/1‑138 (MQ=255)
                                                                                                             cTGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGa                    >  1:582863/1‑140 (MQ=255)
                                                                                                              tGGTTTTCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGaa                   <  2:582863/140‑1 (MQ=255)
                                                                                                                   ttCGGATCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATAtg              <  1:116854/140‑1 (MQ=255)
                                                                                                                        aTCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATgctgct         <  2:185983/140‑1 (MQ=255)
                                                                                                                         tCGGCAAACTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATGctgctg        <  1:550623/140‑1 (MQ=255)
                                                                                                                               aaCTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATGCTGCTGTGGATg  <  2:151577/140‑1 (MQ=255)
                                                                                                                               aaCTCGCGCTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATGCTGCTGTGGATg  <  2:709541/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                   
AATGTCGACATCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATGCTGCTGTGGATG  >  NC_000913/477295‑477561

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