Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 450,764 C→T A283T (GCT→ACT)  cyoA ← cytochrome o ubiquinol oxidase subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913450,7640CT100.0% 221.2 / NA 69A283T (GCT→ACT) cyoAcytochrome o ubiquinol oxidase subunit II
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/34);  total (35/34)

AATTTTCCGAACATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGCCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGTTCGGCGACTGCTTCGC  >  NC_000913/450628‑450899
                                                                                                                                        |                                                                                                                                       
aaTTTTCCGAACATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGTCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCa                                                                                                                                       <  2:678529/139‑1 (MQ=255)
    ttCCGAACATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATaaa                                                                                                                                   >  1:1443976/1‑139 (MQ=255)
     tCCGAACATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACTGTGAGTCATa                                                                                                                                     >  2:685396/1‑136 (MQ=255)
     tCCGAACATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATa                                                                                                                                     <  1:685396/136‑1 (MQ=255)
     tCCGAACATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                                  <  2:135263/139‑1 (MQ=255)
     tCCGAACATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                                  <  2:95093/139‑1 (MQ=255)
                        ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACtt                                                                                                                                <  2:292959/122‑1 (MQ=255)
                        ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACtt                                                                                                                                >  1:292959/1‑122 (MQ=255)
                        ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                               >  2:11433/1‑139 (MQ=255)
                        ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                               >  2:1241804/1‑139 (MQ=255)
                        ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                               >  2:385056/1‑139 (MQ=255)
                        ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                               >  2:923691/1‑139 (MQ=255)
                          tCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAAc                                                                                                             <  1:1241804/139‑1 (MQ=255)
                                        gggCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCAc                                                                                               >  1:936530/1‑139 (MQ=255)
                                        gggCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCAc                                                                                               >  2:193590/1‑139 (MQ=255)
                                         ggCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACg                                                                                              <  2:886953/139‑1 (MQ=255)
                                           cGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGtt                                                                                            >  1:1210812/1‑139 (MQ=255)
                                            ggATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTg                                                                                           <  1:1102777/139‑1 (MQ=255)
                                                tCCGCGTGGCTCATGTCCATGCCTTCCATATCGTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGgaga                                                                                       <  2:1052349/139‑1 (MQ=255)
                                                    cgTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTa                                                                                                                            >  1:326105/1‑98 (MQ=255)
                                                    cgTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTa                                                                                                                            <  2:326105/98‑1 (MQ=255)
                                                    cgTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAata                                                                                   <  1:385056/139‑1 (MQ=255)
                                                      tGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATAtt                                                                                 <  1:11433/139‑1 (MQ=255)
                                                        gCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATAtt                                                                                 >  1:1239731/1‑137 (MQ=255)
                                                        gCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATAtt                                                                                 <  2:1239731/137‑1 (MQ=255)
                                                         cTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCa                                                                              <  2:936530/139‑1 (MQ=255)
                                                          tCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCAc                                                                             <  2:105945/139‑1 (MQ=255)
                                                           cATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCAcc                                                                            <  1:109964/139‑1 (MQ=255)
                                                             tGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATAtt                                                                                 >  1:644636/1‑132 (MQ=255)
                                                             tGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATAtt                                                                                 <  2:644636/132‑1 (MQ=255)
                                                                   tGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCAc                                                                             >  2:575871/1‑130 (MQ=255)
                                                                   tGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCAc                                                                             <  1:575871/130‑1 (MQ=255)
                                                                       ttCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCg                                                                >  2:699420/1‑139 (MQ=255)
                                                                        tCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCa                                                                                                >  2:1116237/1‑106 (MQ=255)
                                                                        tCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCa                                                                                                <  1:1116237/106‑1 (MQ=255)
                                                                                ttCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGcgcg                                                       >  2:1382371/1‑139 (MQ=255)
                                                                                  cGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACAt                                                                                                                     >  2:770674/1‑75 (MQ=255)
                                                                                  cGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACAt                                                                                                                     <  1:770674/75‑1 (MQ=255)
                                                                                   gtgtGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTg                                                                                           >  1:171765/1‑100 (MQ=255)
                                                                                   gtgtGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTg                                                                                           <  2:171765/100‑1 (MQ=255)
                                                                                    tgtgCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                                  <  2:1108664/60‑1 (MQ=255)
                                                                                    tgtgCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                                  >  1:1108664/1‑60 (MQ=255)
                                                                                    tgtgCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCa                                                   >  1:845796/1‑139 (MQ=255)
                                                                                        cgcTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACtt                                                                                                                                >  2:610635/1‑58 (MQ=255)
                                                                                        cgcTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACtt                                                                                                                                <  1:610635/58‑1 (MQ=255)
                                                                                        cgcTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCg                                                                                                                   <  2:227628/71‑1 (MQ=255)
                                                                                        cgcTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCg                                                                                                                   >  1:227628/1‑71 (MQ=255)
                                                                                              gCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACtt                                                                                                                                <  2:1073045/52‑1 (MQ=255)
                                                                                              gCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACtt                                                                                                                                >  1:1073045/1‑52 (MQ=255)
                                                                                                   ccTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAg                                    <  1:902000/139‑1 (MQ=255)
                                                                                                    cTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGc                                   <  2:1393215/139‑1 (MQ=255)
                                                                                                    cTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGc                                   <  2:845796/139‑1 (MQ=255)
                                                                                                       ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTg                                                                                           <  2:913371/80‑1 (MQ=255)
                                                                                                       ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTg                                                                                           >  1:913371/1‑80 (MQ=255)
                                                                                                       ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGcc                                  >  1:1139195/1‑137 (MQ=255)
                                                                                                       ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGcc                                  <  2:1139195/137‑1 (MQ=255)
                                                                                                         ggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGt                              >  2:973894/1‑139 (MQ=255)
                                                                                                           ctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCa                            >  2:323853/1‑139 (MQ=255)
                                                                                                            tggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAg                           >  2:595494/1‑139 (MQ=255)
                                                                                                              ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                               >  2:826495/1‑53 (MQ=255)
                                                                                                              ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                               <  1:826495/53‑1 (MQ=255)
                                                                                                              ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGaa                                         >  1:408944/1‑123 (MQ=255)
                                                                                                              ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGaa                                         <  2:408944/123‑1 (MQ=255)
                                                                                                              ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGAc                         >  2:185475/1‑139 (MQ=255)
                                                                                                               gTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACa                        <  1:828572/139‑1 (MQ=255)
                                                                                                                               tACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGTTCGGCGACTg        >  2:1452809/1‑139 (MQ=255)
                                                                                                                                aCCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGTTCGGCGACTGc       >  1:442618/1‑139 (MQ=255)
                                                                                                                                aCCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGTTCGGCGACTGc       >  2:379728/1‑139 (MQ=255)
                                                                                                                                     gAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGTTCGGCGACTGCTTCGc  <  1:415735/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
AATTTTCCGAACATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGCCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGTTCGGCGACTGCTTCGC  >  NC_000913/450628‑450899

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