Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913266,3490CA44.4% 99.5 / 31.4 55V69F (GTT→TTT) yafYlipoprotein, inner membrane; degP regulator; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base C (4/26);  new base A (24/0);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 1.48e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAT  >  NC_000913/266204‑266499
                                                                                                                                                   |                                                                                                                                                      
tCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GGAAAAACTGc                                                                                                                                                    >  1:347162/1‑150 (MQ=255)
     cACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTg                                                                                                                                              <  1:654718/151‑1 (MQ=255)
     cACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GTAAAAACTGCAAgggg                                                                                                                                              >  1:136460/1‑149 (MQ=255)
            tGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GTAAAAACCGCCAAGGGGGGCCgg                                                                                                                                       >  2:468891/1‑151 (MQ=255)
              gTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGa                                                                                                                                      <  2:735653/150‑1 (MQ=255)
              gTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGGC‑GTAAAAACCGCAAGGTGGGGCCgg                                                                                                                                       >  1:368760/1‑149 (MQ=255)
                aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAAACTGCAAGGGGGGGCCGGATa                                                                                                                                    >  1:532513/1‑150 (MQ=255)
                  tCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GGAAAAACCGCAAAGGGGGGCCgg                                                                                                                                       >  2:777144/1‑145 (MQ=255)
                   ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAagt                                                                                                                                 <  2:150616/150‑1 (MQ=255)
                   ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAagt                                                                                                                                 <  2:945274/150‑1 (MQ=255)
                    cAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GTAAAACCCGCAAGGGGGGGCCggggaa                                                                                                                                   >  1:449654/1‑143 (MQ=255)
                       tttGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTACGGTAAAAGTTGCGGGGGGGGTCCggcgagcg                                                                                                                                 >  2:397377/1‑141 (MQ=255)
                         tGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCaa                                                                                                                           <  1:288156/150‑1 (MQ=255)
                         tGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCaa                                                                                                                           <  1:840724/150‑1 (MQ=255)
                               tATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCtgtg                                                                                                                     <  1:94214/150‑1 (MQ=255)
                                  cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCa                                                                                                                 <  1:257814/151‑1 (MQ=255)
                                  cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GGAAAAACCGCAAGGGGGGGCCGGGTAAGGAGGCAAACTGGGaac                                                                                                                  >  2:612510/1‑148 (MQ=255)
                                     cAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAAACTGCAAGGGGGGGCCGGATAAAGAGTCAACCTGTGACCaat                                                                                                               >  1:495315/1‑148 (MQ=255)
                                      aaTACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCAttct                                                                                                             <  2:669612/151‑1 (MQ=255)
                                                 cATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGAcatc                                                                                                  <  1:202315/151‑1 (MQ=255)
                                                    tcatcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAA‑GTAAAAACCGCCAGGGGGGACCCGGTAAAAAGTCAAACCGTGGCCCATCTTCGCCGAcatca                                                                                                 >  1:657127/1‑149 (MQ=255)
                                                      atcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCt                                                                                             <  2:619993/151‑1 (MQ=255)
                                                               tGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGGA‑GGAAAAACAGCAAGGCGGGGACGGGCAAGGAGTGCACCTGTGCCCATTCTTTACTGTCATAATCATGGGGGCTg                                                                                     >  2:829255/1‑150 (MQ=255)
                                                                ggACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GTAAAAACCGCAAGGGGGGGCCGGGTAAGTAGTCAACCCGTGACCATTCTTTGCTGACATCATCCTGGGGGGTg                                                                                     >  2:783671/1‑149 (MQ=255)
                                                                      ttttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GTAAAAACTTCAAAGTGGGGCCGGACAAAAAGGCAACCCTTGTCCATTCTTTGCTGGTATTATCCCGGTTGCTGa                                                                                    >  1:305332/1‑144 (MQ=255)
                                                                         tttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAA‑AGAAAAACCGCAAGGTGGGGCCGGATAAGGAAACAACCTTTGACCCTTTTTTGCTGTCAATATCCCGGGGGGTGGGCACATc                                                                             >  1:141953/1‑148 (MQ=255)
                                                                          ttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAAACTGCAAGGGGGGGCCGGGTAAAAAGGCAACCTGTGACCATTTTTCGCTGACATCATCCTTGGGGCTGAcccc                                                                                >  1:624432/1‑142 (MQ=255)
                                                                            ccTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GTAAAAACTGCAAGGGGGGGCCGGATAAAGAGGCAAACTGTGACCATTTTTTCCTGACATCATCCTGGTGGCGGGCCACCTCACtta                                                                        >  1:171347/1‑148 (MQ=255)
                                                                              tGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCTGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGaa                                                                       <  2:253538/149‑1 (MQ=255)
                                                                               gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGaaa                                                                      <  2:935844/149‑1 (MQ=255)
                                                                                cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GTAAAAACTGCAAGGGGGGGCCGGGGAAGTAGTCAACCTGTGACCATTCTTCGCCTACATCATCCTCGTGGGTGACCACATCACCGaaa                                                                      >  1:562725/1‑148 (MQ=21)
                                                                                     ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTg                                                               <  2:436110/150‑1 (MQ=255)
                                                                                       aCGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAA‑GTAAAAACCGCCAAGGGGGACCGGGAAAGTAGTCCACCTGTGGCCATTCCTCCCTGACCTCATCCTGGTGGGTGACCCCATAACTCAAACGCGGCgggg                                                            >  2:191775/1‑151 (MQ=255)
                                                                                        cGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGt                                                           <  2:268869/151‑1 (MQ=255)
                                                                                        cGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGt                                                           <  2:851356/151‑1 (MQ=255)
                                                                                              agtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAA‑GTAAAAACGGCAAGGGGGGGCCCGGTAAGAAGTCCAACTGTGACCATTCTTCGCTGACATCCCCCTGGCGGCAGACCACATCACTTAAAAGCGgggggggg                                                          >  1:351168/1‑144 (MQ=12)
                                                                                               gtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTcc                                                    <  2:950856/151‑1 (MQ=255)
                                                                                                  gtaACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAAACCGCAAGGGGGGACCGGATAAGGAGTCAACCCGTGACCATTCTTCGCTGACATCATCCTGGTGGGTGACCACATCACTGAAACGCGGGGGGGGGAGGTCCttt                                                 >  1:581822/1‑150 (MQ=18)
                                                                                                  gtaACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GGAAAAACCGCAAGGTGGGGCCGGGTAAAGAGTCAAACCGTGACCATTCTTCGCTGACATCATCCTGGGGGCGGACCACATCACTGAAACGCGGGGGGGGGGGGGCCt                                                   >  2:172290/1‑149 (MQ=18)
                                                                                                    aaCTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAA‑GTAAAAACCGCAAGGTGGGGCCCGCTAAGAAGTCAACCTTTGACCCTTCTTCTCTTACATCATCCCCGTGGGTGGCCACCACACTTAAAAGCGGTGGGGGGGGGGCCttta                                                >  2:588874/1‑148 (MQ=255)
                                                                                                     aCTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAAACTGCAAAGGGGGACCGGATAAATAATCAACCCGTGAACACTTTTTCCTGACATCATCCTGGGGGGTGACCAAATCACTGAAACGCGGGGGGGGGAGGGCCTTaaa                                               >  1:787370/1‑150 (MQ=255)
                                                                                                       tGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GGAAAAACTGCAAGGGGGGGCCGGATAAATAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGGGGGTGGCCACATCACTGAAACGCGGGGGGGGGGGGGCCTTTAAttt                                            >  2:358009/1‑151 (MQ=18)
                                                                                                         ccATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTAc                                          <  1:231285/151‑1 (MQ=255)
                                                                                                         ccATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTAc                                          <  1:201791/151‑1 (MQ=255)
                                                                                                             ccTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGAt                                      <  2:406526/151‑1 (MQ=255)
                                                                                                             ccTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGAt                                      <  2:788675/151‑1 (MQ=255)
                                                                                                              cTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAAACTGCAAGGTGGGGCCCGATAAGTAGTCAACCTGTGACCATTCTTCTCTGACATCATCCTGGGGGCTGACCACATCACTGAAAAGCGGTGGGGGGAGGGCCTTTAATTTACTGAt                                      >  2:622324/1‑150 (MQ=255)
                                                                                                                   aGCCATCGGGCGACTGGCGGGGCAG‑ATAAAACC‑GCAAAGGGGGGCCGGATAAAAAGTTAACCCGTGTCCACTCTTCGCTGACATCCTTCTGGTGGGTGACCACCTCAACGAAACACGGAGGGGGGAGGGGCTTAAAATTAATGATGGCCtt                                >  2:961345/1‑151 (MQ=255)
                                                                                                                   aGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCtt                                <  2:814675/151‑1 (MQ=255)
                                                                                                                         cGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtc                          <  2:722556/151‑1 (MQ=255)
                                                                                                                                  ggC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTc                  <  2:622846/150‑1 (MQ=255)
                                                                                                                                  ggC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCt                 <  1:603455/151‑1 (MQ=255)
                                                                                                                                    c‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGt               >  2:198156/1‑151 (MQ=255)
                                                                                                                                       ggTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCa             <  2:431756/151‑1 (MQ=255)
                                                                                                                                         tAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATc           >  1:250643/1‑151 (MQ=255)
                                                                                                                                                   ccTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAt  >  2:826254/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAT  >  NC_000913/266204‑266499

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