Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 272,201 A→G 6.7% intergenic (+12/‑379) insN → / → eyeA pseudogene, IS911 transposase A;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related/novel sRNA, function unknown, CP4‑6 propahge

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913272,2010AG 6.7% 177.5 / 3.2 60intergenic (+12/‑379)insN/eyeApseudogene, IS911 transposase A;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related/novel sRNA, function unknown, CP4‑6 propahge
Reads supporting (aligned to +/- strand):  ref base A (24/31);  new base G (2/2);  total (27/33)
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.40e-01

GTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGAT  >  NC_000913/272065‑272349
                                                                                                                                        |                                                                                                                                                    
gTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAgggcggg                                                                                                                                          <  5:126021/149‑1 (MQ=255)
gTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAgggcggg                                                                                                                                          >  7:332310/1‑149 (MQ=255)
 atGCCAGCAGCCGACACACCGGTATAAGCGTGGCGGTCATGAGCACGTGGCTATCACGAATCATATTGAGCGACAGTTCGCCGTAACGGAACCAAAACAGGGGTGGTGCGGGGATGTGACCAATAGTGTGGCCGGCGTTCAgggcggg                                                                                                                                          >  6:210045/2‑148 (MQ=255)
    ccAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATg                                                                                                                                      <  8:41002/149‑1 (MQ=255)
        cagcCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATg                                                                                                                                  <  8:307862/149‑1 (MQ=255)
          gcCGACTCACCGGTATAAGCTGGGCGGTCACGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGCTCGCCGTACCGGAACCACATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCt                                                                                                                                <  6:110153/149‑1 (MQ=37)
              aCTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTaaa                                                                                                                            <  6:52097/149‑1 (MQ=255)
                tCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAAt                                                                                                                           >  5:286577/1‑148 (MQ=255)
                  acaggTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAAACAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGCAGTTCAGGGCGGGCATGGATGCTTAAATGaa                                                                                                                        >  3:269038/4‑149 (MQ=255)
                      gTATAAGCGTGGCGGTCATTAGCACGTTGCTATCCCGAATCATCTTTATCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACcgc                                                                                                                    <  6:286577/149‑1 (MQ=255)
                           aGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAgt                                                                                                                <  7:199558/148‑1 (MQ=255)
                            gCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGACGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAgtct                                                                                                              >  3:97780/1‑149 (MQ=255)
                               tGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAgtctgtc                                                                                                           <  6:107154/149‑1 (MQ=255)
                                        tGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATAtt                                                                                                  <  7:86973/149‑1 (MQ=255)
                                        tGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGACGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATAtt                                                                                                  <  4:97780/149‑1 (MQ=255)
                                           gCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATGGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGaa                                                                                               <  4:78131/149‑1 (MQ=255)
                                           gCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGaa                                                                                               <  6:118693/149‑1 (MQ=255)
                                                  gCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTAAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTATGGAATATTGAACCGGTaa                                                                                        >  7:161511/1‑149 (MQ=255)
                                                        ccGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGCAATATTGAACCGGTAACTCAc                                                                                   >  2:95788/1‑148 (MQ=255)
                                                         cGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGa                                                                                 <  8:332310/149‑1 (MQ=255)
                                                         cGAATCATCTTGAGCGACAGTGCGCCGTAACGGAACCAAATCAGGTGTGGGGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGa                                                                                 <  5:210045/149‑1 (MQ=255)
                                                          gAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGAt                                                                                >  6:34387/1‑149 (MQ=255)
                                                          caaTCATCTTGAGCGACAGTTCGCCGTACCGGAACCAACTCCCGTGTGGTGCGGTGATGTCACCTATAGTGTGCCCGGAGTGCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTCGAATATTGAACCGGTAACTCACGAt                                                                                <  6:253893/148‑1 (MQ=255)
                                                             tcatcTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGAt                                                                                >  5:286892/1‑146 (MQ=255)
                                                                   tGAGCGACTGTTCGACGTAACGTAACCATATCAGTGGTGGTGCGGTGATGTGAACTATTGTGTGCCAGTAGTTGATGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAAcc                                                                        <  8:161511/148‑1 (MQ=255)
                                                                       cGACAGTTGGGCGTAACGGAACCATATCTGGTTTGGTGCGGTGTTGTGACCTAAAGTGTGCGCGGAGTTCAGGGGGGGCATGGTTTCTTTAATGAACCGCGATTCTGTCTGGAATATTGAAGCTGTAACTCTCGATGAGAAACCcaaca                                                                   <  4:269038/149‑1 (MQ=255)
                                                                         aCAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCcaa                                                                     >  1:217858/1‑145 (MQ=255)
                                                                         aCAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCcaa                                                                     <  2:217858/145‑1 (MQ=255)
                                                                             ttCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTCTAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCa                                                             <  6:286892/149‑1 (MQ=255)
                                                                                   gTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCAc                                                                                   >  2:314994/1‑121 (MQ=255)
                                                                                     aaCGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAat                                                                                                     >  6:175319/1‑101 (MQ=255)
                                                                                     aaCGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAat                                                                                                     <  5:175319/101‑1 (MQ=255)
                                                                                        ggAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGAc                                                  >  3:174326/1‑149 (MQ=255)
                                                                                            ccAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTg                                              >  3:186751/1‑149 (MQ=255)
                                                                                                aTCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGgtgt                                                     <  5:81992/138‑1 (MQ=255)
                                                                                                aTCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGgtgt                                                     >  6:81992/1‑138 (MQ=255)
                                                                                                 tCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGaa                                         <  5:34387/149‑1 (MQ=255)
                                                                                                   aGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAAcc                                       >  3:254990/1‑149 (MQ=255)
                                                                                                        tggtgCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGc                                  >  2:358332/1‑149 (MQ=255)
                                                                                                               gTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCAc                                                            <  8:37999/116‑1 (MQ=255)
                                                                                                               gTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCAc                                                            >  7:37999/1‑116 (MQ=255)
                                                                                                                 gATGTGACCTATAGTGTGCCCGGGGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTg                                              <  4:214380/128‑1 (MQ=255)
                                                                                                                 gATGTGACCTATAGTGTGCCCGGGGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTg                                              <  4:209547/128‑1 (MQ=255)
                                                                                                                 gATGTGACCTATAGTGTGCCCGGGGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTg                                              >  3:209547/1‑128 (MQ=255)
                                                                                                                 gATGTGACCTATAGTGTGCCCGGGGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTg                                              >  3:214380/1‑128 (MQ=255)
                                                                                                                  aTGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGAc                         <  1:95788/148‑1 (MQ=255)
                                                                                                                          tataGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCgg                                                         >  6:149849/1‑108 (MQ=255)
                                                                                                                          tataGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCgg                                                         <  5:149849/108‑1 (MQ=255)
                                                                                                                            taGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGt                            >  6:340991/1‑135 (MQ=255)
                                                                                                                            taGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGt                            <  5:340991/135‑1 (MQ=255)
                                                                                                                               tgtgCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGc                                  <  2:204886/126‑1 (MQ=255)
                                                                                                                                gtgCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAgtctgtct                                                                                                          >  1:197058/1‑53 (MQ=255)
                                                                                                                                gtgCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAgtctgtct                                                                                                          <  2:197058/53‑1 (MQ=255)
                                                                                                                                    ccGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGga                                            <  4:115308/111‑1 (MQ=255)
                                                                                                                                    ccGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGga                                            >  3:115308/1‑111 (MQ=255)
                                                                                                                                    ccGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAggg      >  5:279835/1‑149 (MQ=255)
                                                                                                                                    ccGGAGTTCAGGGCGGGCATAGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGa                                    >  4:178200/1‑119 (MQ=255)
                                                                                                                                    ccGGAGTTCAGGGCGGGCATAGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGa                                    <  3:178200/119‑1 (MQ=255)
                                                                                                                                     cGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAATCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGt     >  7:357097/1‑149 (MQ=255)
                                                                                                                                     cGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGt     <  4:254990/149‑1 (MQ=255)
                                                                                                                                        aGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGAt  <  4:174326/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                    
GTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGAT  >  NC_000913/272065‑272349

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: