Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,903,6920CA44.1% 78.2 / 45.6 59intergenic (+318/+1055)ygcE/queEputative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base C (4/29);  new base A (26/0);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 7.84e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGA  >  NC_000913/2903542‑2903842
                                                                                                                                                      |                                                                                                                                                        
gTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGTTTATGTAAAAc                                                                                                                                                          >  1:93485/1‑151 (MQ=255)
      aCAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACa                                                                                                                                                        >  2:475626/1‑147 (MQ=255)
                 tAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAAAACTCATTATCCCCaa                                                                                                                                          >  1:92055/1‑150 (MQ=255)
                 tAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACTCATTAGCCTCaa                                                                                                                                          >  2:37080/1‑150 (MQ=255)
                  aaCTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCTCTCATTAGCCTCaaaa                                                                                                                                        <  2:181819/151‑1 (MQ=255)
                    cTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCCATTTACCTCcaaa                                                                                                                                        >  2:537646/1‑149 (MQ=255)
                      acacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACa                                                                                                                                      <  1:232666/149‑1 (MQ=255)
                       cacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCATTAGCCTCaaaaaaaa                                                                                                                                    >  1:270380/1‑150 (MQ=255)
                           cGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACaa                                                                                                                              <  2:455293/151‑1 (MQ=255)
                            gATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACCCTCATTAACCTCcaaacaaaaaaaa                                                                                                                               >  2:481364/1‑145 (MQ=255)
                            gATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACCCATTTACCTCAAAACaaaaaaa                                                                                                                                >  2:177852/1‑145 (MQ=255)
                            gATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATCTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGGTTTTTTAAAACAAACATTATCCTCaaaaaaaaaaaaa                                                                                                                               >  2:320513/1‑145 (MQ=255)
                                   ttAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCa                                                                                                                      <  1:918607/150‑1 (MQ=255)
                                    tAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACTCCATAACCTCCAaaaaaaaaaaaa                                                                                                                             >  2:887771/1‑143 (MQ=255)
                                       aaTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACTCATTTAACTCCAAACAAAAAAAAAACTCCAAAttt                                                                                                                 >  1:206060/1‑151 (MQ=255)
                                         tAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAAAACTCATTATCCTCAAAAAAAAAAAAAAAATCCAGAatttt                                                                                                               >  1:145676/1‑149 (MQ=255)
                                                cAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACTACTTATGCTTAAAAAAAAAAAAAAAATACAGATttattta                                                                                                            >  1:808252/1‑145 (MQ=255)
                                                cAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACAAACCTTAAACTCAAAAAAAAAACAAAAATACAAATttatttat                                                                                                           >  2:1025686/1‑145 (MQ=255)
                                                       aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGa                                                                                                 <  1:515934/151‑1 (MQ=255)
                                                       aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGa                                                                                                 <  2:354528/151‑1 (MQ=255)
                                                          cGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATa                                                                                              <  2:608538/151‑1 (MQ=255)
                                                           gCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATaa                                                                                             <  2:727435/151‑1 (MQ=255)
                                                              cTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACTCATTTTCCCCAAAAAAAAAAAAAAAATCCAAATTTATTTCTTTAACAATGAATaaaa                                                                                           >  2:749105/1‑148 (MQ=255)
                                                               tttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAACACTCATTAACCCCAAAAAAAAAAAAAAATTCCAGATTTTTTTCTTAAACAATTaaaaaaaa                                                                                          >  2:925530/1‑150 (MQ=255)
                                                                ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACaa                                                                                          <  2:652235/149‑1 (MQ=255)
                                                                ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACaa                                                                                          <  2:175344/149‑1 (MQ=255)
                                                                 tATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGc                                                                                        <  2:423025/150‑1 (MQ=255)
                                                                    cAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTaa                                                                                    <  1:120047/151‑1 (MQ=255)
                                                                       ttttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCATTAACCTTCAAAAAAAAAAAAAATTCCCCATTTATTTATTTACCAAAAAAAAACaaaatttaaa                                                                                  >  2:258932/1‑142 (MQ=255)
                                                                         ttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCCTTATCCTCAAAAAAAAAAAAAACATCCAGATTTTTTAATTTAACAAAGAAAAAAAAAAATaaaaaa                                                                                >  1:761076/1‑150 (MQ=255)
                                                                          tCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCCATTAGCCCCCAAAAAAAAAAAAAATTCCAGATTTTTTTCTTTAACAATGAATAAAAAAATTtaaaaa                                                                                >  2:20803/1‑149 (MQ=255)
                                                                               tAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCCATTATCCTCCAAAAAAAAACAAAAATACCGATTTATTAATTAAACAAAAAATAACAAACTTAAAAAAtataat                                                                          >  2:443902/1‑148 (MQ=255)
                                                                               tAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACt                                                                          <  1:104008/150‑1 (MQ=255)
                                                                                     aTGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTc                                                                   <  1:414483/151‑1 (MQ=255)
                                                                                        aaCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTaa                                                                <  2:295278/151‑1 (MQ=255)
                                                                                         aCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTGTTTTAAAACACTCATTATGCTCAAAAAAAAAAAAAACTTCCAAGTTTTTTTACTTACCAAAGAATAACAAACTTAAAAAAAATACTATttttttta                                                                >  2:676357/1‑146 (MQ=255)
                                                                                          cTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACCCATTAACCCCAAAAAAAAAAAAAAAATTCAGATTTATTTATTTAACAAAGAAAAAAAAAATTAAAAAAaataattttttttttaa                                                               >  1:650704/1‑137 (MQ=255)
                                                                                           tGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACTCATTAACCTCAAAAAAAAAAAAAAATTCCAAAATTTTTATTTAACAAAAAAAAAAAAAACTTAAAAAAAAAAATattatttttaaa                                                              >  1:123515/1‑144 (MQ=255)
                                                                                                   tttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAACAACTCATTATCCTCAAAAAAAAAACAAAAATCCAGATTTATTTCTTTAACAAAGAAAAAAAAAATAAAAAAAAATATTTTTTTTTTaaaaaaaaaa                                                        >  1:729419/1‑148 (MQ=255)
                                                                                                    ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACTCCTTAGCCTCCAAAAAAAAACAAACTTCCCGATTTATTACTTAACCAATAAAAAAAAAACTTAAAAAAAATACTTTTTTTTTTAaaaaaaaa                                                        >  2:760309/1‑147 (MQ=255)
                                                                                                    ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGGTTTTTAAAAACACCCATTTTCCTCCAAAAAAAAAAAAAAATCCCGAATTATTATTTTACCAAAAAAAAAAAAAATTAAAAAAAAAAATATTTTTTTTAaaaaaaaaa                                                       >  1:871076/1‑147 (MQ=255)
                                                                                                           tAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCCATTAGCCTCAAAAAAAAAAAAAAAATTCAGATTTATTTCTTTAACAAAAAATAAAAAAATTAAAAAAAAAAATGTTATTTTAAAAAAAAAAATaaaa                                                  >  1:1020563/1‑144 (MQ=255)
                                                                                                             aGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTa                                           <  2:581422/151‑1 (MQ=255)
                                                                                                              gAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTaa                                          <  2:1020665/151‑1 (MQ=255)
                                                                                                              gAACTATTCATCAACAGCGTATAGAGGCGGCTAATTAAAAACACCCCTTATACTCAAAAAAAAAACAAAACTTCCGCCTTTTTTCTTTACCCAATAAAAAAAAACATAAAAAATATACATTTATTTTTAAGAAAAAAATAACaa                                                 >  2:570050/1‑143 (MQ=255)
                                                                                                                    ttCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTc                                       <  1:574943/148‑1 (MQ=255)
                                                                                                                      catcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAata                                   <  1:481324/150‑1 (MQ=255)
                                                                                                                      catcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAata                                   <  2:131390/150‑1 (MQ=255)
                                                                                                                       atcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATAtt                                 <  1:567533/151‑1 (MQ=255)
                                                                                                                         caACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTa                                <  1:202825/150‑1 (MQ=255)
                                                                                                                            cAGCGTATAGAGGCGGTTATGTAAAAACACTCATTAGCCCCAAAAAAAAAACAAACATCCAGATTTATTACTTAAACAAAGAAAAAAAAACATAAAAAAAAAAATGTTTTTTTAAAAAAAAAAATAAAAAGCATTAatttaaaatttaaa                             >  1:673488/1‑150 (MQ=255)
                                                                                                                             aGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaa                            <  1:364326/150‑1 (MQ=255)
                                                                                                                                    agagGCGTTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAg                    <  2:93494/151‑1 (MQ=255)
                                                                                                                                    agagGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAg                    <  1:320487/151‑1 (MQ=255)
                                                                                                                                         cGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTaa               <  2:650767/151‑1 (MQ=255)
                                                                                                                                         cGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTa                >  2:337387/1‑150 (MQ=255)
                                                                                                                                             tATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAt            <  2:688292/150‑1 (MQ=255)
                                                                                                                                              aTGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAttt          <  1:475586/151‑1 (MQ=255)
                                                                                                                                               tGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAGATAATAGATTAAAAttt          >  1:290490/1‑150 (MQ=255)
                                                                                                                                                 tAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCt        <  1:887688/150‑1 (MQ=255)
                                                                                                                                                      ccACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGa  <  1:104358/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                        
GTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGA  >  NC_000913/2903542‑2903842

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