Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,028,861 G→A A51V (GCG→GTG)  ygfS ← putative 4Fe‑4S ferredoxin‑type oxidoreductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,028,8610GA98.2% 220.1 / NA 57A51V (GCG→GTG) ygfSputative 4Fe‑4S ferredoxin‑type oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/34);  minor base T (0/1);  total (22/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 = 1.00e+00

ATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCGCACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGAT  >  NC_000913/3028714‑3029007
                                                                                                                                                   |                                                                                                                                                  
aTTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCacac                                                                                                                                                 >  2:70678/1‑151 (MQ=255)
     aaaCGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATg                                                                                                                                            >  2:337874/1‑151 (MQ=255)
          ggCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTc                                                                                                                                       <  2:630190/151‑1 (MQ=255)
             acgcacTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAg                                                                                                                                    <  2:552026/151‑1 (MQ=255)
               gcacTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTc                                                                                                                                  >  2:172639/1‑151 (MQ=255)
                cacTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTc                                                                                                                                  <  2:276681/150‑1 (MQ=255)
                    gACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGt                                                                                                                             >  2:47256/1‑151 (MQ=255)
                       gCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCGAGTCGCTGTAc                                                                                                                           >  2:88083/1‑150 (MQ=255)
                       gCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTAcc                                                                                                                          >  2:123035/1‑151 (MQ=255)
                        catctctGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCt                                                                                                                         <  1:228750/148‑1 (MQ=255)
                         aaCTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCtt                                                                                                                        <  2:356222/151‑1 (MQ=255)
                          aCTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTCCCTTc                                                                                                                       >  2:730613/1‑151 (MQ=255)
                               gACAGCCCATACGACGGGCAGAATTGGTCTGCTCCACCTGCTCTACCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCg                                                                                                                  <  2:352572/151‑1 (MQ=255)
                               gACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGGCGCTGTACCTTCAGc                                                                                                                    <  1:298640/149‑1 (MQ=255)
                                 cAGCCAATACAACGGGCAGAATTGGTCTGCGCCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCTGCTCACTGATGCTGTCCAGTCGCTGTACCTTCAGCCgg                                                                                                                 <  2:417411/150‑1 (MQ=255)
                                 cAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCggg                                                                                                                >  2:500211/1‑151 (MQ=255)
                                 cAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCggg                                                                                                                >  2:662201/1‑151 (MQ=255)
                                 cAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCgg                                                                                                                 <  2:26480/150‑1 (MQ=255)
                                  aGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCgggg                                                                                                               >  1:649183/1‑151 (MQ=255)
                                     cAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTaa                                                                                                            <  2:471413/151‑1 (MQ=255)
                                     cAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTaa                                                                                                            <  1:571602/151‑1 (MQ=255)
                                      aaTACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTaaa                                                                                                           <  1:91966/151‑1 (MQ=255)
                                      aaTACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTaaa                                                                                                           <  2:544036/151‑1 (MQ=255)
                                       aTACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTaaaa                                                                                                          <  2:176865/151‑1 (MQ=255)
                                             cGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACAt                                                                                                    <  1:557178/151‑1 (MQ=255)
                                              gggCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATc                                                                                                   <  1:449165/151‑1 (MQ=255)
                                                 cAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGc                                                                                                <  1:252740/151‑1 (MQ=255)
                                                    aaTTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCAtt                                                                                             <  2:42850/151‑1 (MQ=255)
                                                    aaTTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCAtt                                                                                             <  2:607416/151‑1 (MQ=255)
                                                       tGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCa                                                                                           >  1:62503/1‑150 (MQ=255)
                                                         gTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACt                                                                                        >  2:676394/1‑151 (MQ=255)
                                                             gCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTg                                                                                    >  1:729399/1‑151 (MQ=255)
                                                                     tGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGa                                                                              >  1:765739/1‑149 (MQ=255)
                                                                      gCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAgg                                                                           <  1:210775/151‑1 (MQ=255)
                                                                      gCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAg                                                                            <  2:765388/150‑1 (MQ=255)
                                                                         cGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTg                                                                        <  1:472209/151‑1 (MQ=255)
                                                                          gCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTg                                                                        >  1:339104/1‑150 (MQ=255)
                                                                           cccATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAg                                                                      <  2:291856/151‑1 (MQ=255)
                                                                              aTCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCAAACTGATGGCACATCACCGGCACAATGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCg                                                                    >  1:738007/1‑150 (MQ=255)
                                                                                   cAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACcaca                                                              <  1:728892/151‑1 (MQ=255)
                                                                                    aGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCTCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACcacac                                                             <  1:346390/151‑1 (MQ=255)
                                                                                          cccACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTa                                                       <  2:329068/151‑1 (MQ=255)
                                                                                            cACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACt                                                     >  2:419262/1‑151 (MQ=255)
                                                                                                gggCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTcac                                                 >  1:489740/1‑151 (MQ=255)
                                                                                                    aaGCGCCAACACACGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGt                                             >  2:536038/1‑151 (MQ=255)
                                                                                                         ccAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGc                                        >  1:260735/1‑151 (MQ=255)
                                                                                                            acacAAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGc                                     <  2:688694/151‑1 (MQ=255)
                                                                                                             cacaAGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGcc                                    >  1:162346/1‑151 (MQ=255)
                                                                                                                aaGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAAt                                 <  1:699692/151‑1 (MQ=255)
                                                                                                                 aGGGGCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATa                                <  2:686703/151‑1 (MQ=255)
                                                                                                                    ggCGTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAg                             <  1:83883/151‑1 (MQ=255)
                                                                                                                       gTTTTCACACTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCa                          <  2:270527/151‑1 (MQ=255)
                                                                                                                               acTGATGGCACATCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGAtt                  >  2:372876/1‑151 (MQ=255)
                                                                                                                                       cacaTCACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATaa          <  1:272473/151‑1 (MQ=255)
                                                                                                                                            cACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGcc                                    <  2:215857/120‑1 (MQ=255)
                                                                                                                                            cACCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGcc                                    >  1:215865/1‑120 (MQ=255)
                                                                                                                                             aCCGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACg    <  2:724528/151‑1 (MQ=255)
                                                                                                                                               cGGCACACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGAt  <  1:495588/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCATCGTCAGCGCCCCCACGGGGCAAGCGCCAACACAAGGGGCGTTTTCACACTGATGGCACATCACCGGCGCACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGAT  >  NC_000913/3028714‑3029007

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