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

TCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCAT  >  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)
                        tCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTT‑cc                                                                                                                         >  1:603951/1‑151 (MQ=255)
                                 ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCtttt                                                                                                                   <  1:132743/148‑1 (MQ=255)
                                  tCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGcc                                                                                                                <  1:651429/150‑1 (MQ=255)
                                    aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCCACCAGATACAACCCTT‑TCGCTTTTTCCgct                                                                                                             >  1:719906/1‑149 (MQ=255)
                                         gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATTATACCGGGCAGGGGAAACCCCCGCGACCAGATAAAGCCCTT‑CCGCTTTTGTCGtt                                                                                                             >  2:857487/1‑146 (MQ=255)
                                          tAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCCCGCCCCGATACAACCCCT‑CCGCCTTTGCCGTtgctg                                                                                                         >  2:618861/1‑149 (MQ=255)
                                          tAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGGCAGGGGAATCCACCGCCACCAGATACAGCCCTT‑CCCCTTTTTCCGTTgccgcc                                                                                                       >  2:672083/1‑151 (MQ=255)
                                           aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCCT‑CCGCTTTTGCCGGTGCTgcc                                                                                                       >  2:716822/1‑150 (MQ=255)
                                           aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTgcc                                                                                                       <  1:831405/150‑1 (MQ=255)
                                                                aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCCGGGGAAATCACCGCGACCAGATAAAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGATCTTTCCTTTgcc                                                                                  >  1:914526/1‑150 (MQ=255)
                                                                     ttCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCCCCAGCGTTCTTTCCTTTTCCGCCAg                                                                             >  1:567466/1‑150 (MQ=255)
                                                                      tCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCCGCGTTCTGTCCTTTTccgcccgc                                                                            >  1:224207/1‑150 (MQ=255)
                                                                             attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTATACCCGTCAGGGGAAACCACCCCGACCCGATACAGCCCCT‑CCGCTTTTGCCGTTGCTGCCCCCCGCGGTCTGTCCTTTGCCCCCCGCAACTCTg                                                                     >  2:299212/1‑150 (MQ=255)
                                                                             attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGCAGATCATACCCGTTAAGGGGAACCACCGCGCCCAGATCCACGCCCT‑CCGCTTTTTCCAGTGTTGGCGCCACCGCGCTGTCCTTTTCCGCCAGCGCCTCTGc                                                                    >  1:807280/1‑151 (MQ=255)
                                                                              ttattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATACCGGTTAGGGTAAACCCCCACGACCAAAAACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCCCCCTCCTTCTGTCCTTTtccgcccccaa                                                                          >  1:456178/1‑139 (MQ=255)
                                                                              ttattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACCCTGc                                                                    >  1:261126/1‑150 (MQ=255)
                                                                                   ctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACCGTCAGGGGAAATCACCGCGACCAGATACAAGCCTT‑CCCCTTTTGCCGTTGCTGCCCCCAGCGGTCTGTCCCTTTCCCCCCGCAACTCTGCTCTcc                                                               >  2:322921/1‑150 (MQ=255)
                                                                                    tcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCGCGGCCAGATACAGCCCTT‑CCCCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTGCGCCCGCAACTCTGCt                                                                   >  2:474871/1‑145 (MQ=255)
                                                                                    tcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAAATCATACCCGGCAGGGGAATTCACCCCGACCCGATACAGCCCCT‑CCGCTTTTTCCGTTGCTCCCCCCAGCGCTCTGTCCTTTTCCCCCAACCACtcttcttttc                                                               >  2:707252/1‑145 (MQ=255)
                                                                                    tcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAAATAATACCCGGCCGGGGAATCCACCGCGGCCAGATACAGACCTT‑CCCCTTTTTCCGCTGCTGGCCCCAGCGGTTTGTCCTTTGCCGCCACCAACTCTGCTATcc                                                               >  1:867310/1‑149 (MQ=255)
                                                                                      gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCCGGCAGGGGAATCCACCGCGGCCAGATACAGCCCCT‑TCGCCTTTTCCGTTGCTGCCCCCCGCCGTCTGTCCTTTTCCGCCAGCAACCCTGCCATcc                                                               >  2:593574/1‑147 (MQ=255)
                                                                                         aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAACGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTTCCGTTGCTGCCGCCCGCGTTCTGTCCTTTGCCGCCAGCAACCCTGCTATCCAGCcctc                                                        >  1:562139/1‑151 (MQ=255)
                                                                                           tCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGAAGCTTCCGGTTAGGGTAAGCCACCCCCGCCGGATAAAGGCCCT‑CGGCTTTTCCCCCTGCCCCCGCCCGCGGTCTGTCCCTTGCCGCCAGGCACACAGCCGTCCACCTCTcc                                                       >  2:116074/1‑150 (MQ=255)
                                                                                                  gAAATAGCGGGCCCATTGGCGGTTTTGGGAAGTTATCCCGGGGCGGGGAAAACACCCCGGCCAGCAACAACCCCT‑CCCCCTTTGCCGTTGCTGCCGCCCACGGTTCCTGCTTTGACGCCCGCCAAACCGCTATCCCGCTCGCCGGCTggg                                                >  2:665129/1‑149 (MQ=255)
                                                                                                     aTAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGCAATCCACCGCGACCAGATAAAACCCTTCCCCCTTTTCCCGTT‑CTGCCGCCAGCGGTCTGGCCTTTTCCCCCAACAACCCTTCTATTCAGCCCTCCTGCCgggtttg                                            >  2:790943/1‑146 (MQ=255)
                                                                                                     aTAGCGGGCCCATTGGCGGTTTTGGAAGATCATAACGGGCAGGGGAAACCACCGCGACCAGATAAAGCCCTT‑CCGCTTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCTTTGCCGCCAGCAACTCTGCTATTCAGCTCTCCGGCTgggt                                               >  1:619634/1‑146 (MQ=255)
                                                                                                         cGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCa                                        <  2:631115/151‑1 (MQ=255)
                                                                                                          gggCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATCCACCCCGACCCGATACAGCCCTT‑CCCCTTTTGCCGTTGCTGCCCCCCGCGGTCTGTCCTTTGCCGCCAGCAACTCTTCTATCCCGCTCTCCGGCTGGGTGTGGGCAc                                       >  2:378267/1‑151 (MQ=255)
                                                                                                          gggCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTT‑CCCCTTTTGCCGTTGCTGCCGCCCGCGTTCTGTCCCTTTCCGCCCGCAACCCTTCTTTCCCGCTCTCCCGCCGGGTGTgggga                                        >  2:744815/1‑148 (MQ=255)
                                                                                                                   ggCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc                               >  1:673002/1‑150 (MQ=255)
                                                                                                                         tttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAgg                         <  1:64108/150‑1 (MQ=255)
                                                                                                                          ttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt                       >  1:626696/1‑151 (MQ=255)
                                                                                                                           tGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt                       <  1:367449/150‑1 (MQ=255)
                                                                                                                             gAAGATCATACCGGTCGGGGTAATCCCCCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg                      <  2:153352/149‑1 (MQ=255)
                                                                                                                             gAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGcc                    >  1:146369/1‑151 (MQ=255)
                                                                                                                               aGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCg                   >  2:452993/1‑150 (MQ=255)
                                                                                                                                    aTACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat              >  1:543803/1‑150 (MQ=255)
                                                                                                                                      aCCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAatat            <  1:726165/150‑1 (MQ=255)
                                                                                                                                      aCCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATAtt           <  2:786315/151‑1 (MQ=255)
                                                                                                                                          gTCAGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTg       <  1:708237/151‑1 (MQ=255)
                                                                                                                                             aGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc    >  2:571984/1‑151 (MQ=255)
                                                                                                                                             aGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc    >  2:565690/1‑151 (MQ=255)
                                                                                                                                             aGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc    >  1:211994/1‑151 (MQ=255)
                                                                                                                                             aGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc    >  2:844074/1‑151 (MQ=255)
                                                                                                                                               ggTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCat  >  2:641528/1‑151 (MQ=255)
                                                                                                                                               ggTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCat  <  1:95989/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTT‑CCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCAT  >  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: ‑