Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,287,7190TA43.6% 57.8 / 39.8 55intergenic (‑348/‑45)acs/nrfAacetyl‑CoA synthetase/nitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base T (25/6);  new base A (0/24);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 2.42e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

AGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCT  >  NC_000913/4287575‑4287868
                                                                                                                                                |                                                                                                                                                      
agGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAAc                                                                                                                                                 >  2:891732/1‑151 (MQ=255)
    aGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGcc                                                                                                                                             >  2:488518/1‑151 (MQ=255)
              tacctctAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACGTTGTAGGTGCGGTTACGGTCGTAAAAAACAAAAAAAACCGCCATTg                                                                                                                                         <  2:657732/142‑1 (MQ=255)
              tAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGcaaca                                                                                                                                    >  1:547687/1‑150 (MQ=255)
               aaCTCTAAAGTGGTATTTTACATGCACTTACAATTGATAAAAGACAACATTCACAGTGTGGTTATTTGTTACACATGGGGGCGAGCAATGTAATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGcaaca                                                                                                                                    <  2:861967/149‑1 (MQ=255)
               aaCTCAAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAAt                                                                                                                                  <  2:610462/151‑1 (MQ=255)
                      aaGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACACTAAAAACCGCCATTGCAACAATGGCGCa                                                                                                                            >  1:310494/1‑150 (MQ=255)
                       aGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAAt                                                                                                                          >  2:344526/1‑151 (MQ=255)
                          ggTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCg                                                                                                                       >  1:898338/1‑151 (MQ=255)
                             aTTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGAt                                                                                                                    >  1:921975/1‑151 (MQ=255)
                                 tACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGa                                                                                                                  >  2:657261/1‑149 (MQ=255)
                                     tGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAATATCAATAAAAACCGCCATTGCAACAATGGCGCCATTCGGAt                                                                                                                    >  2:81417/1‑143 (MQ=255)
                                       cACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCt                                                                                                           >  2:405085/1‑150 (MQ=255)
                                        aCTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt                                                                                                         >  2:81124/1‑151 (MQ=255)
                                        aCTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt                                                                                                         >  1:416449/1‑151 (MQ=255)
                                           tACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTAACACATAGGGGCGAGTAATGTAAGGACAGTGTAGGTGCGGTTATTGTCGTAAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATg                                                                                                        <  1:116861/149‑1 (MQ=255)
                                              aaTTGATGAAAGCCAACATTCACAGTGTGATTATTTTTAACACATAGGGGCGAGAAAGTTCATGCCAGTGTAGGTGCGGTTCCTTTCGGGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAg                                                                                                   <  2:248982/151‑1 (MQ=255)
                                                     tAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAAAGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaa                                                                                             >  2:582170/1‑150 (MQ=255)
                                                          acaacaTTCACAGTGTGGTTATTGGTTACACATAGGGGCGAGCAATGCCATGACAGTGTAGGTGCGTTTACTGTCGTGAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATaa                                                                                        <  2:603997/150‑1 (MQ=255)
                                                            aacaTTCGCAGTGTGATTATTTGTTCCAAAAGGGGACGATCAACGTAATGACTGTGTATGACCGGTTATTGGCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacg                                                                                     <  1:217112/151‑1 (MQ=255)
                                                                 tCACAGTGTGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgcacg                                                                                 <  2:692588/150‑1 (MQ=255)
                                                                   acaGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCg                                                                              >  1:531357/1‑151 (MQ=255)
                                                                    caGTGTGTTTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGTTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCg                                                                              <  2:802694/150‑1 (MQ=255)
                                                                            ttttttGTTACAAAGAGGGGGGGGAAATTTCAAGACATTGTAGGTCCGGTTACTGTCGGGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGt                                                                             <  2:581838/140‑1 (MQ=255)
                                                                               tttGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCtt                                                                  >  1:117751/1‑151 (MQ=255)
                                                                                   ttACAAATGGGGGCGAGCAATGTAATGACGGGGTGGGTGCGGTTACTGTCGGAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAt                                                                <  1:244624/149‑1 (MQ=255)
                                                                                   ttACAAATGGGGGCGAGCAAGGTAATGACAGTGTAGGTGGGGTTCTTGTGGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAt                                                                <  1:922411/149‑1 (MQ=255)
                                                                                       acaagggggggAGAAATGTCATGACGGTGTGGGGGCGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGa                                                             <  2:358193/143‑1 (MQ=255)
                                                                                        caTGGGGGCGAGAAATGTCATGACAGTGTGGGGGCGTTTACTGTCGTGAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAtt                                                           <  1:146854/149‑1 (MQ=255)
                                                                                         atgggggCGAACAATGTCATGACAGTGTAGGTGGGTTTACTGGCGGAAAAAACAATAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAtt                                                               <  1:637252/141‑1 (MQ=255)
                                                                                         aTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCt                                                        >  2:47434/1‑151 (MQ=255)
                                                                                            ggggggAGAATTTTCATGCCAGTGTGGGTGGGGTTCTGGTGGGAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAtt                                                           <  2:75940/145‑1 (MQ=255)
                                                                                              gggaagaaaTGCCATGACGGGTTGGGGGCGTTTACTGTGGGAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttttt                                                   <  2:523561/144‑1 (MQ=255)
                                                                                                    cAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt                                             >  1:10117/1‑151 (MQ=255)
                                                                                                     aaTTTCATGACAGTGGAGGTGGTTTTCCTTTCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttc                                                  <  2:547262/145‑1 (MQ=255)
                                                                                                     aaTTTCATACCAGTGGAGGGGCGTTTACTGTGGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttt                                              <  2:414793/149‑1 (MQ=255)
                                                                                                     aaTGTCATGACGGTGTGGGGGGGTTTCCTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt                                             <  1:728837/150‑1 (MQ=255)
                                                                                                           aTGACGGTGTAGGTGGGTTTATTGTCGTAAAAAACAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttt                                               <  1:925508/142‑1 (MQ=255)
                                                                                                            tGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCt                                      >  1:612753/1‑150 (MQ=255)
                                                                                                            tGACAGGGTGGGTGCGTTTCCTGTCGTGAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCa                                           <  1:618468/145‑1 (MQ=255)
                                                                                                             gACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCCCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTCTTTTTCACTTCTGtt                                   >  1:758981/1‑151 (MQ=255)
                                                                                                              aCAGTGTAGGTGGGTTTCCTGCGGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGt                                    <  1:725623/150‑1 (MQ=255)
                                                                                                               cAGTGTAGGTGCGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTc                                  <  1:213688/151‑1 (MQ=255)
                                                                                                                 gtgtGGGTGGGTTTACTGTGGAAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGtt                                   <  1:224574/148‑1 (MQ=255)
                                                                                                                 ggtagggTGCGGTAAGTGTCGAAAAAAAAAAAAAAAACCGCCATTGCAACAAGGGCGCAATTCGGATGAAGCCCCTAGGAAAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGtt                                   <  2:41817/143‑1 (MQ=255)
                                                                                                                   gtgggTGCGGTTACTGTCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTg                                     <  2:513653/141‑1 (MQ=255)
                                                                                                                    tggggggggTTGCGGTAGTCAAAAAAAAAAAAAGCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGAAAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTg                                     <  1:224561/136‑1 (MQ=255)
                                                                                                                          ggtttactgtgtaaaaaaaaaaaaaaaaCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaa                       <  2:178750/137‑1 (MQ=255)
                                                                                                                           ggTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaa                       >  1:798708/1‑150 (MQ=255)
                                                                                                                             ttctggccgaaaaaaaaaaaaaaaaCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAAc                     <  2:517299/139‑1 (MQ=255)
                                                                                                                               actgtggaaaaaaaaaaaaaaaaCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACgg                   <  2:179021/141‑1 (MQ=255)
                                                                                                                                    cGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTc             >  1:671175/1‑151 (MQ=255)
                                                                                                                                      tGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTccc           >  1:175345/1‑151 (MQ=255)
                                                                                                                                       gAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTccc           >  2:351337/1‑150 (MQ=255)
                                                                                                                                        aaaaaaaaaaaaaaCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgct                 <  1:145400/143‑1 (MQ=255)
                                                                                                                                        aaaaaaaaTAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGc         <  2:189514/151‑1 (MQ=255)
                                                                                                                                         aaaaCAAAAAAAATCCGCCATTGCAACAATGGCGCATTTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGc         <  2:276959/151‑1 (MQ=255)
                                                                                                                                              aaTAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAAcc   >  1:531298/1‑151 (MQ=255)
                                                                                                                                               aTAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCt  >  2:351023/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
AGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCT  >  NC_000913/4287575‑4287868

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