Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,723,459 T→C T91T (ACT→ACC U712_19195 → putative transporter ywcJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,723,4590TC96.9% 230.2 / ‑3.8 65T91T (ACT→ACCU712_19195putative transporter ywcJ
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (27/35);  total (28/37)
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

AACAGGCAGTTATTTCTCTATGGCGGATTCTCCGTTTGCTTTTCCGGCAGCCGCTGTCACATTCGGGGCCGCTATTCTGCTGATTGCTTACGGAGGCGGAGATTTATTTACTGGCAACACTTTTTATTTCACTTATACCGCGCTCCGAAAAAAAATTCGCTGGCGAGACACCCTGTACCTGTGGATGTCAAGCTATGCCGGCAATTTGATCGGCGCCATTCTATTTGCCATTCTGATCAGCGCGACAGGGCTTTTTGAG  >  CP006881/3723327‑3723585
                                                                                                                                    |                                                                                                                              
aaCAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCTGCAGCCGATGTCACTTTCGGTGCCGCTATTATGATTATTGCTTACGGAGGCGGAGATTTTTTTACCGGCAACACCTTTTATTTCACCTATAcc                                                                                                                          >  8:376342/1‑139 (MQ=255)
            tttCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCgga                                                                                                                <  3:204015/137‑3 (MQ=255)
             ttCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGGGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGgaaa                                                                                                              >  5:435801/1‑138 (MQ=255)
                   aTGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTGTTGTGATGGTTGCTTTAGGCGGGGGGGTTTTATTTTACGGCACCACCGTTTTTTTAACCGATACCACCCTCCGgcaaaaaaaaa                                                                                                       >  6:321135/1‑136 (MQ=255)
                       ctgATTCTCCGATTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTg                                                                                                   <  3:158018/137‑1 (MQ=255)
                          aTTCTCCGTTTGCCTTTCCGGCATCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCg                                                                                                >  4:438710/1‑139 (MQ=255)
                          aTTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCg                                                                                                >  1:263775/1‑139 (MQ=255)
                          aTTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCg                                                                                                >  6:313081/1‑139 (MQ=255)
                          aTTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCg                                                                                                >  3:153992/1‑139 (MQ=255)
                          aTTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCg                                                                                                >  7:64528/1‑139 (MQ=255)
                          aTTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCAGGCg                                                                                                >  8:344208/1‑139 (MQ=255)
                             ctcCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTACTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCgcg                                                                                              <  7:256991/138‑3 (MQ=255)
                               ccGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGaca                                                                                           <  6:157342/139‑1 (MQ=255)
                                cGTTTGCCTTTCCGGCTGCCGCTGTCACTTTGGGGGCAGCTATTCTGATGATTGATTACGGTGGCGGAGATTTATTTACCGGCTACACCTTTTATTTCACCTATACCGCGCTCGTGAAAAAAATCAGCTGGCggga                                                                                             <  2:181420/136‑1 (MQ=255)
                                cGTTTGCCTTTCCGGCAGCCGCTGTCAGTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGaca                                                                                           <  6:427746/138‑1 (MQ=255)
                                cGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCgcga                                                                                             <  8:328226/136‑1 (MQ=255)
                                cGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCgcga                                                                                             <  3:201328/136‑1 (MQ=255)
                                cGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCgcga                                                                                             <  2:12872/136‑1 (MQ=255)
                                cGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCgcga                                                                                             <  3:417452/136‑1 (MQ=255)
                                cGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCgcga                                                                                             <  8:37439/136‑1 (MQ=255)
                                cGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCgcga                                                                                             <  4:9567/136‑1 (MQ=255)
                                cGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGaca                                                                                           <  7:205484/138‑1 (MQ=255)
                                cGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGaca                                                                                           <  5:313081/138‑1 (MQ=255)
                                cGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGaca                                                                                           <  3:86584/138‑1 (MQ=255)
                                cGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCGCCGGAAAAAAATCAGCTGGCgcga                                                                                             <  2:16690/136‑1 (MQ=255)
                                 gTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACAcc                                                                                         <  8:325521/139‑1 (MQ=255)
                                 gTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACAcc                                                                                         <  5:321135/139‑1 (MQ=255)
                                 gTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACAcc                                                                                         <  3:390333/139‑1 (MQ=255)
                                      cctttCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCtata                                                                                    <  6:435801/137‑1 (MQ=255)
                                       ctttCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCgcga                                                                                             <  8:468520/128‑1 (MQ=255)
                                       ctttCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCgcga                                                                                             >  7:468520/2‑129 (MQ=255)
                                                               cGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCgg                                                            <  4:447937/138‑1 (MQ=255)
                                                               cGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCgg                                                            <  3:369859/138‑1 (MQ=255)
                                                                      gCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGaca                                                                                           >  8:431551/1‑100 (MQ=255)
                                                                      gCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGaca                                                                                           <  7:431551/100‑1 (MQ=255)
                                                                          ttCTGATGATTGCTTACGGAGGCGGATATTTATTTACCGGCAACACCTTTTATTTCACATATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCgg                                                <  8:301753/139‑1 (MQ=255)
                                                                          ttCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCgg                                                <  8:343386/139‑1 (MQ=255)
                                                                           tCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTATATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCgg                                                <  2:112389/138‑1 (MQ=255)
                                                                              gatgATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGgcgc                                             >  1:232002/3‑138 (MQ=255)
                                                                               atgATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACCCCTTTTATTTAACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACCCCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCa                                           >  6:375870/2‑139 (MQ=255)
                                                                               atgATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCgcga                                                                                             >  4:327766/2‑89 (MQ=255)
                                                                               atgATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCgcga                                                                                             <  3:327766/88‑1 (MQ=255)
                                                                               atgATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCa                                           >  4:435247/2‑139 (MQ=255)
                                                                               atgATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATATAATCCGCGcc                                            >  5:136227/2‑138 (MQ=255)
                                                                               atgATTGCTTACGGAGGCGGAGATTTATTTACCGGAAACACCGTTTATTTAACCGATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACCCCCGATACTTGTGGATGGCAAGCTATTCCGGCAATTTAATCCGCGCCa                                           >  6:100709/2‑139 (MQ=255)
                                                                                 gATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCAtt                                         >  3:432583/1‑139 (MQ=255)
                                                                                   ttGCTTACGGAGGCGTAGATTTATTTTCCGGCAACACATTTTATTTCACCTATACCGCGAGCCGGAAATAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCATTCt                                       <  4:18548/139‑1 (MQ=255)
                                                                                    tGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACCCCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCATTCtg                                      >  8:213511/1‑138 (MQ=255)
                                                                                    tGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCgg                                                            >  7:104407/1‑117 (MQ=255)
                                                                                    tGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCgg                                                            <  8:104407/117‑1 (MQ=255)
                                                                                    tGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCATTCtg                                      >  3:93081/1‑138 (MQ=255)
                                                                                    tGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCATTCtg                                      >  3:403364/1‑138 (MQ=255)
                                                                                    tGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTGCACCTATACCGCGCTGCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCGGGCAATTTAATCGCCGCAATACtg                                      >  6:163689/1‑138 (MQ=255)
                                                                                         aCGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCATTCTGTTTg                                  <  7:66151/138‑1 (MQ=255)
                                                                                         aCGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCATTCTGTTTg                                  <  6:238048/138‑1 (MQ=255)
                                                                                            gAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACTTATACCGCGCTCCGGAAAAAAATCAGCTGGCgcga                                                                                             >  2:66046/1‑76 (MQ=255)
                                                                                            gAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACTTATACCGCGCTCCGGAAAAAAATCAGCTGGCgcga                                                                                             <  1:66046/76‑1 (MQ=255)
                                                                                             aGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCATTCTGTTTGCCAtc                             <  7:225436/139‑2 (MQ=255)
                                                                                               gCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCATTCTGTTTGCCATCCt                           >  4:438317/1‑139 (MQ=255)
                                                                                                cGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCATTCTGTTTGCCATCCt                           >  6:47623/1‑138 (MQ=255)
                                                                                                     atttatttACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCATTCTGTTTGCCATCCTGATCAg                     >  3:256691/1‑139 (MQ=255)
                                                                                                              ccggCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCATTCTGTTTGCCATCCTGATCAGCGCGACggg            <  3:438317/137‑1 (MQ=255)
                                                                                                               cggCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCATTCTGTTTGCCATCCTGATCAGCGCGACggga           >  4:249910/2‑138 (MQ=255)
                                                                                                                        cttttATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCATTCTGTTTGCCATCCTGATCAGCGCGACGGGACTTTTTgag  <  2:263775/138‑1 (MQ=255)
                                                                                                                        cttttATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGACACCCTATACTTGTGGATGTCAAGCTATGCCGGCAATTTAATCGGCGCCATTCTGTTTGCCATCCTGATCAGCGCGACGGGACTTTTTgag  <  8:64528/138‑1 (MQ=255)
                                                                                                                                    |                                                                                                                              
AACAGGCAGTTATTTCTCTATGGCGGATTCTCCGTTTGCTTTTCCGGCAGCCGCTGTCACATTCGGGGCCGCTATTCTGCTGATTGCTTACGGAGGCGGAGATTTATTTACTGGCAACACTTTTTATTTCACTTATACCGCGCTCCGAAAAAAAATTCGCTGGCGAGACACCCTGTACCTGTGGATGTCAAGCTATGCCGGCAATTTGATCGGCGCCATTCTATTTGCCATTCTGATCAGCGCGACAGGGCTTTTTGAG  >  CP006881/3723327‑3723585

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: