Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 78,524 A→C N302H (AAT→CAT)  setA → sugar exporter SetA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091378,5240AC94.4% 101.2 / ‑4.9 36N302H (AAT→CAT) setAsugar exporter SetA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (12/22);  minor base T (2/0);  total (14/22)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGACTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTCTTTAATAGCCGTATGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGT  >  NC_000913/78279‑78762
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               
tATTTGTCGCCTCGACGTTAATGTGGTCCTGCAACACCATGTACATTATTGATATGTCGTTGTGGATCATTAGCGAGTTAGGATTGCCAGACAAACGGGCGGGTTTCCTGATGGGGGCGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGGGGCGGCAGGAGTACTGTTTTACACCGGATTGATGTTATTTCATAg                                                                                                                                                                                                                                             <  1:185387/250‑1 (MQ=255)
tATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAg                                                                                                                                                                                                                                             <  2:227862/250‑1 (MQ=255)
         ccTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCg                                                                                                                                                                                                                                   <  1:120850/251‑1 (MQ=255)
                                 acacCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAAc                                                                                                                                                                                                          <  1:500516/251‑1 (MQ=255)
                                                 tGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGc                                                                                                                                                                                           <  1:11098/250‑1 (MQ=255)
                                                  gATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGGTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGGTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGGGGCCGGGGGACCGTTTTTTACCGGGGTGTTTTTTTTTTATATCCGTCTGGCGTTGTTGGCGCTGCACCttttttttgatgtttttttctggct                                                                                                                                                                                         >  2:142713/1‑231 (MQ=255)
                                                   atatGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGCGTTCCTTATGGGGACGGCAGCTGGTATTGTAACACCCGCATTGTTTCTGGCTAGCTCATTTGTCAAAAGTTATGGGCAGCGGCGAATGGTGGTTATTGCCGTGGCGGCCAGGGGACTGTTTTATACCCGATTTTTTTTATTTTTTAACCCTCCGGCGTTGATGGAGCTGCAACCTTTTTAACCCGTATTTTTCGGcctt                                                                                                                                                                                        >  1:689017/1‑248 (MQ=255)
                                                   aaatggcggtttggaataGTAGCGAGTTAGGATTGCCATGCAAACTGGTCGGTTTCCTGATGGGGTCGGCAGGTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCCAACGTTATGGTAAGCGGAGAATGATGGTCATAGCAGTGGCGGCAGGAGGACTGTTTTACCCCGTATTGATTTTATTTCATAGCCGTCTGGCGTTGAGGACGCGGCAACTTGTTAACTCTGTATGTATCGGCAt                                                                                                                                                                                         <  2:668647/233‑1 (MQ=255)
                                                    tatGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCCGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAAAGTTATGGTAAGCGGCGAATGATGGTCATAGCCGTGGCGGCAGGAGTAATGGTTTACACCGGGTTGATTTTATTTCATAGCCGGCTGGCGTTGATGACGCTGCAACTTTTTAAAGCTGTATTTATTGGCAttg                                                                                                                                                                                       >  1:530143/1‑251 (MQ=255)
                                                                gatccgtcgcgaggtgggATTGCCAGCCAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCGGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGTTAGGCGGCGAATGATGGTCATATCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATGGCATAGCCGTCTGGCGTTGATGACGCTGCAACTGTTTAACGCTGTATTTATCGTCATTGTTGCGGGTATTg                                                                                                                                                                           <  2:638322/235‑1 (MQ=255)
                                                                   cAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGGCATAGCAGGGGCGGCCGGAGGACTGTTTTACAACGGATTGATTTTTTTTCCTAGCCGTCTGGCGTTGATGAAGCGGCAACTTTTTTAAGATGTATTTATCGGGATTGTTGCggggtttgggg                                                                                                                                                                        >  1:530821/1‑243 (MQ=255)
                                                                          gAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATg                                                                                                                                                                 <  1:70197/251‑1 (MQ=255)
                                                                                      ccAGACAAACTGGCGGGGTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCCGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGGATGATGGTCATAGCAGTGGCGGCAGGAGTACTGGTTTTCCACGGATTTGTTTTATTTCATAGCCGTCTGGGGTTGGTGACGATGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGGATTTGGGTGCTATGGTTTTAGGAttt                                                                                                                                                     >  1:632108/1‑251 (MQ=255)
                                                                                                 ggCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGa                                                                                                                                           <  1:142713/250‑1 (MQ=255)
                                                                                                         tCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCggg                                                                                                                                   <  1:372415/250‑1 (MQ=255)
                                                                                                                ggggACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCt                                                                                                                            >  1:244571/1‑250 (MQ=255)
                                                                                                                     cGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTaccacc                                                                                                                      <  2:244571/251‑1 (MQ=255)
                                                                                                                         aGCTGGTCTGGAAATACCCGCACTGATTCTGTCTGGCTACTATGTCACACGTTATGGTACGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGGTTTTGTTTCCTAGCCGTCTGGCGTTGATGACGCTGCAAATTTTTTACGCTTTTTTTATCGGCATTGTTGCGGGTATTGGGATGCTTTGGTTTCAGGTTTTAATTCCTGGAAGAGCGGGGGCAGCCACCACCttct                                                                                                                  >  2:650345/1‑249 (MQ=255)
                                                                                                                                         ccAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATg                                                                                                                                                                 <  2:154266/188‑1 (MQ=255)
                                                                                                                                         ccAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATg                                                                                                                                                                 >  1:154266/1‑188 (MQ=255)
                                                                                                                                                               cctaTGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGGCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGGGCGGGGGCAGCTACCACCTTATTTACTCACAGTATTTATAACGGGGGAATTctggctggg                                                                            >  2:254903/2‑250 (MQ=255)
                                                                                                                                                                     tCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCCGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTAtt                                                                      >  1:463164/1‑251 (MQ=255)
                                                                                                                                                                         aCGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTa                                                                                                                                                    <  1:542167/169‑1 (MQ=255)
                                                                                                                                                                         aCGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCGGGCGTTGATGACGCTGCAACTTTTTAAAGCTGTATTTATCGGCATTGTTGCGGGGATTGGGTTGCTATGGTTTCATGATTTa                                                                                                                                                    >  2:542167/1‑169 (MQ=255)
                                                                                                                                                                            ttATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAggg                                                                 >  1:455920/1‑249 (MQ=255)
                                                                                                                                                                                   aaGCGTCGAATGAGGATCATAGCAGTGGCGGCAGGAGTACGGTTTTGCACCGTAGTGAGTTTATTTCATAGCCGTGAGGCGTTGATGACGCTGCAGCTTTTTAACGCTGTATTTATCGGCATTGTTGCGGTGATTGGGATGCTATGGTTTCATGAGTTAATGCCTGGAAGAGCGGGGGCAGCGACCACCTTATTGACTAACAGTATGTCTACCGGGGTAAGTCTGGCTGGCGTTATTCAGGGAGCAATTGc                                                        <  2:632116/251‑1 (MQ=255)
                                                                                                                                                                                   aaGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGc                                                        <  2:632108/251‑1 (MQ=255)
                                                                                                                                                                                        gcgAATGATGGTCATAGCAGTGGCTGCAGGAGTACTGTTTTTCACCGGATTGATTTTATTTCATAGCCGTCTGGCGGTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCATGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACaaa                                                   <  2:612309/251‑1 (MQ=255)
                                                                                                                                                                                                    cATAGCAGGGGCGGCAGGAGTACTGTTTTACACCGGGTTGATTTTATTTCATAGGCGTCGGGCGGTGATGACGCGGCAACGTTTTAACGCTGTATTTAGCGGCATTGTTGCGGGTAGTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTAGCACCTTATTTGCTAACAGGATTTCTACCGGGGTAAGTCTGGCTGGCGTTAGTCAGGGAGCAATTGCACAAAGGTGGGGGcgct                                       <  2:68166/251‑4 (MQ=255)
                                                                                                                                                                                                             ggcggcAGGATGTCTGGTTTACTCCGGATTGAGTTGATTTCATAGCCGTCTGGCTTTGATGACTCTGCAACTTTTTAACGCTGTTTTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGGGCGGGGGCAGCTACCACCTTATTTTCTAACAGTATTTCTACCGGGGTAATTCTGGCTGTCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTc                               <  2:643339/250‑1 (MQ=255)
                                                                                                                                                                                                                  cAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTgg                         >  1:570446/1‑251 (MQ=255)
                                                                                                                                                                                                                    gcaggaatactttgtttccccGGATTGTTTTTATTTCATAGCTGTCTGGCGTTGAGTACGCTTCAACTTTTTAACGCTGTATTTATCGTCTTTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGTCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGTGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACATTGCTGTCTACTg                          <  2:571465/232‑1 (MQ=255)
                                                                                                                                                                                                                                 tttACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGtt             >  1:512694/1‑251 (MQ=255)
                                                                                                                                                                                                                                   tACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTa            <  1:254903/250‑1 (MQ=255)
                                                                                                                                                                                                                                          gattgattTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGTTGCACTTTGCTGTCTACTTGGTAATTGCTGGTATTGCtgtt    <  2:665113/251‑1 (MQ=255)
                                                                                                                                                                                                                                             tgattTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCtgttat  <  1:533641/250‑3 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               
TATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGACTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTCTTTAATAGCCGTATGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGT  >  NC_000913/78279‑78762

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