Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,059,867 A→C 9.1% S323A (TCC→GCC)  PP_4458 ← opine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,059,8670AC 9.1% 158.9 / 6.8 55S323A (TCC→GCC) PP_4458opine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (29/21);  new base C (2/3);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

AATTGCGGGGCCGCGTTGGAAACAACCGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACGGAATTCGGCAGGC  >  NC_002947/5059719‑5060009
                                                                                                                                                    |                                                                                                                                                 
aaTTGCGGGGCCGCGTTGGAAACAACCGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGa                                                                                                                                                   <  6:35823/149‑1 (MQ=255)
             cgTTGGAACCAACCGATTTGAGTTTTAACCCGCTGGGATAACCCGCGCTGGCCAACAAGGCGTTCGCCTGCGCAGGATCGTACGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGCGGATGC‑‑‑TCCATCa                                                                                                                                   <  4:121681/149‑1 (MQ=255)
              gTTGGAAACAACCGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAg                                                                                                                                  >  2:169517/1‑149 (MQ=255)
                 ggAAACAACCGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTa                                                                                                                               >  2:122123/1‑149 (MQ=255)
                  gAAACAACCGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTa                                                                                                                               >  1:61299/1‑148 (MQ=255)
                  gAAACAACCGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTa                                                                                                                               >  7:4602/1‑148 (MQ=255)
                        acCGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCttt                                                                                                                        <  4:148715/149‑1 (MQ=255)
                         ccGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGCGGATGC‑‑‑TCCATCatac                                                                                                                                >  4:263210/1‑137 (MQ=255)
                         ccGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGCGGATGC‑‑‑TCCATCatac                                                                                                                                <  3:263210/140‑4 (MQ=255)
                            aTTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTAccc                                                                                                                    <  4:99540/149‑1 (MQ=255)
                                 aGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTa                                                                                                               <  1:255470/149‑1 (MQ=255)
                                        gcccGTTGGGATAAACCGCGCTGGACAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATAGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCACCGAACCGAAt                                                                                                         >  2:306772/2‑148 (MQ=255)
                                        aCCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAAt                                                                                                         >  7:18440/1‑148 (MQ=255)
                                         cccGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAg                                                                                                              <  8:117119/142‑1 (MQ=255)
                                         cccGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAg                                                                                                              >  7:117119/1‑142 (MQ=255)
                                           cGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGAtc                                                                                                     <  8:4602/149‑1 (MQ=255)
                                                 gATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGACCGGCACTGCAAGAAAGCCCCTGGTAGCCAGTGGGTATGATGGa                                                                                                                                                   >  4:248701/1‑100 (MQ=255)
                                                     aCCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATg                                                                                            <  4:226937/148‑1 (MQ=255)
                                                      cccGCGCTGGCCAAGAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCATTCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGtt                                                                                          <  5:288624/149‑1 (MQ=255)
                                                      cccGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGCGGATGC‑‑‑TCCATCatac                                                                                                                                >  6:266436/1‑108 (MQ=255)
                                                      cccGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGCGGATGC‑‑‑TCCATCatac                                                                                                                                <  5:266436/111‑4 (MQ=255)
                                                       ccGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGtt                                                                                          >  3:290778/1‑148 (MQ=255)
                                                         gcgcTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGAc                                                                                       <  7:10208/149‑1 (MQ=255)
                                                                          tttCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCtt                                                                      >  7:229886/1‑149 (MQ=255)
                                                                               cgTGCGCAGGATCGTTAGCGTATTGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCCCTTTCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGc                                                                    <  1:122123/149‑1 (MQ=255)
                                                                                    gcAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGAGCTCATTGATGGTGACCGACGCAGCGATAGCGTGCCGAACCtt                                                            >  4:305916/1‑149 (MQ=255)
                                                                                       tgATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGAc                                                         <  1:367298/148‑1 (MQ=255)
                                                                                       ggATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGAc                                                         <  2:61299/149‑1 (MQ=255)
                                                                                       ggATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGAc                                                         <  1:361565/149‑1 (MQ=255)
                                                                                        gATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGATACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGAcg                                                        >  5:207816/1‑149 (MQ=255)
                                                                                          tCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGAcgtc                                                      >  1:197160/1‑149 (MQ=255)
                                                                                             tAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCATAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACgtcgtc                                                   >  4:383046/1‑149 (MQ=255)
                                                                                                    aTGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGCCAGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTACCGGc                                            >  2:248243/1‑149 (MQ=255)
                                                                                                       gCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCtt                                          <  2:197160/148‑1 (MQ=255)
                                                                                                        cccGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACGTCTTTACCCGCGTAGCGAATGATCTCATTGATGTTTACCGACGCAGCGATAGCTTGc                                                                    <  2:200801/121‑1 (MQ=255)
                                                                                                        cccGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACGTCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGc                                                                    >  1:200801/1‑121 (MQ=255)
                                                                                                         ccGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGa                                        >  7:255338/1‑148 (MQ=255)
                                                                                                         ccGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGa                                        >  5:77345/1‑148 (MQ=255)
                                                                                                         ccGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGa                                        >  2:282927/1‑148 (MQ=255)
                                                                                                               cTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCg                                 <  1:185828/149‑1 (MQ=255)
                                                                                                                 gCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGt                               <  8:229886/149‑1 (MQ=255)
                                                                                                                      cAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTaag                           >  6:356054/1‑148 (MQ=255)
                                                                                                                          ccccTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTaagaaaga                      <  8:18440/149‑1 (MQ=255)
                                                                                                                            ccTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTaagaaaga                      >  2:45507/1‑147 (MQ=255)
                                                                                                                            ccTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTaagaaaga                      >  7:16083/1‑147 (MQ=255)
                                                                                                                            ccTGGTAGCCATTGGGTATGATGGAGAATCC‑‑‑GCCATCAGCTACATCTTTACCCTCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAagaga                    <  3:333568/149‑1 (MQ=255)
                                                                                                                             cTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAagagag                   >  6:11206/1‑149 (MQ=255)
                                                                                                                             cTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACAGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTaagaaaga                      >  5:33844/1‑146 (MQ=255)
                                                                                                                                 tAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACg               >  5:123722/1‑149 (MQ=255)
                                                                                                                                  aGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGATTGATATCATTGATGTTGACCGACGCAGAGATAGCTTTCCGAACCTTGACGTCGTATAACGTCTTGATATTCAGGTTAAGAAATAGATTAct               >  5:352134/1‑147 (MQ=255)
                                                                                                                                          ggTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGGGTACGGAATTCgg       >  7:91462/1‑148 (MQ=255)
                                                                                                                                          ggTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACGGAATTCgg       >  3:2924/1‑148 (MQ=255)
                                                                                                                                           gTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACGGAATTCGGCa     <  5:11206/149‑1 (MQ=255)
                                                                                                                                            tATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACGGAATTCGGCAg    >  2:176569/1‑149 (MQ=255)
                                                                                                                                            tATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACGGAATTCGGCAg    <  7:265067/149‑1 (MQ=255)
                                                                                                                                            tATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACGGAATTCGGCAg    >  5:198080/1‑149 (MQ=255)
                                                                                                                                              tgatgGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACGGAATTCGGCAGGc  <  8:16083/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AATTGCGGGGCCGCGTTGGAAACAACCGATTTGAGTTTTAACCCGTTGGGATAACCCGCGCTGGCCAACAAGGCTTTCGCGTGCGCAGGATCGTAAGCGTATGGCCCGGCACTGCAATCAAGCCCCTGGTAGCCATTGGGTATGATGGAGCATCC‑‑‑GCCATCAGCTACATCTTTACCCGCGTAGCGAATGATCTCATTGATGTTGACCGACGCAGCGATAGCTTGCCGAACCTTGACGTCGTCTAACGGCTTGATATTCCGGTTAAGAAAGAGAGTACGGAATTCGGCAGGC  >  NC_002947/5059719‑5060009

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