Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4410AC42.1% 96.3 / 20.1 57G384G (GGT→GGGyqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (32/1);  new base C (0/24);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 3.32e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACC  >  NC_000913/2984292‑2984585
                                                                                                                                                      |                                                                                                                                                
gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac                                                                                                                                                 >  2:978742/1‑151 (MQ=255)
 aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCGATGGTCAATGCCACAcc                                                                                                                                                >  2:1086597/1‑151 (MQ=255)
 aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAcc                                                                                                                                                >  2:1232977/1‑151 (MQ=255)
 aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAcc                                                                                                                                                >  2:1129929/1‑151 (MQ=255)
  atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAcc                                                                                                                                                >  1:511604/1‑150 (MQ=255)
    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCt                                                                                                                                              >  2:846074/1‑150 (MQ=255)
              aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACCCCCTGGCCCccgccg                                                                                                                                   <  1:862474/151‑1 (MQ=255)
                    aaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACa                                                                                                                              >  2:655139/1‑150 (MQ=255)
                      aaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGCTACaa                                                                                                                             >  1:989496/1‑149 (MQ=255)
                          aGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGtt                                                                                                                       >  2:796002/1‑151 (MQ=255)
                          aGGGAGCCAAGAAAATGTATGTCGACTTGCCAGCCATAATTATTCAACAGATTAATTACGGCAGATATTTATCGGTGAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACACTAGCCCCCGCCGATACAAAGCg                                                                                                                         >  1:512065/1‑149 (MQ=255)
                           gggccccaaaaaaaTTTATGTCCACTTCCCACCCATTATTATTCAAAAGATAAATGCCGGAAAATTTTTAGCGGAAAGGAGGAAATGCAGATCTTTTCACCACGTTAAATGTTAAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGtt                                                                                                                       <  1:889891/140‑1 (MQ=255)
                           gggAGCCAAAAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATGGGTAAGGAGAAAATGCAGAGTTATTCGTCCCGTTAAATGGTCAATCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTg                                                                                                                      <  2:1093220/151‑1 (MQ=255)
                             taGCCCAGAAAATGTATGTCCACCTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGc                                                                                                                     >  1:578277/2‑150 (MQ=255)
                                   aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                              >  1:369792/1‑151 (MQ=255)
                                     aaaaTGTATGTCGATTTCCCAGCAATTATTATTAAAAAGATTAATGCGGGCAGATTTTTATGGGAAGGGAGAAAATGCAGAGTTTTTCGTCACTTTAAATGGTAAAGGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTc                                                                                                                <  2:734372/147‑1 (MQ=255)
                                       aaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                          >  2:1009642/1‑151 (MQ=255)
                                                 gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                >  1:156786/1‑151 (MQ=255)
                                                    ttGCCAGCCATCTTTATACAACAGATAAAGGACGGCAGATTTTTAGGGGTAGGGAGAAAAGGGAGAGCTTTTCGCCACGTTAAATGGTAAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTa                                                                                                  <  1:194366/146‑1 (MQ=255)
                                                     tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                            >  2:392824/1‑151 (MQ=255)
                                                      gccagcctTCATTATCCAACAGATTAATGACGGCAGATATTTTTCGGTAGGGAGGAAAGGCAGAGTTTTTCGTCACGTTAATGGGTAATGCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                            <  2:576358/150‑1 (MQ=255)
                                                       ccagcaaTCATTATAAAACAAATTAATGGCGGGAGATATTTATCGGAAGGGAGAAAATGCAGAGCTATTCGTCACGTTAAATGGTCAGTCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                            <  2:334255/149‑1 (MQ=255)
                                                       ccacccaTAATTATCCAACAGATTAAGGACGGCATATTTTTATCGGTAGGGAAGAAAGCCAGCGCTATTCGTCCCGTTCAATGGTCATTCCCCCCCCCCGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa                                                                                           <  2:25043/150‑1 (MQ=255)
                                                           ccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt                                                                                      >  2:131485/1‑151 (MQ=255)
                                                                attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                 >  2:794048/1‑151 (MQ=255)
                                                                 ttattaAACAGATAAATGAGGGAAGATTTTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGTCCCTTTAAATGGTCAACGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                    <  2:703999/147‑1 (MQ=255)
                                                                   attAAACAGATTAATGACGGAAAATTTTTAGGGGTAAGGAGAAAATGCAGAGTTTTTCGCCCCTTTAAAGGGAAATGCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa                                                                               <  2:1030670/150‑1 (MQ=255)
                                                                       aaaagttaaaggagggCAAATTTTTAGCGGAAAGGAGAAAAGGAAGAGTTTTTCTTCACTTTAAATGGTAAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                                <  1:899407/131‑1 (MQ=255)
                                                                        aCAGATAAAGGACGGAAATTATTTACGTGTAGGGAGGAAATGCAGAGCTTTTCGTCCCGTTAAATGGTCAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                  <  1:1301399/142‑1 (MQ=255)
                                                                            aTTAATGCCGGAAAATTTTTATCGGAAGGGAGGAAAGGCAGTGCTTTTCGTCACGTTAAATGGTCATTCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaaa                                                                      <  2:149064/150‑1 (MQ=255)
                                                                             ttAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                    >  1:1236285/1‑151 (MQ=255)
                                                                                 tGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                 >  1:223945/1‑150 (MQ=255)
                                                                                  gACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt                                                               >  1:239441/1‑151 (MQ=255)
                                                                                    cGGCAGTTATTTAGCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGAAAAGCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                             <  2:263278/151‑1 (MQ=255)
                                                                                       cAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                          >  2:343215/1‑151 (MQ=255)
                                                                                        atatttttATAGGAAAGGAGAAAAGCCAGATTTTTTCGTCCCTTCCATTGGTCATGCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                            <  1:1159899/143‑1 (MQ=255)
                                                                                        aGATTTTTATGGGAAAGGAGAAAATGCAGAGTTTTTCGTCACGTTCGATGGTAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg                                                           <  2:1005767/149‑1 (MQ=255)
                                                                                           tatTTAGCGAAAGGGAGGAAATGAGGACTTTTTCTCCACGTTAAATGGTAAAGGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg                                                           <  1:75787/146‑1 (MQ=255)
                                                                                              ttATCGGTAAGGAGAAAATGCAGAGCTTTTCGCCACTTTAAATGGTAAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                    <  2:875736/150‑1 (MQ=255)
                                                                                               tATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                   >  2:1052659/1‑150 (MQ=255)
                                                                                                      aaGGAGGAAATGGATATTTTTTCGTCACGTTCAATGTTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTa                                           <  1:1214760/151‑1 (MQ=255)
                                                                                                      aaGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                            >  2:1147736/1‑150 (MQ=255)
                                                                                                        ggaggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCTGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATc                                         <  1:639200/151‑1 (MQ=255)
                                                                                                         gagaaaggcaagtgCTTTTCGTCACGTTAAAGGGTCAAGCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg                                        <  1:585544/138‑1 (MQ=255)
                                                                                                          aggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                       >  2:569433/1‑151 (MQ=255)
                                                                                                                  cAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                >  2:1688/1‑150 (MQ=255)
                                                                                                                   agggCTTTTCGTCCCTTTAAATGGTCAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGa                               <  2:476770/147‑1 (MQ=255)
                                                                                                                     agCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                            >  2:1110838/1‑151 (MQ=255)
                                                                                                                        tATTCTTCCCTTCCAATGGTAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                            <  1:795972/148‑1 (MQ=255)
                                                                                                                         aTTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGa                         >  2:202883/1‑150 (MQ=255)
                                                                                                                            cGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                     >  2:618805/1‑151 (MQ=255)
                                                                                                                               cACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc                  >  2:597617/1‑151 (MQ=255)
                                                                                                                               cACGTTCAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                     <  2:1033228/148‑1 (MQ=255)
                                                                                                                                aCGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc                  >  1:743130/1‑150 (MQ=255)
                                                                                                                                  gTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc               >  1:216441/1‑151 (MQ=255)
                                                                                                                                        tGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgac         >  2:405982/1‑151 (MQ=255)
                                                                                                                                        tGGTCAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattga          <  2:592331/150‑1 (MQ=255)
                                                                                                                                               tGCCACACCCTGTCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAcc  >  1:363170/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACC  >  NC_000913/2984292‑2984585

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: