Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,626,456 A→C 7.7% E37A (GAG→GCG)  PP_2298 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,626,4560AC 7.7% 154.1 / 3.1 52E37A (GAG→GCG) PP_2298hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/27);  new base C (2/2);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

GTGGCACCCCAATAAAAGCTGTTAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCGGAGCAGCTGCAGGCCATG  >  NC_002947/2626311‑2626598
                                                                                                                                                 |                                                                                                                                              
gtggCACCCCAATAAAAGCTGTTAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCgagg                                                                                                                                              <  2:190480/148‑1 (MQ=255)
gtggCACCCCAATAAAAGCTGTTAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCgag                                                                                                                                               <  5:373153/147‑1 (MQ=255)
gtggCACCCCAATAAAAGCTGTTAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCgag                                                                                                                                               <  8:210482/147‑1 (MQ=255)
       cccAATAAAAGCTGTTAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCgcggatctcctc                                                                                                                                      <  3:346860/149‑8 (MQ=255)
              aaaGCTGTTAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGCGGATCTCCTCGAGcaa                                                                                                                                <  5:140430/148‑3 (MQ=255)
              aaaGCTGTTAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTgg                                                                                                                               <  2:52807/149‑1 (MQ=255)
                  cTGTTAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATgc                                                                                                                           <  3:328964/149‑1 (MQ=255)
                      tAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGTGAGCTGGATGCGCGt                                                                                                                       >  2:266890/1‑149 (MQ=255)
                      tAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGt                                                                                                                       >  6:163562/1‑149 (MQ=255)
                       aGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCGCGAGGAGATCCGCGAGCTGGATGcgcggg                                                                                                                      >  6:127938/1‑147 (MQ=255)
                         aTTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAg                                                                                                                    <  4:311538/149‑1 (MQ=255)
                         aTTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAg                                                                                                                    >  2:40659/1‑149 (MQ=255)
                           ttCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAg                                                                                                                    >  2:346126/1‑147 (MQ=255)
                                   gTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAg                                                                                                           >  1:104135/1‑148 (MQ=255)
                                   gTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCAGCGAGCTGGATGCGCGTGAGCAGGCCCAg                                                                                                           >  7:148570/1‑148 (MQ=255)
                                     tGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCa                                                                                                         >  1:325897/1‑148 (MQ=255)
                                       aaaCACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCaaa                                                                                                       <  5:163562/148‑1 (MQ=255)
                                            ccGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGTGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGt                                                                                                 <  1:266890/149‑1 (MQ=255)
                                            ccGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGt                                                                                                 <  4:292935/149‑1 (MQ=255)
                                              ggCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTgg                                                                                               <  1:346126/149‑1 (MQ=255)
                                               gCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTgg                                                                                               <  4:382886/148‑1 (MQ=255)
                                               gCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTgg                                                                                               >  3:382886/1‑148 (MQ=255)
                                                  cGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTg                                                                                                <  6:72363/144‑1 (MQ=255)
                                                  cGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTg                                                                                                >  5:72363/1‑144 (MQ=255)
                                                   gTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGCGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGtt                                                                                          >  8:353917/1‑149 (MQ=255)
                                                   gTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGtt                                                                                          <  2:310301/149‑1 (MQ=255)
                                                    tttGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTc                                                                                         <  7:171728/149‑1 (MQ=255)
                                                       ggCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGa                                                                                       <  1:40659/148‑1 (MQ=255)
                                                       ggCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGa                                                                                       <  7:167861/148‑1 (MQ=255)
                                                              aaaGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATgccg                                                                               >  2:225355/1‑149 (MQ=255)
                                                               aaGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATgccgc                                                                              <  1:225355/149‑1 (MQ=255)
                                                                aGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATgccgcc                                                                             <  5:130525/149‑1 (MQ=255)
                                                                          gatgatTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAgg                                                                   <  7:281904/149‑1 (MQ=255)
                                                                          gatgatTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAgg                                                                   <  2:325897/149‑1 (MQ=255)
                                                                              atTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTa                                                                <  1:322477/148‑1 (MQ=255)
                                                                                tGCAAAAGAAGCCCGTCAGGCCCTGGCGGTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTAtt                                                              >  1:223788/1‑148 (MQ=255)
                                                                                  cAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGa                                                            <  7:353917/148‑1 (MQ=255)
                                                                                   aaaaGAAGCCCGTCAGGCCCTGGCGCTTCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCTTGAGCAggaa                                                                                                              <  3:67356/97‑3 (MQ=255)
                                                                                   aaaaGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGcc                                                                                                              >  4:67356/1‑97 (MQ=255)
                                                                                   aaaaGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGc                                                          >  3:270509/1‑149 (MQ=255)
                                                                                         agCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTggg                                                    <  1:177753/149‑1 (MQ=255)
                                                                                          gCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTgg                                                     >  7:211535/1‑147 (MQ=255)
                                                                                              gTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCCAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTg                                               <  2:104135/149‑1 (MQ=255)
                                                                                                    cccTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTg                                         <  2:17784/149‑1 (MQ=255)
                                                                                                     ccTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGa                                                                           >  7:265331/1‑114 (MQ=255)
                                                                                                     ccTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGa                                                                           <  8:265331/114‑1 (MQ=255)
                                                                                                         gcgcTGCCTCGGTTTGCCGAATACAATCAGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGTGCGTTCGAGGATGCCGCCGAGGAACAGTGTATTGAGCCTTGGGAGCTGGCGCTGCaa                                      <  8:328344/147‑1 (MQ=255)
                                                                                                          cgcTGCAGCGCTTTGACGAAGCCAATCCGCAGTTGCTCGATTAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGTCGTTCGAGGATGCCTCCGAGGAACAGGGTATTGAGCCTTGTGAGCTGGCTCTGCAattg                                   <  8:211535/149‑1 (MQ=255)
                                                                                                                         ccGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCg                    >  3:113104/1‑149 (MQ=255)
                                                                                                                          cgcaGCCACTCCGCAGTTGCTCGAGGAGCTGCGCGAGCTGGATGGGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCgg                   <  4:301254/146‑1 (MQ=255)
                                                                                                                                 aaTCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCGGAGCAGCt            >  7:147456/1‑149 (MQ=255)
                                                                                                                                        aGTTGCTCGTGTACATCCGAGAGCTGGATGCGCGTGAGCAGGCCCAGCTAATCGTGTGGGCGTTCGAGGATGCCTCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCTGAGCAGCTGCAgg       <  8:148570/147‑1 (MQ=255)
                                                                                                                                          ttGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCGGAGCAGCTGCAGGCCAt   >  3:274019/1‑149 (MQ=255)
                                                                                                                                          cagCTCGCGGATCTCCTCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTCGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAattgatt                                >  7:180753/3‑120 (MQ=255)
                                                                                                                                          cagCTCGCGGATCTCCTCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTCGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAattgatt                                <  8:180753/118‑1 (MQ=255)
                                                                                                                                           agCTCGCGGATCTCCTCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCGGAGCAGCTGCAGGCCATg  >  8:47375/2‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GTGGCACCCCAATAAAAGCTGTTAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCGGAGCAGCTGCAGGCCATG  >  NC_002947/2626311‑2626598

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