Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,312,2130GA40.7% 74.9 / 33.4 59L12F (CTT→TTT) yciAacyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base G (9/26);  new base A (24/0);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 2.15e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTC  >  NC_000913/1312065‑1312347
                                                                                                                                                    |                                                                                                                                         
ccTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGAt                                                                                                                                         <  2:602061/151‑1 (MQ=255)
      aCCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCt                                                                                                                                   <  2:33870/151‑1 (MQ=255)
          gCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGCTGGGTAAAAAAAATAAAAAAAAAAAGCCCTGAgg                                                                                                                               >  1:686290/1‑151 (MQ=255)
          gCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGTAAAGAACGGAAAAAAAAAAAGCCCCGAgg                                                                                                                               >  2:103034/1‑151 (MQ=255)
                 cACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGAAGGGTAAAAGAAGGAAAAAAAGaaaggccctgggg                                                                                                                              >  2:599338/1‑133 (MQ=255)
                        cgACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGtgttg                                                                                                                <  2:683524/151‑1 (MQ=255)
                             gTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGCCGGGTAAAAAAAGTAAAAAAAAAAAGCCCCTAGGGGCGTTAtggttttttg                                                                                                              >  1:750483/1‑141 (MQ=255)
                             gTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAAAGGAAAAAAAAAACGCCCCGCgggggggtgtttgtt                                                                                                                  >  1:403573/1‑133 (MQ=255)
                                   aaTTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTa                                                                                                    <  2:46920/151‑1 (MQ=255)
                                      ttCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGCTAAAATAAGTAAAAAAAAAAAGCCCCCAGGGGGGGTATGTGTTTTTGACATGGGaa                                                                                                   >  1:466756/1‑149 (MQ=255)
                                             gCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGAcc                                                                                           <  1:736213/150‑2 (MQ=255)
                                               cAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGGTAGGGTAAGAAAGCAAAACAAAATAGCCCTGGGGGACGGGATTTATTTTTGACATTGTTAAAAccaagt                                                                                           >  2:292972/1‑145 (MQ=255)
                                                aGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGTAAAATAAGTAAAAAAAAAAAGCCCTGAGGGGCGGTTTGTGTTGTAGACATGGTAAAACAGACTTaa                                                                                       >  1:435307/1‑150 (MQ=255)
                                                  aaTAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGAACGGATAAAAAAACTAAAACAAAAACGCCCCGAGGGGGGTTTTGGTTTGGAGACATGGTGAAAACGACTTaa                                                                                       >  1:92459/1‑148 (MQ=255)
                                                     aGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGa                                                                                  <  1:518287/150‑1 (MQ=255)
                                                          cgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGCTAAAAAAAGTAAAAAAAAAACGCCCTGGGGGGCGGTATGTGGTGTGGAAATAGTTAAAACGACCTTAAAAGaaaa                                                                               >  1:650769/1‑148 (MQ=255)
                                                           gccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATcac                                                                           <  2:904685/151‑1 (MQ=255)
                                                             cAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAg                                                                         <  1:205394/151‑1 (MQ=255)
                                                               atatCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGgc                                                                       <  2:344262/151‑1 (MQ=255)
                                                                 atCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCGGGGCTAAAAGAAGTAAAACAAAAATCCCCTCCGGGGGGGTATTTGTTGTTGACATGGTAAAACCCACATTAAAAGGaaaaaaaacaagc                                                                      >  2:753476/1‑142 (MQ=255)
                                                                 atCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGgcg                                                                      <  2:342436/150‑1 (MQ=255)
                                                                    cATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGAt                                                                  <  1:399810/151‑1 (MQ=255)
                                                                        tgtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAATACGTAAAAAAAAAACGCCCCGAGGGGCGGTATGTGTTTTTGAAAAGGTAAAAAACACTTTAAAAGAAAAACAAAGGCGCGaaggt                                                               >  2:556343/1‑146 (MQ=255)
                                                                        tgtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCCCAAAAAAGGAAAAAAAAAACGGCCCGAGGGGGGGTATGTGTTGTAGAAATGGTAAAAACGAACTAAAAAAAAAAAAACAAGCGCGGTg                                                                 >  1:932780/1‑148 (MQ=255)
                                                                         gtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGcta                                                              <  2:934465/150‑1 (MQ=255)
                                                                           gACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGCTAAAAAAAGTCAAACAAAATCGCCCCGAGGGGCGGTATGTGTTGTAGGAAAGGTAAAAAAGACTTAAAAAACAAAACACAAGGGCGGGGctac                                                             >  2:464746/1‑149 (MQ=255)
                                                                           gACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGctacta                                                           <  2:889000/151‑1 (MQ=255)
                                                                            aCATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTGTCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTAt                                                          >  2:704081/1‑151 (MQ=255)
                                                                             cATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAAAAAAAACGCCCTGAGGGGCGGTATGTGTTGTTGACAAGGTAAAAACCACTTAAAAAGAAAATCACAGGCGCGGTGct                                                               >  2:662931/1‑145 (MQ=255)
                                                                             cATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGTAAAAAAAAAACGCCCCGGGGGGCGGTTTTTGTTTTAGAAAAGGGAAAAACGACTTAAAAAGAAAAAAAAAGGCGGGGTGcga                                                              >  2:613566/1‑144 (MQ=255)
                                                                              aTTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGTAATGCTAAAGAACAAAAAACAAGAACGGCCCGTGGGGGGTTTTGTGGTGGCAGGAACGGAAAAACCAACTAAAAAGAAAAACACAGGGGGGAAGatacta                                                           >  2:331560/1‑148 (MQ=255)
                                                                                 aaccaaccACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGACAGCCAAAGAACGTAAAAAAAAAACACCCTCTGGGGGGGTATTTGGTTTTGAAAATGTAAAAAAGAATTTAAAAGAAAACAAAAGGCGCGAttccccttttatat                                                     >  1:674335/1‑138 (MQ=255)
                                                                                  accaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATTGCGAAAGGAAGAAAAACAAGAAAGCCCTGAGGGGCGGTATGTTGTGTGGACAAAGAAAAAACCGCTTTAAAAGAGAAAAACAGGGGCGagtgcagtattagttt                                                    >  2:346680/1‑136 (MQ=255)
                                                                                    caaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCTGGGTAAAAAACCGAAAAAAAAAACGCCCCGAGGGGGGGTTTTTGTTTTTGAAAAGGGAAAACCCACCTTAAAAGAAAAACACAAGGGCGGTGGTACTATTAtttt                                                    >  2:48012/1‑149 (MQ=255)
                                                                                      accaccAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTg                                                <  1:605337/151‑1 (MQ=255)
                                                                                            aaaGATGTCACCATTGGCATTGGTATCGGCGGGGAGGGCGAAAGAACGGAAAAAAAAATTGCCCCGAGGGGCGGTTTGTGTTGTTGACATGGTAAAAACGACTTAAAAAGAAAATAACAGGGGCGATGGTCCTATCTTTTTTCTtaggag                                           >  2:408955/1‑150 (MQ=255)
                                                                                            aaaGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGTAAAAAAAAAACGCCCTGAGGGGCGGTTTGTGTTGTTGAAATGGGAAAAACCACCTAAAAAGAAAAACACAAGGGCGGGGCTACTATTATTTTTCTtaggag                                           >  2:221384/1‑150 (MQ=255)
                                                                                             aaGATGTCCCCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGaa                                         <  1:734011/151‑1 (MQ=255)
                                                                                             aaGATGTCACCATTGGCATTGGTATCGGCGGGGATGGCTAAAGAAAGGAAAAAAAAAACGGCCCGAGGGGGGGTATGTGTTGTAGACATGGGAAAAACGCCTTAAAAAGAAAAACACAAGGGCGGTGCTAATATGATTTTTCTgagg                                             >  1:225255/1‑147 (MQ=255)
                                                                                             aaGATGTCACCATTGGCATTGGTATCGGCGGGGACGGCTGAAATAAGTAAAAAAAAATAGCCCCGGGGGGAGTTATTTGTTGTTGAAATGGTAAAAACGACTTTAAAGGAAAAACACAGGCGGGAGGGTACTATTATTTTTCTgcgg                                             >  2:51129/1‑147 (MQ=255)
                                                                                              aGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAAc                                        <  2:573520/151‑1 (MQ=255)
                                                                                              aGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAAc                                        <  2:774590/151‑1 (MQ=255)
                                                                                              aGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAAc                                        <  1:622938/151‑1 (MQ=255)
                                                                                                    cACCATTGGCATTGGTATCGGCGGGCATTGCCAAAGAACCTAAAAAAAAAATGCCCCGCGGGGCGGTATTTGTTGTTGACATGGGAAAAAACACTTAAAAAGGAAAACACAAGGGGGGTGCTACTATGGTTTTTCTtaggaggaaaaagg                                   >  2:454564/1‑142 (MQ=255)
                                                                                                    cACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGTAAAACAAAAACGCCCCCAGGGAGGGTTTGTGTTGTTGAAATGGGAAAAACCACTTAAAAAGAAAAAAACAAGGGCGGTGCCACTATGATTTTTCTGAGGGGAAAAaagg                                   >  2:782382/1‑150 (MQ=255)
                                                                                                     aCCATTGGCATTGGTATCGGCGGGCCTGGCTTAAAGAAGTAAAAAAAAAACGCCCTTAGGGGGGTTATTTGTTGTAGATATGGTAAAAAAGACTTAAAAGAAAAAACAAAAGCGCGAGGCTACTATTATTTTTCTTAGGAGAAAaga                                     >  1:194850/1‑147 (MQ=255)
                                                                                                           ggcattggTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACa                            <  2:561761/150‑1 (MQ=255)
                                                                                                             cattggTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAgcg                         <  2:668756/151‑1 (MQ=255)
                                                                                                                 ggTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAg                     <  1:752453/151‑1 (MQ=255)
                                                                                                                    aTCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGccc                  <  1:51118/151‑1 (MQ=255)
                                                                                                                      cggcggGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCg                <  2:294813/151‑1 (MQ=255)
                                                                                                                       ggcggGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCg                >  1:235418/1‑150 (MQ=255)
                                                                                                                          gggCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAgg            >  1:492655/1‑151 (MQ=255)
                                                                                                                             cATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACt         >  1:161672/1‑151 (MQ=255)
                                                                                                                                ggCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTgg       <  2:900463/150‑1 (MQ=255)
                                                                                                                                 gCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCa     <  2:219417/151‑1 (MQ=255)
                                                                                                                                   tAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGGTACAGCGCCAGCCCCGAAGGACTGGCAg    >  1:306003/1‑150 (MQ=255)
                                                                                                                                   tAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGt   >  1:727662/1‑151 (MQ=255)
                                                                                                                                    aaaGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTc  <  1:414098/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                         
CCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTC  >  NC_000913/1312065‑1312347

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