Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,969,7940AC43.9% 101.7 / 11.5 57G397G (GGT→GGGtapmethyl‑accepting protein IV
Reads supporting (aligned to +/- strand):  ref base A (32/0);  new base C (0/25);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.01e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTG  >  NC_000913/1969645‑1969943
                                                                                                                                                     |                                                                                                                                                     
cgGGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCAc                                                                                                                                                      >  1:709900/1‑151 (MQ=255)
 gggTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCAcc                                                                                                                                                     >  1:495612/1‑151 (MQ=255)
      aCCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGcc                                                                                                                                                 >  2:435339/1‑150 (MQ=255)
       ccGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCCCCTGCCAc                                                                                                                                               <  2:350506/151‑1 (MQ=255)
          aaCTGACAATATCGATAATGGTCGCGGGGGCGTTATTCACCAGTTCCGACCCCTGCGGGACACGTTTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCGGCCTGCGCGCTACGGCGGGAAAGATTGCGCACTTCCCCTGCCACt                                                                                                                                              <  2:888010/149‑1 (MQ=255)
           aCTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTAcc                                                                                                                                           >  2:363593/1‑151 (MQ=255)
                cAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCaa                                                                                                                                       >  2:814025/1‑150 (MQ=255)
                 aaTATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCaaa                                                                                                                                      >  1:628747/1‑150 (MQ=255)
                 aaTATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCaaa                                                                                                                                      >  2:22212/1‑150 (MQ=255)
                       gATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCACCCCTTTGACCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCa                                                                                                                                  <  2:789739/148‑1 (MQ=255)
                             ggTCGCGGCGGCGTTATTCCCCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTCCGGCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGc                                                                                                                               <  1:569166/145‑1 (MQ=255)
                                cgcgGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCGGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGctc                                                                                                                       <  1:703411/150‑1 (MQ=255)
                                  cggcggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTcc                                                                                                                    >  1:228784/1‑151 (MQ=255)
                                         gTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGcgcg                                                                                                              >  1:806153/1‑150 (MQ=255)
                                              tCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGc                                                                                                        >  1:372170/1‑151 (MQ=255)
                                                   aGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTcc                                                                                                    >  2:660682/1‑150 (MQ=255)
                                                     tttCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACt                                                                                                 >  2:370359/1‑151 (MQ=255)
                                                     tttCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACt                                                                                                 >  1:378421/1‑151 (MQ=255)
                                                      ttCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTg                                                                                                >  2:405161/1‑151 (MQ=255)
                                                            cccTGCTGGACACGTTTGACTGACTCTTCGATCAGCCCTTTGATCTTTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCCCCTGCCACTACCACAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAt                                                                                          <  1:305621/151‑1 (MQ=255)
                                                            cccTGCTGGACACGATTGACTGCCTCTTCGATCACCCCTTTGATCTTTTTTGCTGCCTGCGCGCTACGGCGGGCAAGTTTGCGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAt                                                                                          <  2:632484/151‑1 (MQ=255)
                                                             ccTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt                                                                                         >  1:666689/1‑151 (MQ=255)
                                                             ccTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt                                                                                         >  2:750823/1‑151 (MQ=255)
                                                                gctgGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGTGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCACTTCCACTGCCGCATTc                                                                                        >  2:361462/1‑149 (MQ=255)
                                                                  tgGACACGATTGCCTGCCTCTCCGATCAGCCTTTTGATCTCTTTGGCTCCCTCCGCGCTCCGGCGGGCAAGATGGCCCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                                        <  2:141731/147‑1 (MQ=255)
                                                                       acGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGcc                                                                                  >  2:485037/1‑148 (MQ=255)
                                                                           ttGCCTGACTCTCCGACCACCCTTTTGTTCTTTTTTGCTCCCGGCGCGCGCCGGCGGGCAAGTTTCCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAata                                                                            <  2:643041/150‑1 (MQ=21)
                                                                           ttGACGGCCTCTCCGTTCACCCTTTTGATCTTTTTTCCTCCCGGCGCGCTCCGGCTGGCAAGATTCCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAata                                                                            <  1:273177/150‑1 (MQ=21)
                                                                              actgactCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTgg                                                                        >  1:815207/1‑151 (MQ=255)
                                                                               ctccctctTAGAGCGGCCCTTTACTCTTTTTTGCTCCCCCCGCGCTGGGGCGGGGAAGTTTCCGCACTCCCCCTGCCACTACCGCAAAGCCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATAtt                                                                          <  2:556816/144‑1 (MQ=16)
                                                                                tgactCTTCGATCAGCCTTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTgg                                                                        <  1:132282/149‑1 (MQ=255)
                                                                                tgactCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTc                                                                      >  2:646098/1‑151 (MQ=255)
                                                                                   ctctTCCATCACCCTTTTGTTCTTTTTGGCTCCCTCGGCGCTAGGTCTGGCAAGTTTGCCCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTg                                                                         <  2:889938/145‑1 (MQ=21)
                                                                                    tctTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTgg                                                                   >  1:398817/1‑150 (MQ=255)
                                                                                          aTCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCa                                                             >  1:705271/1‑150 (MQ=255)
                                                                                               ccttttgtccttttttGCGCCCTGCCCGTTCCGGCGGGAAAGATTCCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGc                                                              <  2:173562/132‑1 (MQ=18)
                                                                                                catttGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGt                                                      >  2:861180/3‑151 (MQ=255)
                                                                                                           ttttGCTGCCTGCGCGCTAGGGCTGGAAAGTTTGCCCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGa                                                   <  2:759547/143‑1 (MQ=37)
                                                                                                           ttttGCTGCCTGCGCGCGACGGCGGGAAAGTTTGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATaa                                                <  2:169033/146‑1 (MQ=37)
                                                                                                           ttttGCCCCCTGCGCGCTACGGCTGGCAAGATTGCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGc                                             <  1:555650/149‑1 (MQ=255)
                                                                                                            tttGCGCCCGGCGCGCTAGGGCTGGCAAGATTGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTg                                           <  1:645079/150‑1 (MQ=37)
                                                                                                             ttGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGa                                          >  1:128664/1‑150 (MQ=255)
                                                                                                                ctgcctgcGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAAt                                      >  2:492784/1‑151 (MQ=255)
                                                                                                                 tgcctgcGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATg                                     >  1:848405/1‑151 (MQ=255)
                                                                                                                   cctgcGCGCTACGGCTGGCAAGATCGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATa                                        <  1:673465/146‑1 (MQ=255)
                                                                                                                       cgcgcTCCGGCGGGCAAGTTTGCGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGcc                                <  1:58020/150‑1 (MQ=255)
                                                                                                                               ggCTGGAAAGTTTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttttt                        <  1:659907/150‑1 (MQ=25)
                                                                                                                               ggCTGGAAAGATTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttttt                        <  1:405211/150‑1 (MQ=255)
                                                                                                                                gCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGc                      >  1:389148/1‑151 (MQ=255)
                                                                                                                                      agttttGCGCCCTCCCCCTGCCACTACCGCAAACCCTCGCCCCTGCTCTCCGGCGCGGGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGc                  <  1:539033/145‑1 (MQ=25)
                                                                                                                                      aaGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCtg                >  1:838761/1‑151 (MQ=255)
                                                                                                                                        gATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCtggt              >  1:709379/1‑151 (MQ=255)
                                                                                                                                         aTTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCtggtg             >  1:726457/1‑151 (MQ=255)
                                                                                                                                             cgcCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGa         <  1:95868/151‑1 (MQ=255)
                                                                                                                                               cccTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGAtc       <  1:802152/149‑1 (MQ=255)
                                                                                                                                                   tCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCt   >  1:713938/1‑151 (MQ=255)
                                                                                                                                                    cACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTg  >  1:731225/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTG  >  NC_000913/1969645‑1969943

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