Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913847,4440GT44.6% 72.0 / 34.4 56I187I (ATC→ATAglnHglutamine transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (5/26);  new base T (25/0);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 2.73e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGTAATCAACGGAGCCAGTACCG  >  NC_000913/847299‑847591
                                                                                                                                                 |                                                                                                                                                   
tGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTtgatga                                                                                                                                                 <  2:336063/150‑1 (MQ=255)
 gTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGaa                                                                                                                                                <  1:950274/150‑1 (MQ=255)
     gTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGTAc                                                                                                                                            >  2:6431/1‑150 (MQ=255)
           ttCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGGACAGaa                                                                                                                                        >  2:973567/1‑148 (MQ=255)
            tctcGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGGACAGAAtagt                                                                                                                                    >  1:356765/1‑148 (MQ=255)
            tctcGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGGAGTACAGAATGtt                                                                                                                                    >  1:541080/1‑151 (MQ=255)
               cgcgCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGGAAAAAAAGttt                                                                                                                                   >  1:923468/1‑149 (MQ=255)
                 cgcAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGGTTTTTTGAAGTAAAGAAttgtt                                                                                                                                   >  1:41781/1‑143 (MQ=255)
                    aGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGGTTTTTTGAAGTACAGAAAGttt                                                                                                                                   >  2:164380/1‑144 (MQ=255)
                     gggTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATc                                                                                                                           <  1:716570/151‑1 (MQ=255)
                        ttttCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTGTTAAGTACAGAATGttttggggattgt                                                                                                                         >  2:460956/1‑140 (MQ=255)
                        ttttCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTGTGAAGTACATAAAGTTTGGGGGATAGTg                                                                                                                        >  2:1061750/1‑151 (MQ=255)
                         tttCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTg                                                                                                                        <  2:950320/150‑1 (MQ=255)
                            cAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAg                                                                                                                     <  1:883329/150‑1 (MQ=255)
                              aCGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGAAGGAAAGAATGTTTTGCGGATCGTGCAGaa                                                                                                                   >  2:256032/1‑150 (MQ=255)
                               cgcgCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACg                                                                                                                 <  1:613791/151‑1 (MQ=255)
                                  gcCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCg                                                                                                              <  2:790726/151‑1 (MQ=255)
                                     gTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGGAAAGAAAGTTTTGGGGATCGTGCAGaa                                                                                                                   >  1:994348/1‑143 (MQ=255)
                                     gTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTGTGAAGGAAAAAATGTTTTGCCTATTCTGCAggaacg                                                                                                                >  1:81008/1‑141 (MQ=255)
                                     gTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCt                                                                                                           <  1:950942/151‑1 (MQ=255)
                                       tGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGTACAGAATGTTTTGGGGATCGTGCAGAAAGGCGTCTg                                                                                                          >  2:1078241/1‑150 (MQ=255)
                                                  cGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGc                                                                                              <  1:975038/151‑1 (MQ=255)
                                                      gCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTTAAGTACAGAATGTTTTGCGGATCGTGGAGAAACGGGCCTGCGCGGGTgggggcc                                                                                            >  1:497180/1‑144 (MQ=255)
                                                       cTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGTACAGCATGTTTGGCGGGTCTTGCCGAACGGGGTTTGCGCGGTTGGTGccccct                                                                                         >  1:10006/1‑148 (MQ=255)
                                                       cTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGATTTGATTAAGTATAGAAAAGTTTGCGTAACGTGCCAGACGGCGGCTGCGCGGgtggtgggcag                                                                                          >  1:744929/1‑145 (MQ=255)
                                                          gtcgCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGTAAAGAATTTTTGGCGGATCGTGCAGAAACGCGGCTGCGCGGGTGGGGCCCCGtttc                                                                                      >  1:245389/1‑149 (MQ=255)
                                                          gtcgCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGAAGGACAGAATGTTTTGCGGATCGTGCAGAACGGCGGCTGCGCGGGTGGTGCCCAGGTcc                                                                                      >  2:991380/1‑151 (MQ=255)
                                                               tACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGGAAAGAAAGTTTTGGGTATCGTGCAGAAAGGCGGCTTCGCGGGTGGGGGCCAGtttcct                                                                                    >  2:511976/1‑144 (MQ=255)
                                                                 ccTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAgg                                                                               <  2:798635/151‑1 (MQ=255)
                                                                   tttCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTTAAGTAAAGAATGTTTGGGGTATCGTGCAGAACGGGCGCTGCGCGGGTGGGGcccccgtccc                                                                                     >  1:209724/1‑136 (MQ=255)
                                                                    ttCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGc                                                                              <  1:719048/149‑1 (MQ=255)
                                                                            cGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATg                                                                    <  2:509854/151‑1 (MQ=255)
                                                                                     gTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGAAGAAAAAAATGTTTTGCGGATCGTGCAGAAAGGGGTCTGCGCGGGTGGGGCCCCGGTCCCTATAGGCGTTATCGAAGTTCgggga                                                            >  2:1069103/1‑148 (MQ=255)
                                                                                              cgcTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGTTGAAGTAAAGAATGTTTTGCGGATCGTGCAGAACGGCGTCTGCGCGGGTGGTGCCCCGTTCCATTTAGGCGGTATCGATGTTTGGGGAATGACg                                                      >  1:420195/1‑147 (MQ=255)
                                                                                              cgcTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCa                                                    <  2:116199/149‑1 (MQ=255)
                                                                                               gcTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGCTGAAGGAAAGAAAAATTTGCGGATTGTGGAGAACGGCGGCTTCGCGGTTGGGGGCCAGTTCCACAAAAGCGTTTTTGTAGTTTGGGGACTg                                                         >  2:907259/1‑143 (MQ=255)
                                                                                                cTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTATGAAGCACAGAAAGTTTTGCGGGTCGGGCCGAAAGGCGGCTGCGCGGGTGGTGGCCACTTTCCTATATGCGGTATCCATGTTCGGGGACTg                                                         >  1:853157/1‑142 (MQ=255)
                                                                                                   ccAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCtt                                              <  1:475145/150‑1 (MQ=255)
                                                                                                     agagagTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAg                                           <  2:617098/151‑1 (MQ=255)
                                                                                                      gagagTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAGCGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAg                                           <  1:7382/150‑1 (MQ=255)
                                                                                                        gagTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGttt                                        <  2:368459/151‑1 (MQ=255)
                                                                                                         agTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGGACAAAAATTTTTGCGTATCGGGCAGAACGGCGTCTGCGCGGGTGGGGCCCCGGTCCCTATAAGGGGTATTGATGGTCGGGAAATGACGCCAATCTTTTGtttt                                       >  2:747605/1‑151 (MQ=255)
                                                                                                            cACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGAt                                    <  2:171636/151‑1 (MQ=255)
                                                                                                            cACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGAt                                    <  2:90260/151‑1 (MQ=255)
                                                                                                             aCCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGTAGTACAAGATGTTTTGCGGATCGTGCAGAACGGCGTCCGCGGGGGTGGGGCCCCGTTCCCTATAAGGGGTTTCGATGTTCGGGGACTGACGCAGATCTTTTGttttttt                                    >  1:666300/1‑147 (MQ=255)
                                                                                                                tACCGCTTTGAACTGACCGTTACCGGCGGTTTTTGTGAAGGAAAGAAAGTTTTGCGGATCGTGCAGAACGGGGGCTGCGCGGGTGGGGCCCCGTTCCCTATAGGCGGTTTCGATGTTTGGGGAATGACGCAGGTCttttgtttttttgttt                                >  2:893707/1‑151 (MQ=255)
                                                                                                                 aCCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTg                               <  2:565283/151‑1 (MQ=255)
                                                                                                                   cGCTTTGAACTGACCGTTACCGGCGGTTTTTTTGAAGTACAAAATGTTTTGCGGATCCTGGAGGACGGCGTCTGCGCGGGTGGTGCCCAGTTCCATATAGGGGTTATCGATGTTCGGGGACTGGCGCAGAACTTTTTTTTTGATGTTTGCt                             >  2:852400/1‑151 (MQ=255)
                                                                                                                    gcgttGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCtt                            <  2:492193/148‑1 (MQ=255)
                                                                                                                             tGACCGTTACCGGCGGTTTTTTTGAAGGACAGAAAGGTTTGCGTGGTGGGCAGGACGGCGTCTGCTCGGGTGGGGACCAACTTCCTATTGGCGGTTTCGGTGTTTCGGGACTGACCCAGAAACTTTGTTTTTTTGTTTGCTTTCGCGTAAt                   >  1:133660/1‑151 (MQ=255)
                                                                                                                                    tACCGGCGGTTTTGATGAAGTACATAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGTAATCAACgg             <  1:423698/150‑1 (MQ=255)
                                                                                                                                     aCCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGTAATCAACGGa            <  2:527439/150‑1 (MQ=255)
                                                                                                                                     aCCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGTAATCAACGGAg           <  1:255975/151‑1 (MQ=255)
                                                                                                                                     aCCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGTAATCAACGGAg           <  1:13707/151‑1 (MQ=255)
                                                                                                                                        ggcggTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGTAATCAACGGAGCCa        <  2:824084/151‑1 (MQ=255)
                                                                                                                                          cggTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATATGCGTTATCGATGTTCGGGATCTGACGCAGATCGTTAGTTTTGATGTTTGCTTTCGCGTAATCAACGGAGCCAg       <  1:846387/150‑1 (MQ=255)
                                                                                                                                              tttGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGTAATCAACGGAGCCAGTACCg  >  2:1052599/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGCGTCTGCGCGGTTGGTGCCCAGTTCCATATAGGCGTTATCGATGTTCGGGAACTGACGCAGATCTTTAGTTTTGATGTTTGCTTTCGCGTAATCAACGGAGCCAGTACCG  >  NC_000913/847299‑847591

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