Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,124,5680GT34.6% 34.5 / 13.1 26T192N (ACC→AAC) gcvT‑Iaminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (8/9);  new base T (9/0);  total (17/9)
Fisher's exact test for biased strand distribution p-value = 9.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCTCATTG  >  AE015451/1124306‑1124867
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                            
gAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGAACGGTACCGAGAGTTCGTCGCCGTCTTCGCCGGTCTAGCCCGACGGGCTGAccaaaccgtccaaccccagc                                                                                                                                                                                                                                                                        >  2:742956/1‑282 (MQ=255)
             ggacgcGGGGAGCCGGCACAACGAGCACCGACGGCACGGCGCGCCTTGGAGATTGACGACACGAGGATGCCATCGACCGGCGTAGTCTCTATGACCAGCTCGTGGCCATGAAGGAACAGGCCAGCTTCCAGGCGCAGCGATGTGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCATTCCACACCCAGCAGTTTTACCcgga                                                                                                                                                                                                                                                           <  2:297050/296‑2 (MQ=255)
                     gAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGGCGGCCTTCTGCTTGTAGAACAAGAGGCTGACGAAGCAAGCCACACCCAGCCGTTTTACCCGGCGGAACTGc                                                                                                                                                                                                                                                   >  4:566352/1‑301 (MQ=255)
                                     gaaccgtagacacggctcaacgtGGAAAGTGCCGACAGCAGGCTGTCTGCGATCGGCGTAGTCTCGCTGTTCATGGCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGGGGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTGTTACCCGGCGGAACTGCATGAAGGTCATGCCgg                                                                                                                                                                                                                                   <  4:743216/279‑1 (MQ=255)
                                           tCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGCTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCCGATGACGAAGCAGTCGCCACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATggggagagggg                                                                                                                                                                                                                             >  2:77140/1‑291 (MQ=255)
                                                                         ccgctcttccgcacGATCGGCGTAGTCTCGCTGTTAATGTCGGGTCCATACAGGCACAGGCCAGCTTCCAGGCGCATCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTg                                                                                                                                                                                               <  1:614527/288‑1 (MQ=255)
                                                                                        aTCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCTTCTTCTGCTTGTAAAAAGAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTgagatcggaagagcac                                                                                                                                                                                >  3:614529/1‑286 (MQ=255)
                                                                                                 gTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGACAGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGAATGAAGGTCATGCCGGCCACTGCTGGAGCCAGGCGTGACAGCACCTTGACTGCAGCCGGGcgtatcagggcggtcc                                                                                                                                                                       >  4:526895/1‑286 (MQ=255)
                                                                                                                      ccATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCGTGAAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTCTCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGTGCAGGGCACGCTCCTCTAACatcg                                                                                                                                                  >  4:688360/1‑298 (MQ=255)
                                                                                                                                        gCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCCGGTTGAAGAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGc                                                                                                                                >  4:565009/1‑301 (MQ=255)
                                                                                                                                                      cGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCGTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCCGCTCCACCTCgatcggaagagctccctgtacg                                                                                                                  >  2:26213/1‑280 (MQ=255)
                                                                                                                                                        ctcagcgacgcttttcagatatacgggtGTACTGCCGGCTGGGACATCAGGCGGCGGGCCAAGGGTTCGGCATCGGTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGttt                                                                                                                 <  2:348263/274‑1 (MQ=255)
                                                                                                                                                                           ccggtcctggTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTg                                                                                             <  4:726847/292‑1 (MQ=255)
                                                                                                                                                                            gATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGCCTTCTGCTTGAAAAAAAAAAGGCAGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGTAACTGCATGAAGGTCATGCCGGCCACTTCTGTAGCCAGGCGTTCCAGCACCTTTACTGCAGCCGGGCCTTGCATGGCGGTCAGGGCACCTTCCTCGACCAGCGGCTTCACCTCCCAACGCCGTACGATGTGGGTTTGCATGTGAGCCAGGTCCTGc                                                                                            >  4:623167/1‑301 (MQ=255)
                                                                                                                                                                              tCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGTAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGtttagatcggaagagcacacgtctg                                                                                           >  4:348347/1‑279 (MQ=255)
                                                                                                                                                                                  cTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAACGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGAAGGTGAGCCAGGTCCTGCTCCttt                                                                                      >  3:66567/1‑300 (MQ=255)
                                                                                                                                                                                   tGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTGGCACGGTAGCCAGGTCCTgagagcggg                                                                                     >  2:726587/1‑293 (MQ=255)
                                                                                                                                                                                              cTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGAAAAACAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGAGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGGTCCAGCACCTTAACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGAGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCATGCACCTTGCAGGCGGCGTTg                                                                          >  4:779416/1‑301 (MQ=255)
                                                                                                                                                                                                     aTCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTTACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGttaaatcggg                                                                   >  3:621671/1‑293 (MQ=255)
                                                                                                                                                                                                                          cagacgtgtgctcgtccgatctagaTCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGACCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGATCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCa                                              <  4:264418/279‑1 (MQ=255)
                                                                                                                                                                                                                               ggttcagacgtgtgctcttccgatctAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTg                                         <  3:88560/276‑1 (MQ=255)
                                                                                                                                                                                                                                         ttcagacgtgtgctctgccgatctttCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATc                                <  4:368020/277‑1 (MQ=255)
                                                                                                                                                                                                                                           acgtgtgctcttccgatctgTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAg                              <  3:300714/282‑1 (MQ=255)
                                                                                                                                                                                                                                              cccctacacgacgctcttccgatctggTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGgtc                           <  2:128837/276‑1 (MQ=255)
                                                                                                                                                                                                                                                agaTCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTACTTGCAGGCGTCGTTGACGACCAGGAACAGGGTGTCCTCGccaaagatcggaagagcggcgtggag                        >  2:264352/1‑277 (MQ=255)
                                                                                                                                                                                                                                                         aGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTgagatcggaagagcgtcgtgtggggaa               >  1:88557/1‑275 (MQ=255)
                                                                                                                                                                                                                                                              gTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAgagatcggaagagcgtcgtg           >  1:300715/1‑282 (MQ=255)
                                                                                                                                                                                                                                                                       ggTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGgtcagatcggaagagcacacgtctgaac  >  4:128866/1‑276 (MQ=255)
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                            
GAAGATCGCCTCGGCACCCGGGAAGCCGGCAGCACGTGCACCGTCGGCACGGCGCACCTTGGAAATTGCCCACAGCAGGCTGGCTTCGATCGGCGTAGTCTCGCTGTTCATGTCGTGGCCATACAGGCACAGGCCAGCTTCCAGGCGCAGCGAGTCGCGTGCGCCCAAGCCGATCGGCTGTACTTCCGGCTCGGCCATCAGGCGGCGGGCCAAGGCGTCGGCAGCGTTGACCGGTACCGAGATCTCGTAGCCGTCTTCGCCGGTGTAGCCCGAGCGGCTGACGAAGCAGTCCACACCCAGCAGTTTTACCCGGCGGAACTGCATGAAGGTCATGCCGGCCACTTCTGGAGCCAGGCGTTCCAGCACCTTGACTGCAGCCGGGCCTTGCAGGGCGAGCAGGGCACGCTCCTCGAACAGCGGCTGCACCTCGCAACGGCTGCCGATGTGGGTTTGCAGGTGAGCCAGGTCCTGCTCCTTGCAGGCGGCGTTGACGACCAGGAACAGGGTGTCCTCGCCCAGGTTGGCAACCATCAGGTCGTCGAGAATGCCGCCTTGCTCATTG  >  AE015451/1124306‑1124867

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: