Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,718,463 T→G 34.2% Q65P (CAG→CCG)  CJ019_01614 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,718,4630TG34.2% 43.4 / 11.1 60Q65P (CAG→CCG) CJ019_01614hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (7/30);  new base G (14/8);  total (22/38)
Fisher's exact test for biased strand distribution p-value = 7.87e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGG  >  ACB122/1718335‑1718592
                                                                                                                                |                                                                                                                                 
tgatgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGgcgcg                                                                                                                             >  3:281360/1‑135 (MQ=255)
tgatgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGgcgc                                                                                                                              <  7:200792/134‑1 (MQ=255)
    gaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAg                                                                                                                         <  4:305230/135‑1 (MQ=255)
          cAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCAc                                                                                                                   >  7:92066/1‑135 (MQ=255)
          cAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCggggccc                                                                                                                   >  5:283245/1‑131 (MQ=255)
          cAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGTcgc                                                                                                                   >  1:187350/1‑133 (MQ=255)
              ttGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCa                                                                                                               <  8:317357/135‑1 (MQ=255)
                  cATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGGCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAgggg                                                                                                           <  3:299517/135‑1 (MQ=255)
                       aGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGCGGGCGCGCCCGGGGGGCGGCAGGGGCGTcc                                                                                                      >  1:12390/1‑135 (MQ=255)
                             ctGCGCGCACTACCCGCATTGCGGTGGCCAGGCTTGCCCTGCGCGGGGTCATAGGGGAGTTGGAACGGCAGCGGCACGTTGCCTTGGACCGGGCCGGCCTGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGc                                                                                                <  2:65935/134‑1 (MQ=255)
                             aTGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGGGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGc                                                                                                <  7:296059/135‑1 (MQ=255)
                             aTGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGc                                                                                                <  3:373363/135‑1 (MQ=255)
                             aTGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGAGGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGc                                                                                                <  7:202106/135‑1 (MQ=255)
                                cgcgcACTTACCGCATAGCGGTGGGCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGac                                                                                             <  7:24075/135‑1 (MQ=255)
                                cgcgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGGGAGAGGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGac                                                                                             <  7:79312/135‑1 (MQ=255)
                                cgcgcACTTACCGCATAGCGGTGGCCAGGCGTGACCGGCGCGGGGGCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGac                                                                                             <  5:86/135‑1 (MQ=255)
                                cgcgaAATTACCGCATAGGGGGGGCCCGGCTTGACCGGGGGGGGGTCATCGGGGCGAGGCACCGGCCGCGGCACGTGGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGac                                                                                             <  4:23378/135‑1 (MQ=255)
                                 gcgcACTTACCTCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTTGCGCGCCAGGGTCACCGCCGGGGCGTGCAGCcgcgcc                                                                                             >  6:267052/1‑132 (MQ=255)
                                 gcgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGGGCGCCCGGGTCCCGGCAGGGGCGGCCAACAGCgccc                                                                                            >  6:22373/1‑132 (MQ=255)
                                 gcgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGTCACGGCCGGGGCGGCCCTCAGCGcca                                                                                            >  4:153257/1‑135 (MQ=255)
                                 gcgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGGCCCGGCCGGGGCGGCCCTCCGCGaca                                                                                            >  5:66975/1‑135 (MQ=255)
                                  cgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGacac                                                                                           <  8:57311/135‑1 (MQ=255)
                                     aCTTACCGATTAGCGGTGGCCAGGCTGGACGTGCGCGGGGTCATTGGTGACAGGGAACGGCAGCGGCAGGTTGCGTTTGACGGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATc                                                                                        <  4:327797/135‑1 (MQ=255)
                                      cTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCGGGGGACGGCCGGGGCGGGCCGCAGCGACCCCTcc                                                                                       >  8:16813/1‑135 (MQ=255)
                                           cGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCGGGTCACGGCCGGGGCGGCCCTCCGCGACCCCGCCTGCAg                                                                                  >  8:202106/1‑135 (MQ=255)
                                                  cGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCAc                                                                            <  4:332367/134‑1 (MQ=255)
                                                  cGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGGCCCGGCCGGGGCGGCCCTCCCCGCCACATCCGTCCGTTccccc                                                                           >  2:346922/1‑132 (MQ=255)
                                                  cGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGGGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACg                                                                           <  2:281646/135‑1 (MQ=255)
                                                          aGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGGCCCGGCGGGGGCGGCCCTCAGCGCCCCATCCTGCGGGGCCGCGCGCAGGGt                                                                   >  3:325256/1‑135 (MQ=255)
                                                              ttGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGCCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcg                                                               <  7:130837/135‑1 (MQ=255)
                                                              gtgcccTGCGGGGGGGGATGGGTGAGCTGCACCGGACGCGCCACGTGGCGGTTGCCCGGGCGGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGGCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcg                                                               <  7:86637/131‑1 (MQ=255)
                                                              gtGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcg                                                               <  4:59614/134‑1 (MQ=255)
                                                               tGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGTCACGGCCGGGGCGTCCCTCAGCGCCACCTCCTGCAGTTCCCCGCGCCGGGGggcgg                                                              >  6:77998/1‑135 (MQ=255)
                                                               tGACCTGCGCGGGGTCATAGGTGAGATGCAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcgg                                                              <  5:267052/135‑1 (MQ=255)
                                                                   gtGCGCGGGGTCTTAGGTGACATGTAACGCCAGCGGCAGGTTGCCTTTGACCGGGCCGTCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcagg                                                          <  3:183853/134‑1 (MQ=255)
                                                                   cTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcagg                                                          <  8:265502/135‑1 (MQ=255)
                                                                   cTGCGCGGGGGCATAGGTGAGATGGACCGGCCGCGGCACGTTGCCTGTGACCGGGCCGGCCGGGCGCGCCCGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcagg                                                          <  1:44959/135‑1 (MQ=255)
                                                                    tGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGCGGGCGCGCCCGGGGCGCGGCCGGGGCGGCCCTCAGCGCCCCCGCCTGCAGGTCCCCCCTCCGGGGggcgggcggc                                                         >  5:187422/1‑130 (MQ=255)
                                                                         ggggTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCAGGCGCGCCACGGGCACCGCCGGGGCGTCCATCCGCGACACCTCCAGCAGTTACACGCTCACGGGggcggccggcggccg                                                    >  3:107889/1‑133 (MQ=255)
                                                                             tCATAGGTTAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCa                                                <  3:129386/135‑1 (MQ=255)
                                                                              cATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGGGGGCGCGCCCGGGGCCCGGCCCGGGGGGGGCGCAGCGACACATCCTGCCGGTACACGCGGAGGGTGGCGGCCGGCGGCCGGGcca                                               >  8:244836/1‑133 (MQ=255)
                                                                               aTAGGTGAGAGGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAAt                                              <  3:184402/135‑1 (MQ=255)
                                                                                   gggagaGGGAACGGCAGCGCCACGTTGCCTTTGCCCGGGCCGGCCTGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgcgc                                          <  2:187350/133‑1 (MQ=255)
                                                                                   gTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCGGGGGACCGCAGGGGCGGCCCTCAGCGCCACAACCAGCCCTTCCCCGCGCAGGGTGGCGGCCCGCGGCAGGGCCATgcgc                                          >  3:262796/1‑135 (MQ=255)
                                                                                   gTGAGAGGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgcgc                                          <  7:244836/135‑1 (MQ=255)
                                                                                         tGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGTCACGGCAGGGGCGGCCATCAGCGCCACAACCTGCACGTCCCCGCGCATGGTGGCGGCCGGCGGCCGGGCCATGCGCGGCAgg                                    >  6:249962/1‑135 (MQ=255)
                                                                                                gCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGGGCGGCCGGGGCCCGGCCGGGGGGGCCATCCGCGCCACAGCCGGCAGATCCCCGCGCCGGGGGGCGGGCGGCGGCCGGGCCATGCGCGGCAgggggcag                              >  8:59111/1‑128 (MQ=255)
                                                                                                  cggggCAGGTTGCCTTTGCCGGGGCGGGCGGGGGGCGCCTGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACt                           <  2:180218/132‑1 (MQ=255)
                                                                                                  cggggCACGTGGCCTTTGACGGGGCCGGCGTGGTGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGTCAATGCGCTGCAGGTAGAAGACt                           <  1:158365/132‑1 (MQ=255)
                                                                                                  cgCGGCAGGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACt                           <  5:344722/134‑1 (MQ=255)
                                                                                                  aGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACt                           <  5:189187/135‑1 (MQ=255)
                                                                                                  aGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACt                           <  2:31027/135‑1 (MQ=255)
                                                                                                  aGCGGCACGTTGCCGTTGACCGGGCCGGCCTGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACt                           <  4:223204/135‑1 (MQ=255)
                                                                                                    cGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCGGGTCACGGCAGGGGCGTCCATCAGCGCCCCAACCAGCAGTTCCACGCGCAGGGTGGCGGCCGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTc                         >  6:186783/1‑135 (MQ=255)
                                                                                                                tttGCCCGGGCCGGCTGGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCt             <  5:372116/135‑1 (MQ=255)
                                                                                                                   gACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTg            <  3:243589/133‑1 (MQ=255)
                                                                                                                    aCCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGcc         <  7:19258/135‑1 (MQ=255)
                                                                                                                     ccGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGGAGAAGACTTCACCGTCCAGGCTGGcc         <  7:19263/134‑1 (MQ=255)
                                                                                                                         gccggccgGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGtt    <  5:34358/135‑1 (MQ=255)
                                                                                                                         gccggccTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGtt    <  8:92066/135‑1 (MQ=255)
                                                                                                                           cggccGGGCGCGCCAGGGTCACGGCCGGGGCGTCCATAAGCGACCCAACCTGGAGTTCCACGCGCCGGGTGGCGGCCGGCGGCCGGGCAATGCGCGGCCGGGTGAACACTTCCACGTCCAGGGGGGCCTGGGtgg  >  6:207347/1‑135 (MQ=255)
                                                                                                                                |                                                                                                                                 
TGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGG  >  ACB122/1718335‑1718592

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