Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,493,695 A→C 10.3% L184R (CTG→CGG) ‡ PP_3097 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,493,6950AC10.3% 153.0 / 8.0 58L184R (CTG→CGG) ‡PP_3097hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/23);  new base C (3/3);  total (32/26)
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.65e-01

TTCAGGGCATCGGCCGGCAGCGCCACCTCGCGGCTGCGGCCCTGGTCGTCGAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCATTTCTGTGAACGAGGCCTGCTGCACGGCAATCGGC  >  NC_002947/3493557‑3493838
                                                                                                                                          |                                                                                                                                               
ttCAGGGCATCGGCCGGCAGCGCCACCTCGCGGCTGCGGCCCTGGTCGTCGAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATg                                                                                                                                      <  1:1376383/150‑1 (MQ=255)
 tCAGGGCATCGGCCGGCAGCGCCACCTCGCGGCTGCGGCCCTGGTCGTCGAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCCTGCGGTTcttcc                                                                                                                                     <  1:796728/150‑5 (MQ=255)
 tCAGGGCATCGGCCGGCAGCGCCACCTCGCGGCTGCGGCCCTGGTCGTCGAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCCTGCGGTTcttcc                                                                                                                                     <  1:467220/150‑5 (MQ=255)
         tCGGCCGGCAGCGCCACCTCGCGGCTGCGGCCCTGGTCGTCGAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGt                                                                                                                                     >  1:389068/1‑142 (MQ=255)
         tCGGCCGGCAGCGCCACCTCGCGGCTGCGGCCCTGGTCGTCGAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGt                                                                                                                                     <  2:389068/142‑1 (MQ=255)
           ggccggcAGCGCCACCTCGCGGCTGCGGCCCTGGTCGTCGAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGccc                                                                                                                           >  2:731033/1‑150 (MQ=255)
                 cAGCGCCACCTCGCGGCTGCGGCCCTGGTCGTCGAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAg                                                                                                                     <  2:743117/150‑1 (MQ=255)
                   gcgcCACCTCGCGGCTGCGGCCCTGGTCGTCGAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAggg                                                                                                                   >  2:836151/1‑150 (MQ=255)
                       cACCTCGCGGCTGCGGCCCTGGTCGTCGAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATc                                                                                                               >  1:99147/1‑150 (MQ=255)
                          cTCGCGGCTGCGGCCCTGGTCGTCGAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGcc                                                                                                            >  1:557013/1‑150 (MQ=255)
                                        ccTGGTCGTCGAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGc                                                                                              <  2:288547/150‑1 (MQ=255)
                                        ccTGGTCGTCGAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAAAAGTTCATCCAGCTGCAGCATCAGGGTTGCGTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGc                                                                                              >  1:174915/1‑150 (MQ=255)
                                          tGGTCGTCGAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCgg                                                                                            >  1:1344806/1‑150 (MQ=255)
                                           ggTCGTCGAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGa                                                                                           <  2:557013/150‑1 (MQ=255)
                                               gtcgAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATgcg                                                                                       >  2:1279372/1‑150 (MQ=255)
                                                   aaGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGa                                                                                   <  1:1002554/150‑1 (MQ=255)
                                                    aGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGAc                                                                                  <  1:473122/150‑1 (MQ=255)
                                                      cTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTNCAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGa                                                                                >  1:1388898/1‑150 (MQ=255)
                                                      cTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGa                                                                                <  1:625693/150‑1 (MQ=255)
                                                       tCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAg                                                                               <  1:1055836/150‑1 (MQ=255)
                                                        cAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGc                                                                              >  1:350998/1‑150 (MQ=255)
                                                          gCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTgg                                                                                                       <  1:911902/123‑1 (MQ=255)
                                                          gCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTgg                                                                                                       >  2:911902/1‑123 (MQ=255)
                                                          gCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGCTGCGCAGACGAGCGa                                                                            >  2:1076717/1‑150 (MQ=255)
                                                                ccTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAgg                                                                                                                    >  2:754129/1‑104 (MQ=255)
                                                                ccTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAgg                                                                                                                    <  1:754129/104‑1 (MQ=255)
                                                                 cTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTc                                                                     >  1:496138/1‑150 (MQ=255)
                                                                       gAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCt                                                               <  1:836151/150‑1 (MQ=255)
                                                                             gttgAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGccc                                                         >  1:345554/1‑150 (MQ=255)
                                                                                    agcagTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAg                                                  >  2:417959/1‑150 (MQ=255)
                                                                                      cagTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGc                                                >  1:756943/1‑150 (MQ=255)
                                                                                        gTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCg                                              >  1:977535/1‑150 (MQ=255)
                                                                                         ttCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGa                                             <  1:1200372/150‑1 (MQ=255)
                                                                                          tCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAg                                            <  2:74527/150‑1 (MQ=255)
                                                                                                aGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTcc                                      >  2:1265723/1‑150 (MQ=255)
                                                                                                   tGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAgg                                                                                                                    >  2:393129/1‑69 (MQ=255)
                                                                                                   tGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAgg                                                                                                                    <  1:393129/69‑1 (MQ=255)
                                                                                                      agcaTCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCATTTCt                                >  1:1406118/1‑150 (MQ=255)
                                                                                                        caTCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCATTTCtgt                              <  2:345554/150‑1 (MQ=255)
                                                                                                                 cGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCATTTCTGTGAACGAGGc                     <  1:417959/150‑1 (MQ=255)
                                                                                                                    ttCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATgc                                                                                        >  1:911307/1‑80 (MQ=255)
                                                                                                                    ttCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATgc                                                                                        <  2:911307/80‑1 (MQ=255)
                                                                                                                    ttCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCATTTCTGTGAACGGGGCctg                  >  1:601693/1‑150 (MQ=255)
                                                                                                                          gTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAg                                                  <  2:989000/112‑1 (MQ=255)
                                                                                                                          gTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAg                                                  >  1:989000/1‑112 (MQ=255)
                                                                                                                          gTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCATTTCTGTGAACGAGGCCTGCTGCAc            >  1:508386/1‑150 (MQ=255)
                                                                                                                          gTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCATTTCTGTGAACGAGGCCTGCTGCAc            >  1:210115/1‑150 (MQ=255)
                                                                                                                              gacatgaaCCGCCTGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCATTTCTGTGAACGAGGCCTGCTGCACGGCa        >  1:345865/6‑150 (MQ=255)
                                                                                                                                gaagaaCCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCATTTCTGTGAACGAGGCCTGCTGCACGGCAAt      <  1:638522/150‑1 (MQ=255)
                                                                                                                                catgaaCCGCCTGCGGTTCATGTCCATCTGCCCGAACAGTGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCa                                     <  2:524484/116‑1 (MQ=255)
                                                                                                                                catgaaCCGCCTGCGGTTCATGTCCATCTGCCCGAACAGTGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCa                                     >  1:524484/4‑119 (MQ=255)
                                                                                                                                catgaaCCGCCTGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCATTTCTGTGAACGAGGCCTGCTGCACGGCAAt      >  1:717232/4‑150 (MQ=255)
                                                                                                                                  agaaCCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCATTTCTGTGAACGAGGCCTGCTGCACGGCAATCg    <  1:417058/150‑1 (MQ=255)
                                                                                                                                  agaaCCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCATTTCTGTGAACGAGGCCTGCTGCACGGCAATCg    >  2:609400/1‑150 (MQ=255)
                                                                                                                                  agaaCCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCATTTCTGTGAACGAGGCCTGCTGCACGGCAATCg    >  1:338424/1‑150 (MQ=255)
                                                                                                                                    aaCCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCATTTCTGTGAACGAGGCCTGCTGCACGGCAATCGGc  <  1:731033/150‑1 (MQ=255)
                                                                                                                                      ccGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGc                                           <  2:714401/107‑1 (MQ=255)
                                                                                                                                      ccGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGc                                           >  1:714401/1‑107 (MQ=255)
                                                                                                                                          |                                                                                                                                               
TTCAGGGCATCGGCCGGCAGCGCCACCTCGCGGCTGCGGCCCTGGTCGTCGAAGCTCAGCGTGGCCTTGGCGAGGTTGTTGAACAGCAGTTCGTACAGCTGCAGCATCAGGGTCGCTTCGCCGTCGAGGAAGAACCGCAGGCGGTTCATGTCCATCTGCCCGAACAGGGCATCGCCAGTGGCTGCCAGGCGGATGCGCAGACGAGCGACCAGGTCGGCGCTGTGCCCGTTGAAGGCCGAGCGTTCCATTTCTGTGAACGAGGCCTGCTGCACGGCAATCGGC  >  NC_002947/3493557‑3493838

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