Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,434,8820TG50.9% 66.5 / 90.6 56Y252S (TAC→TCC) truAtRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (12/15);  new base G (28/0);  total (41/15)
Fisher's exact test for biased strand distribution p-value = 1.46e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-09

TCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCAT  >  NC_000913/2434737‑2435030
                                                                                                                                                 |                                                                                                                                                     
tCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGATCATACCGGTCAGGGGAATcc                                                                                                                                                  >  1:196392/1‑151 (MQ=255)
         gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATCCACCGCGAcc                                                                                                                                         >  2:422242/1‑151 (MQ=255)
         gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGAc                                                                                                                                          <  1:230761/150‑1 (MQ=255)
           cTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCCGGGGAATCCACCGCGAcc                                                                                                                                         >  2:485980/1‑149 (MQ=255)
            tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCGCGAcc                                                                                                                                         >  1:88052/1‑148 (MQ=255)
             aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGAt                                                                                                                                     <  1:861727/151‑1 (MQ=255)
              tGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATAACCGGCCGGGGAATCCCCCGCGGcc                                                                                                                                         >  1:690904/1‑146 (MQ=255)
              tGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGAt                                                                                                                                     <  1:663339/150‑1 (MQ=255)
               gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATAc                                                                                                                                   <  1:950693/151‑1 (MQ=255)
                    gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGAGCATAACTGGCAGGGGAATCCACCGCGACCCGATACAGccc                                                                                                                              >  1:751922/1‑151 (MQ=255)
                        tCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTcc                                                                                                                         >  1:603951/1‑151 (MQ=255)
                                 ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCtttt                                                                                                                   <  1:132743/148‑1 (MQ=255)
                                  tCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGcc                                                                                                                <  1:651429/150‑1 (MQ=255)
                                    aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCCACCAGATACAACCCTTTCGCTTTTTCCgct                                                                                                             >  1:719906/1‑149 (MQ=255)
                                         gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATTATACCGGGCAGGGGAAACCCCCGCGACCAGATAAAGCCCTTCCGCTTTTGTCGtt                                                                                                             >  2:857487/1‑146 (MQ=255)
                                          tAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCCCGCCCCGATACAACCCCTCCGCCTTTGCCGTtgctg                                                                                                         >  2:618861/1‑149 (MQ=255)
                                          tAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGGCAGGGGAATCCACCGCCACCAGATACAGCCCTTCCCCTTTTTCCGTTgccgcc                                                                                                       >  2:672083/1‑151 (MQ=255)
                                           aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCCTCCGCTTTTGCCGGTGCTgcc                                                                                                       >  2:716822/1‑150 (MQ=255)
                                           aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTgcc                                                                                                       <  1:831405/150‑1 (MQ=255)
                                                                aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCCGGGGAAATCACCGCGACCAGATAAAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGATCTTTCCTTTgcc                                                                                  >  1:914526/1‑150 (MQ=255)
                                                                     ttCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCCCCAGCGTTCTTTCCTTTTCCGCCAg                                                                             >  1:567466/1‑150 (MQ=255)
                                                                      tCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCCGCGTTCTGTCCTTTTccgcccgc                                                                            >  1:224207/1‑150 (MQ=255)
                                                                             attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTATACCCGTCAGGGGAAACCACCCCGACCCGATACAGCCCCTCCGCTTTTGCCGTTGCTGCCCCCCGCGGTCTGTCCTTTGCCCCCCGCAACTCTg                                                                     >  2:299212/1‑150 (MQ=255)
                                                                             attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGCAGATCATACCCGTTAAGGGGAACCACCGCGCCCAGATCCACGCCCTCCGCTTTTTCCAGTGTTGGCGCCACCGCGCTGTCCTTTTCCGCCAGCGCCTCTGc                                                                    >  1:807280/1‑151 (MQ=255)
                                                                              ttattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATACCGGTTAGGGTAAACCCCCACGACCAAAAACAGCCCTTCCGCTTTTGCCGTTGCTGCCCCCCTCCTTCTGTCCTTTtccgcccccaa                                                                          >  1:456178/1‑139 (MQ=255)
                                                                              ttattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACCCTGc                                                                    >  1:261126/1‑150 (MQ=255)
                                                                                   ctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACCGTCAGGGGAAATCACCGCGACCAGATACAAGCCTTCCCCTTTTGCCGTTGCTGCCCCCAGCGGTCTGTCCCTTTCCCCCCGCAACTCTGCTCTcc                                                               >  2:322921/1‑150 (MQ=255)
                                                                                    tcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCGCGGCCAGATACAGCCCTTCCCCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTGCGCCCGCAACTCTGCt                                                                   >  2:474871/1‑145 (MQ=255)
                                                                                    tcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAAATCATACCCGGCAGGGGAATTCACCCCGACCCGATACAGCCCCTCCGCTTTTTCCGTTGCTCCCCCCAGCGCTCTGTCCTTTTCCCCCAACCACtcttcttttc                                                               >  2:707252/1‑145 (MQ=255)
                                                                                    tcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAAATAATACCCGGCCGGGGAATCCACCGCGGCCAGATACAGACCTTCCCCTTTTTCCGCTGCTGGCCCCAGCGGTTTGTCCTTTGCCGCCACCAACTCTGCTATcc                                                               >  1:867310/1‑149 (MQ=255)
                                                                                      gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCCGGCAGGGGAATCCACCGCGGCCAGATACAGCCCCTTCGCCTTTTCCGTTGCTGCCCCCCGCCGTCTGTCCTTTTCCGCCAGCAACCCTGCCATcc                                                               >  2:593574/1‑147 (MQ=255)
                                                                                         aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAACGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCCGCGTTCTGTCCTTTGCCGCCAGCAACCCTGCTATCCAGCcctc                                                        >  1:562139/1‑151 (MQ=255)
                                                                                           tCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGAAGCTTCCGGTTAGGGTAAGCCACCCCCGCCGGATAAAGGCCCTCGGCTTTTCCCCCTGCCCCCGCCCGCGGTCTGTCCCTTGCCGCCAGGCACACAGCCGTCCACCTCTcc                                                       >  2:116074/1‑150 (MQ=255)
                                                                                                  gAAATAGCGGGCCCATTGGCGGTTTTGGGAAGTTATCCCGGGGCGGGGAAAACACCCCGGCCAGCAACAACCCCTCCCCCTTTGCCGTTGCTGCCGCCCACGGTTCCTGCTTTGACGCCCGCCAAACCGCTATCCCGCTCGCCGGCTggg                                                >  2:665129/1‑149 (MQ=255)
                                                                                                     aTAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGCAATCCACCGCGACCAGATAAAACCCTTCCCCCTTTTCCCGTTCTGCCGCCAGCGGTCTGGCCTTTTCCCCCAACAACCCTTCTATTCAGCCCTCCTGCCgggtttg                                            >  2:790943/1‑146 (MQ=255)
                                                                                                     aTAGCGGGCCCATTGGCGGTTTTGGAAGATCATAACGGGCAGGGGAAACCACCGCGACCAGATAAAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCTTTGCCGCCAGCAACTCTGCTATTCAGCTCTCCGGCTgggt                                               >  1:619634/1‑146 (MQ=255)
                                                                                                         cGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCa                                        <  2:631115/151‑1 (MQ=255)
                                                                                                          gggCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATCCACCCCGACCCGATACAGCCCTTCCCCTTTTGCCGTTGCTGCCCCCCGCGGTCTGTCCTTTGCCGCCAGCAACTCTTCTATCCCGCTCTCCGGCTGGGTGTGGGCAc                                       >  2:378267/1‑151 (MQ=255)
                                                                                                          gggCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCCCTTTTGCCGTTGCTGCCGCCCGCGTTCTGTCCCTTTCCGCCCGCAACCCTTCTTTCCCGCTCTCCCGCCGGGTGTgggga                                        >  2:744815/1‑148 (MQ=255)
                                                                                                                   ggCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc                               >  1:673002/1‑150 (MQ=255)
                                                                                                                         tttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAgg                         <  1:64108/150‑1 (MQ=255)
                                                                                                                          ttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt                       >  1:626696/1‑151 (MQ=255)
                                                                                                                           tGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt                       <  1:367449/150‑1 (MQ=255)
                                                                                                                             gAAGATCATACCGGTCGGGGTAATCCCCCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg                      <  2:153352/149‑1 (MQ=255)
                                                                                                                             gAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGcc                    >  1:146369/1‑151 (MQ=255)
                                                                                                                               aGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCg                   >  2:452993/1‑150 (MQ=255)
                                                                                                                                    aTACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat              >  1:543803/1‑150 (MQ=255)
                                                                                                                                      aCCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAatat            <  1:726165/150‑1 (MQ=255)
                                                                                                                                      aCCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATAtt           <  2:786315/151‑1 (MQ=255)
                                                                                                                                          gTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTg       <  1:708237/151‑1 (MQ=255)
                                                                                                                                             aGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc    >  2:571984/1‑151 (MQ=255)
                                                                                                                                             aGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc    >  2:565690/1‑151 (MQ=255)
                                                                                                                                             aGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc    >  1:211994/1‑151 (MQ=255)
                                                                                                                                             aGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc    >  2:844074/1‑151 (MQ=255)
                                                                                                                                               ggTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCat  >  2:641528/1‑151 (MQ=255)
                                                                                                                                               ggTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCat  <  1:95989/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCAT  >  NC_000913/2434737‑2435030

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