Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,428,977 C→G 8.6% P360A (CCA→GCA) ‡ mqo‑II → malate:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,428,9770CG 8.6% 179.6 / 7.2 58P360A (CCA→GCA) ‡mqo‑IImalate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (27/26);  new base G (3/2);  total (30/28)
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.22e-01

GACGGCAAGCGCGTCATCCTGTTTGGCCCGTTCGCTACCTTCTCGACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGGCGCCTGTGGCAGGCCGGCCAGCGCGTGCAAATCATCAAACGTGATGCC  >  NC_002947/1428836‑1429123
                                                                                                                                             |                                                                                                                                                  
gACGGCAAGCGCGTCATCCTGTTTGGCCCGTTCGCTACCTTCTCGACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTg                                                                                                                                               <  2:334304/147‑1 (MQ=255)
 aCGGCAAGCGCGTCATCCTGTTTGGCCCGTTCGCTACCTTCTCGACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGCGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTc                                                                                                                                            >  6:120446/1‑149 (MQ=255)
   ggCAAGCGCGTCATCCTGTTTGGCCCGTTCGCTACCTTCTCGACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGa                                                                                                                                          <  5:120446/149‑1 (MQ=255)
               aTCCTGTTTGGCCCGTTCGCTACGTTCTCGACCAAGTTCCTGTAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTg                                                                                                                               <  5:191035/148‑1 (MQ=255)
                 ccTGTTTGGCCCGTTCGCTACCTTCTCGACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGc                                                                                                                             >  8:140136/1‑148 (MQ=255)
                 ccTGTTTGGCCCGTTCGCTACCTTCTCGACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGc                                                                                                                             >  3:63530/1‑148 (MQ=255)
                   tGTTTGGCCCGTTCGCTACCTTCTCGACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTgg                                                                                                                                              <  6:123334/129‑1 (MQ=255)
                   tGTTTGGCCCGTTCGCTACCTTCTCGACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTgg                                                                                                                                              >  5:123334/1‑129 (MQ=255)
                           ccGTTCGCTACCTTCTCGACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTAc                                                                                                                                      >  7:56520/1‑129 (MQ=255)
                           ccGTTCGCTACCTTCTCGACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTAc                                                                                                                                      <  8:56520/129‑1 (MQ=255)
                                 gCTACCTTCTCGACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGt                                                                                                                <  5:170035/145‑1 (MQ=255)
                                 gCTACCTTCTCGACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGt                                                                                                                >  6:170035/1‑145 (MQ=255)
                                           cGACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTGGGTACTGGTCGatac                                                                                                                                       >  3:261146/1‑109 (MQ=255)
                                           cGACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTGGGTACTGGTCGatac                                                                                                                                       <  4:261146/112‑4 (MQ=255)
                                            gACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCt                                                                                                 >  5:366317/1‑149 (MQ=255)
                                            gACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACAATGACCg                                                                                                   <  8:47131/147‑1 (MQ=255)
                                                aaGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCg                                                                                              >  6:358766/1‑148 (MQ=255)
                                                 aGTTCCTGAAGAACTGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCTACCAGTACCCACTTGTCGAGTACCTCGCTGGCCAGCTTATGCTGTCTGACGATGACCGCTTCGaa                                                                                            <  1:92912/149‑1 (MQ=255)
                                                        gaagaaCGGCTCGTACCCGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCAACTGGTCGAGTACCTCGTTGGCCAGCTGCTGCTGTCTGACGATGACCGCTTCGACGCCCTgc                                                                                     >  8:327262/1‑149 (MQ=255)
                                                            aaCGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCAc                                                                                 <  6:366317/149‑1 (MQ=255)
                                                                    gTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTcc                                                                         <  5:358766/149‑1 (MQ=255)
                                                                    gTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTcc                                                                         <  5:335719/149‑1 (MQ=255)
                                                                    gTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTc                                                                          >  3:142137/1‑148 (MQ=255)
                                                                    gTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTc                                                                          <  4:142137/148‑1 (MQ=255)
                                                                       ccTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCg                                                                       >  5:186909/1‑148 (MQ=255)
                                                                        cTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGACGCCCTGCGCACCTCCGTCCCGaa                                                                     >  8:322427/1‑149 (MQ=255)
                                                                         tGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAAc                                                                    >  7:46491/1‑149 (MQ=255)
                                                                              ctgctgAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGt                                                                                                                <  8:298427/100‑1 (MQ=255)
                                                                              ctgctgAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGt                                                                                                                >  7:298427/1‑100 (MQ=255)
                                                                                gctgAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCaa                                                               <  7:242735/147‑1 (MQ=255)
                                                                                gctgAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCaa                                                               >  8:242735/1‑147 (MQ=255)
                                                                                     gcagcaCCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCaagaaag                                                          >  4:98487/1‑147 (MQ=255)
                                                                                     gcagcaCCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAagaaagaa                                                        <  8:46491/149‑1 (MQ=255)
                                                                                      cagcaCCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCaagaaag                                                          >  2:357932/1‑146 (MQ=255)
                                                                                         caccaccacCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCtt                                                                                                <  7:4128/105‑1 (MQ=255)
                                                                                         caccaccacCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCtt                                                                                                >  8:4128/1‑105 (MQ=255)
                                                                                          accaccaccCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATg                                                   <  3:98487/149‑1 (MQ=255)
                                                                                           ccaccaccCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATgg                                                  <  7:140136/149‑1 (MQ=255)
                                                                                            caccaccCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGgc                                                 <  4:63530/149‑1 (MQ=255)
                                                                                            caccaccCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGgc                                                 >  8:334200/1‑149 (MQ=255)
                                                                                             accaccCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGgcg                                                >  3:95902/1‑149 (MQ=255)
                                                                                                  ccacAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGTTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGGCGCCtgt                                           >  6:133181/1‑149 (MQ=255)
                                                                                                  ccacAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCaaga                                                             >  3:109885/1‑131 (MQ=255)
                                                                                                  ccacAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCaaga                                                             <  4:109885/131‑1 (MQ=255)
                                                                                                    acaACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGGCGCCTGTgg                                         >  5:328242/1‑149 (MQ=255)
                                                                                                      aaCGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGGCGCCTGTGGCa                                       <  6:328242/149‑1 (MQ=255)
                                                                                                      aaCGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGGCGCCTGTGGCa                                       <  6:186909/149‑1 (MQ=255)
                                                                                                                    gACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGt                                                                                                                                        <  6:81512/38‑1 (MQ=255)
                                                                                                                    gACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGt                                                                                                                                        >  5:81512/1‑38 (MQ=255)
                                                                                                                      ccAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGGCGCCTGTGGCAGGCCGGCCAGCGCGTg                       >  8:89721/1‑149 (MQ=255)
                                                                                                                      ccAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGGCGCCTGTGGCAGGCCGGCCAGCGCGTg                       <  1:357932/149‑1 (MQ=255)
                                                                                                                       cAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGGCGCCTGTGGCAGGCCGGCCAGCGCGTGc                      <  7:334200/149‑1 (MQ=255)
                                                                                                                                  tCGACCAGTGGGTACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAAt                                                    <  2:305848/108‑1 (MQ=37)
                                                                                                                                  tCGACCAGTGGGTACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAAt                                                    >  1:305848/1‑108 (MQ=37)
                                                                                                                                  tCGACCAGTGGGTACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGGCGCCTGTGGCAGGCCGGCCAGCGCGTGCAAATCATCaaa           >  2:259518/1‑149 (MQ=255)
                                                                                                                                     aCCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTACGAAGCCCTGCGCACCTa                                                                              <  4:233767/79‑1 (MQ=255)
                                                                                                                                     aCCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTACGAAGCCCTGCGCACCTa                                                                              >  3:233767/1‑79 (MQ=255)
                                                                                                                                      ccAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTTGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGGCGCCTGTGGCAGGCCGGCCAGCGCGTGCAAATCATCAAACGt        >  5:115968/1‑148 (MQ=255)
                                                                                                                                            cccACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGGCGCCTGTGGCAGGCCGGCCAGCGCGTGCAAatca                <  4:130693/134‑1 (MQ=255)
                                                                                                                                            cccACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGGCGCCTGTGGCAGGCCGGCCAGCGCGTGCAAatca                >  3:130693/1‑134 (MQ=255)
                                                                                                                                            cccACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGGCGCCTGTGGCAGGCCGGCCAGCGCGTGCAAATCATCAAACGTGATGcc  >  8:118545/1‑148 (MQ=255)
                                                                                                                                             ccACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGGCGCCTGTGGCAGGCCGGCCAGCGCGTGCAAATCATCAAACGTGATGcc  >  3:80428/1‑147 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GACGGCAAGCGCGTCATCCTGTTTGGCCCGTTCGCTACCTTCTCGACCAAGTTCCTGAAGAACGGCTCGTACCTGGACCTGCTGAGCAGCACCACCACCCACAACGTGTGGCCGATGACCAAAGTCGGTATCGACCAGTACCCACTGGTCGAGTACCTCGCTGGCCAGCTGATGCTGTCTGACGATGACCGCTTCGAAGCCCTGCGCACCTACTTCCCGAACGCCAAGAAAGAAGAATGGCGCCTGTGGCAGGCCGGCCAGCGCGTGCAAATCATCAAACGTGATGCC  >  NC_002947/1428836‑1429123

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