Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,512,090 G→A 6.9% A206V (GCG→GTG)  mntH ← manganese/divalent cation transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,512,0900GA 6.9% 180.1 / 2.9 59A206V (GCG→GTG) mntHmanganese/divalent cation transporter
Reads supporting (aligned to +/- strand):  ref base G (29/25);  new base A (2/2);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

TCATCGCCAGATTGACAAAACCGGCAATAGTCATGGCGATAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCAACAAACAACAGTAACCC  >  NC_000913/2511947‑2512232
                                                                                                                                               |                                                                                                                                              
tcatcGCCAGATTGACAAAACCGGCAATAGTCATGGCGATAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGcccc                                                                                                                                            >  4:288552/1‑148 (MQ=255)
  atcGCCAGATTGACAAAACCGGCAATAGTCATGGCGATAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTa                                                                                                                                          <  5:226654/148‑1 (MQ=255)
     gCCAGATTGACAAAACCGGCAATAGTCATGGCGATAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAaca                                                                                                                                       <  6:345115/148‑1 (MQ=255)
       cAGATTGACAAAACCGGCAATAGTCATGGCGATAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCATATAAATCACATGCGGCATAATCGTCGCCCCTACCACGc                                                                                                                                    >  2:111406/1‑149 (MQ=255)
         gATTGACAAAACCGGCAATAGTCATGGCGATAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCct                                                                                                                                  >  4:111628/1‑149 (MQ=255)
         gATTGACAAAACCGGCAATAGTCATGGCGATAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCct                                                                                                                                  <  7:300080/149‑1 (MQ=255)
             gACAAAACCGGCAATAGTCATGGCGATAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCctgct                                                                                                                               <  3:109544/148‑1 (MQ=255)
             gACAAAACCGGCAATAGTCATGGCGATAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCctgct                                                                                                                               <  5:309707/148‑1 (MQ=255)
               cAAAACCGGCAATAGTCATGGCGATAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCct                                                                                                                                  >  8:282885/1‑143 (MQ=255)
               cAAAACCGGCAATAGTCATGGCGATAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCct                                                                                                                                  <  7:282885/143‑1 (MQ=255)
                     cGGCAATAGTCATGGCGATAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAgg                                                                                                                         >  1:44789/1‑146 (MQ=255)
                      ggCAATAGTCATGGCGATAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGa                                                                                                                     >  7:358485/1‑149 (MQ=255)
                       gCAATAGTCATGGCGATAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGAc                                                                                                                    <  3:200611/149‑1 (MQ=255)
                                 tGGCGATAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGaa                                                                                                          <  3:111628/149‑1 (MQ=255)
                                       tAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGAAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGt                                                                                                    <  2:81035/149‑1 (MQ=255)
                                            acaTCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACt                                                                                               <  8:358485/149‑1 (MQ=255)
                                                   aTTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCACCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGAt                                                                                         >  4:158540/1‑148 (MQ=255)
                                                   aTTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCACCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGAt                                                                                         >  4:152942/1‑148 (MQ=255)
                                                         tggCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCAtt                                                                                  <  2:69550/149‑1 (MQ=255)
                                                                    cGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCCTAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGCCCGCTTCCGAAGTAGGGAAACTCGGGATCCCCATTCCTTTACCCAg                                                                       >  4:252852/1‑149 (MQ=255)
                                                                    cGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGa                                                                                          >  1:331457/1‑130 (MQ=255)
                                                                    cGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGa                                                                                          <  2:331457/130‑1 (MQ=255)
                                                                         ttgACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTgcg                                                                  >  7:346492/1‑149 (MQ=255)
                                                                           gACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCCGGATCACCATTCCTTTACCCACCTgcgc                                                                 <  4:222145/148‑1 (MQ=255)
                                                                              gcgAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAgg                                                             >  8:22165/1‑149 (MQ=255)
                                                                                    cGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGt                                                            >  1:19080/1‑144 (MQ=255)
                                                                                    cGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGt                                                            <  2:19080/144‑1 (MQ=255)
                                                                                    cGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGc                                                       <  7:22165/149‑1 (MQ=255)
                                                                                       cATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTgg                                                    <  8:363661/149‑1 (MQ=255)
                                                                                          gTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTgg                                                    >  4:366861/1‑146 (MQ=255)
                                                                                            aaaTGCTGAGTGAGCGAGGAGAGCAAATAACTCAAATGCGGCATAATCCTCCCCCCTAACACTCCTGCTGCCAGGAATACCCCTTCCGAACTAGGAAAACTCGGCATCACAAATCCTTTACCCAGCAGCGCCAGCTTAGGCTGGcaga                                                >  8:30630/1‑148 (MQ=255)
                                                                                             aaTGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGaaa                                              <  3:255132/149‑1 (MQ=255)
                                                                                             aaTGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGaa                                               >  7:322399/1‑148 (MQ=255)
                                                                                             aaTGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGaa                                               >  2:326944/1‑148 (MQ=255)
                                                                                              aTGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGaaaa                                             <  2:44789/149‑1 (MQ=255)
                                                                                              aTGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGaaaa                                             >  5:150954/1‑149 (MQ=255)
                                                                                               tGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAAt                                            <  5:341513/149‑1 (MQ=255)
                                                                                                    agtgagCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCaa                                         <  3:288552/147‑1 (MQ=255)
                                                                                                         gCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGAc                                   >  5:143977/1‑148 (MQ=255)
                                                                                                          cGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACaa                                 >  2:318192/1‑149 (MQ=255)
                                                                                                           gaggagTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCCCCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACAAt                                >  3:255479/1‑149 (MQ=255)
                                                                                                                   cAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACAATGTAAGCcg                        >  4:238249/1‑149 (MQ=255)
                                                                                                                     aataaatCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGTGAAAATCAACTCGACAATGTAAGCcgcg                      >  5:138932/1‑149 (MQ=255)
                                                                                                                     aataaatCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACAATGTAAGCcgcg                      >  2:3006/1‑149 (MQ=255)
                                                                                                                          atCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACAATGTACGCCGCGGcac                  >  8:246727/1‑147 (MQ=255)
                                                                                                                           tCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTgg                                                    >  3:23025/1‑113 (MQ=255)
                                                                                                                           tCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTgg                                                    <  4:23025/113‑1 (MQ=255)
                                                                                                                           tCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCAGAAGTAGGTTAACTCGGGATCACCATTCCTTTACCCCGCTGCGCCAGGTTAGGCTTGGAGAAACTCAACGCGACAATGTATGCCGCGTcaac                 >  4:197246/1‑148 (MQ=255)
                                                                                                                              caTGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCaacaaac             <  3:366861/149‑1 (MQ=255)
                                                                                                                              caTGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCaacaaa              >  4:149212/1‑148 (MQ=255)
                                                                                                                               aTGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCaacaaaca            >  8:193000/1‑149 (MQ=255)
                                                                                                                                tGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGATAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCAacaaacaa           >  4:254152/1‑149 (MQ=255)
                                                                                                                                tGCGGCATAATCGTCACCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCAacaaacaa           <  3:152942/149‑1 (MQ=255)
                                                                                                                                tGCGGCATAATCGTCACCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCAacaaacaa           <  3:158540/149‑1 (MQ=255)
                                                                                                                                  cGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCAacaaacaa           >  4:161788/1‑147 (MQ=255)
                                                                                                                                  cGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCAACAAacaaca         <  3:238249/149‑1 (MQ=255)
                                                                                                                                   ggCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCAACAAACAACAg        <  1:111406/149‑1 (MQ=255)
                                                                                                                                        aaTCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCGGGGAGAAAATCGACTCGACTATGTAAGCCGCGGCAACAAACAACAGTAccc   >  8:83368/1‑149 (MQ=255)
                                                                                                                                         aTCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGATAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCAACAAACAACAGTAAccc  <  3:254152/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
TCATCGCCAGATTGACAAAACCGGCAATAGTCATGGCGATAGCCACATCCCATTTGGTGGCGGAATAACGTTGTTGACGCGAACCGCCATGTAAATGCTGAGTGAGCGAGGAGTGCAAATAAATCACATGCGGCATAATCGTCGCCCCTAACACGCCTGCTGCCAGGAAGACCGCTTCCGAAGTAGGTAAACTCGGGATCACCATTCCTTTACCCAGCTGCGCCAGGTTAGGCTGGGAGAAAATCAACTCGACAATGTAAGCCGCGGCAACAAACAACAGTAACCC  >  NC_000913/2511947‑2512232

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