Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 816,752 T→A 11.6% intergenic (‑105/‑292) ybhK ← / → moaA putative CofD superfamily transferase/molybdopterin biosynthesis protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913816,7520TA11.6% 276.3 / 4.5 86intergenic (‑105/‑292)ybhK/moaAputative CofD superfamily transferase/molybdopterin biosynthesis protein A
Reads supporting (aligned to +/- strand):  ref base T (51/24);  new base A (2/9);  total (53/33)
Fisher's exact test for biased strand distribution p-value = 2.38e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAG‑AG‑AACG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAAT  >  NC_000913/816517‑816976
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                
ttATCGGTGGTTGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAG‑AG‑AACG‑CGTCCCAGGCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTaaa                                                                                                                                                                                                                               <  1:1120287/239‑1 (MQ=255)
  ttCGGTGGGTGTGAAGATCCCCGTAAAACGCAACCCCAAAGACGAAAG‑GTAGAGAAG‑AACG‑CCCCCCCGCCCGTCCCCTCCCCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAgg                                                                                                                                                                                                                             <  2:714384/238‑1 (MQ=255)
   tCGGTGGTGGTGACGATACCCGTTAAACGAACACCCAAAGGCGAAAGTGATGAG‑AG‑AACG‑CGTCCCAGTCCATGCCCCCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAACCAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGa                                                                                                                                                                                                                            <  1:646442/239‑1 (MQ=255)
   gcGGTGGTGGTGACGATCCCCTTTAAACGAAACCCCAAAGAGGAAAGGGAGGAG‑AG‑AACG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGa                                                                                                                                                                                                                            <  2:118053/238‑1 (MQ=255)
   gcGGTGGGGGTAACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAG‑AG‑AACG‑CGGCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAAATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGa                                                                                                                                                                                                                            <  1:452799/238‑1 (MQ=255)
       tggtggtgACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAG‑AG‑AACG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTTCCAATATTGAGGATACAAACTATTTTATTCTACCTCAAAAAGACGaa                                                                                                                                                                                                                        >  2:1305099/1‑238 (MQ=255)
           ggtgACGATCCCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAG‑AG‑AACG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCa                                                                                                                                                                                                                     <  1:264690/238‑1 (MQ=255)
            gtgACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAG‑AG‑AACG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCAc                                                                                                                                                                                                                    >  1:118053/1‑238 (MQ=255)
               aCGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAG‑AG‑AACG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATA‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGCGGATACAAACTATTTTATTTTACCTCAAAAGAAAGAAGAACGcag                                                                                                                                                                                                                >  2:716824/1‑237 (MQ=255)
                cGATACCCGTTAAACGAGAACCCAAAGCCGAAAGGGATGAG‑AG‑AACG‑CGTCCCAGCCCAGGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTa                                                                                                                                                                                                                <  1:917270/238‑1 (MQ=255)
                 gATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAG‑AG‑AACG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGTTACAAACCATTTTTTTCTACCTTTAAAGGACAATgaaaccaat                                                                                                                                                                                                               >  2:1202001/1‑230 (MQ=255)
                 gATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAG‑AG‑AACG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGGACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCAACCTCTAAAGGACGATGCAAGCTAt                                                                                                                                                                                                               >  2:134557/1‑238 (MQ=255)
                       ccttAAACGGGACCCCAAAGACGAAAGTGA‑GAG‑GGAAACG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGACTccc                                                                                                                                                                                                        <  2:1137345/237‑1 (MQ=255)
                         ttAAACGAGAACCCAAAGACGAAAGTGATGAG‑AG‑AACG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCtg                                                                                                                                                                                                      >  1:305465/1‑239 (MQ=255)
                         ttAAACGAGAACCCAAAGACGAAAGTGATGAG‑AG‑AACG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGAAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGAAACAAACAATTTTCTTCTACCTCTAAAGAACGATGAACGCTATGCCTCCCtg                                                                                                                                                                                                      >  2:1273357/1‑239 (MQ=255)
                           aaaCGAGAACCCAAAGACGAAAGTGATGAG‑CG‑AACG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCATCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGGGGATACAAACTATTTTTTTCTACCTCTAAAggc                                                                                                                                                                                                                            >  1:1352414/1‑214 (MQ=255)
                            aaCGGGCACCCAAAGGCGAAAGTGATGAG‑A‑‑AACGCCGTCCCAGCCCAGGCCCCGCGCCGAGAG‑AAACGACCCGAAC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATGATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCtgatg                                                                                                                                                                                                   <  2:290560/239‑1 (MQ=255)
                            aaCCAGAACCCAAAGACGAAAGGGATGAG‑AA‑CACC‑GCTCCCAGTCCAGGCCCTCCGCCG‑GAGCCAACGAAACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCtgatg                                                                                                                                                                                                   <  2:58926/239‑1 (MQ=255)
                              agagaACCACAAGCAGAACGTGATGAG‑AA‑AACG‑CGCCCCAGTCAATGCACTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTgatgat                                                                                                                                                                                                 <  1:1012652/238‑1 (MQ=255)
                                  aaCCCAAAGACGAAAGTGATGAG‑AG‑AACG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGACCGCTATGCCTCCCTGATGATGtat                                                                                                                                                                                             >  1:335042/1‑239 (MQ=255)
                                       aaaGACGAAAGTGATGAG‑AG‑AACG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCAGCTGATGTATATCaa                                                                                                                                                                                        >  1:807001/1‑239 (MQ=255)
                                               aaGTGATGAG‑AG‑AACG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTaaaa                                                                                                                                                                                 >  1:731026/1‑238 (MQ=255)
                                                              aCG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATTAATATCAAAGAAAAACCGCCATATTCCCt                                                                                                                                                                   >  2:18570/1‑239 (MQ=255)
                                                                   gtccaacgcaaTGCCCTCCGCCGAGA‑‑CAACGCACCGGTCAAAGATCATCCAGCAAAGGATGCCTCATATAGTTCCCGGGAGTCAGATTACCCGCTATACAGGAGCACAATTTG‑GGG‑GAAAAAGAAATTAACCCTGCAAATATGGAGGAAACAAATTATTTTCTTCTACCTCAAAGGGACGAGGCAGGCAATCCCTCCGTGATGGTGTATATAAAAGTAAAACCGCCATTTTCCCTTAtt                                                                                                                                                               <  1:434395/229‑1 (MQ=255)
                                                                    tCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCATTATTGCGGATCCAAACTATTTTCTCCTACATCTAAAGGACGATGCACGCTATGCACACCTAATGATGTATATCAAATTAAAACCCCCCATTTCCCTTAttt                                                                                                                                                              >  2:434395/1‑238 (MQ=255)
                                                                      ccAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCGCATGATGTAAATCAAAGTAAAACACCCATTTTCCCTTAATCTg                                                                                                                                                            >  2:673936/1‑239 (MQ=255)
                                                                           ccATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCG‑CATTTTCCCTTATTCTGTAg                                                                                                                                                         <  1:1305099/236‑1 (MQ=255)
                                                                               gCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATg                                                                                                                                                                                                              >  1:378921/1‑180 (MQ=255)
                                                                               gCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATg                                                                                                                                                                                                              <  2:378921/180‑1 (MQ=255)
                                                                               gCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCCTTTTCCCTTTTTTTTTAGCGaaaa                                                                                                                                                   >  2:391070/1‑238 (MQ=255)
                                                                                 ccaccgccgAGAG‑AAACGTCAGGATA‑AATACCAGCCAGCGTCGGTTTGCCCATATATTTCCCTGGAGCCAGATTACCCGCGCTACAGAAGCGCAATTTGCGGG‑AAAACAGCAATAAA‑CCTCCAAATATGGAGGATACAACTTATTTTCTTCTACCTCAAAAGGACGAGGCACGTTTTGCCTCCCTGATGATGTATATAAAAGAAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAg                                                                                                                                                 <  1:18570/236‑1 (MQ=255)
                                                                                    ccgccgAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCACTATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGAACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc                                                                                                                                              >  1:61445/1‑239 (MQ=255)
                                                                                    ccgccgAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc                                                                                                                                              >  1:157806/1‑239 (MQ=255)
                                                                                    ccgccgAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc                                                                                                                                              >  1:90305/1‑239 (MQ=255)
                                                                                    ccgccgAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc                                                                                                                                              >  1:1070663/1‑239 (MQ=255)
                                                                                    ccgccgAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc                                                                                                                                              >  1:1101244/1‑239 (MQ=255)
                                                                                    ccgccgAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc                                                                                                                                              >  1:422662/1‑239 (MQ=255)
                                                                                    ccgccgAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTAATTCTTTAGCGTAATAGCAc                                                                                                                                              >  1:372002/1‑239 (MQ=255)
                                                                                            gaaaaaCCACACGAAA‑AAGATCAGCCAGCGTACGATTGCTCATAGAGTTGCCGGGAGTCAAATTACCCGCGCAACGGTAGCGAAAATTGGGGG‑AAAACAGAAATTAC‑CCTGCAAATATTGAGGATACAAACTTTTTTCTTCTACCTCTAAAGGACGATGCACGTTATGCCTCCCTGATGATGTATATCAAAGAAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCAt                                                                                                                                        <  1:673936/234‑1 (MQ=255)
                                                                                              caaaCGACAAGCAC‑AAGCTCACCAAGCGTACTATTGCGCATAGTGTTTACGGGAGTCAGATTACCCGCGACACAGTAGCCCAAATTGGGGG‑AAAAAAGAAATTAC‑CCTGCCAATATTGAGGATACAAACTTTTTTCTTCTACCTCTAAAGGACGTTGCACGCTAGGCCTCCCTGAGGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGa                                                                                                                                      <  1:716824/236‑1 (MQ=255)
                                                                                               cAAAGACACGAAA‑AAGATCGGCCAGCGAACAATTGCCAAAATATTTTCCGGGAGCAAGATTACCCGGGCTAAAGTAGCGAAATTTGGGGG‑AAAACAGAAATAAA‑CCTGCCAATATTGAGGTTACAAACTTTTTTTTTCTCCCTCAAAAGGACGATGCACGTTATGCCTCCCTGATGATGTATATAAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACg                                                                                                                                    <  1:1202001/239‑1 (MQ=255)
                                                                                                  cGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCta                                                                                                                                 >  1:762644/1‑239 (MQ=255)
                                                                                                  cGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACAATCATGACCCta                                                                                                                                 >  2:802798/1‑239 (MQ=255)
                                                                                                       cGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGAAAAACCGCCATTTTCCCTTATTCTGTAGCGAACAAGCCACATCATGACACTATATAc                                                                                                                            >  2:86458/1‑239 (MQ=255)
                                                                                                                   agccagcGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATa                                                                                                                 >  1:1221095/1‑239 (MQ=255)
                                                                                                                     ccagcGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GGAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAAAAGCACAATCATGACCCAATAAAAATGATAACAAAgg                                                                                                               >  2:552114/1‑238 (MQ=255)
                                                                                                                     ccagcGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAAAAGAACGACCATGAAGCAATATACAAAATTaaaaaga                                                                                                               >  2:487154/1‑233 (MQ=255)
                                                                                                                         cGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGcaa                                                                                                           >  1:354431/1‑239 (MQ=255)
                                                                                                                                             ttCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACg                                                                                                                                                                                                                   <  1:1291433/115‑1 (MQ=255)
                                                                                                                                             ttCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACg                                                                                                                                                                                                                   >  2:1291433/1‑115 (MQ=255)
                                                                                                                                              tCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCa                                                                                      >  1:1040948/1‑239 (MQ=255)
                                                                                                                                               ccgggatccAGTTTATCCGGGCTAAAGTAGCAAAAATTGCGGG‑GAAACAGAAATTAA‑CCGGCAAATTGTGGGGATAAAACCTATTTCTTTCTACCTCAAAAGGGCGGTGCAAGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCAt                                                                                     <  1:309411/231‑1 (MQ=255)
                                                                                                                                                  ggttgcattttACCCGCGTTACAGTGGCAAAATTTGGGGG‑AAAACAGAAATTAC‑CCTGCCAATATTGATGATAAAAATTTTTTTCTTCTCCCTCTAAAGGACGATGCACGTTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTa                                                                                  <  2:1040948/230‑1 (MQ=255)
                                                                                                                                                   gAGTCAGATTACCCGCGCTACAGTAGCGAAATTTGGGGG‑AAAACAGAATTGAC‑CCTGCAAATATGGAGGATACAAACTTTTTTTTTCTACCTCAAAAGGACGATGCACGCTATGCCTCCCGGATGATGTATATCAAAGAAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTa                                                                                  <  1:1140507/238‑1 (MQ=255)
                                                                                                                                                      tCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGAACATGACGCTATTTACATGATTAAATAGCGAAAGTGGGGAAGGGAAAAATACCATTTaaac                                                                               >  2:449659/1‑238 (MQ=255)
                                                                                                                                                       gaatttccccGCTCTAGAGTA‑CGAAAATTGGGGGAAAAACAAAAATTAC‑CCGGCAAATTTGGGGGAACAAAACTTTTTTCTTCTACCACAAAAGGAAGCTCCCCTCTATCCCCCCCTGATGATGTAAAAAAAAATAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGc                                                                             <  2:157806/231‑1 (MQ=255)
                                                                                                                                                        aGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGCAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGGGTGGATGGTAAAAATCTCATTTACACGCt                                                                            >  1:1112611/1‑239 (MQ=255)
                                                                                                                                                        aGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTCCATAGCGAAAGGGTGGATGGGAAAAATCTCATTTACACGCt                                                                            >  1:305793/1‑239 (MQ=255)
                                                                                                                                                          aTTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTa                                                                                                                     >  1:696649/1‑196 (MQ=255)
                                                                                                                                                                      aCAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACAt                                                                                                                          >  1:709427/1‑179 (MQ=255)
                                                                                                                                                                        aGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAACTCTCATTTACACGCTAGTATCGGCATAACCa                                                            >  1:661419/1‑239 (MQ=255)
                                                                                                                                                                           agccaaaaTTGCGGA‑AAAACAAAATTTAA‑CCCACCAAGTTGGGGGATAAAATCTATTTTCTCCAACCTCAAAAGGCAGACGCCCGTGATCCCTCCCTTATGATGTATAAAAAATTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACt                                                          <  1:86458/233‑1 (MQ=255)
                                                                                                                                                                             cgcAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAGAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATCGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGCTGGTAAAAATCTCCTTTACACGCTAGAATCGGCATAACCACTaaa                                                       >  1:1291944/1‑239 (MQ=255)
                                                                                                                                                                             cgcAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTcc                                                                                                                                                                                                         <  1:160945/93‑1 (MQ=255)
                                                                                                                                                                             cgcAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTcc                                                                                                                                                                                                         >  2:160945/1‑93 (MQ=255)
                                                                                                                                                                             cgcAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTatat                                                                                                                                                                                           >  2:484619/1‑107 (MQ=255)
                                                                                                                                                                             cgcAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTatat                                                                                                                                                                                           <  1:484619/107‑1 (MQ=255)
                                                                                                                                                                             cgcAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCtatata                                                                                                                             >  1:1193559/1‑169 (MQ=255)
                                                                                                                                                                             cgcAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACAt                                                                                                                          >  1:1191696/1‑172 (MQ=255)
                                                                                                                                                                             cgcAAAGTGGGGG‑AAAACGACAATTAA‑CCGCCAAATTTAGGAGATCAAAATTATTTTCTTCTACCTTAAAAGGACGATGCACGTTATGCCTCCTTGATGAGGTATAACAAAGAAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTaa                                                        <  1:134557/238‑1 (MQ=255)
                                                                                                                                                                               cAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGAACGATCATGACGCTATATACATTATTACATAGCGAACGCGTGGCTTGTAAAACTCTCGTTTACACGAAAGAATCGGCAGAACCACAAAaca                                                     >  2:1025781/1‑239 (MQ=255)
                                                                                                                                                                                      cGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGGGTGGATGGTAAAAATCTCATTTACACGCTAGAATCGGCATAACCACTAAAAACCCTaa                                               >  2:1071928/1‑237 (MQ=255)
                                                                                                                                                                                           aaaaaAGAAATTAC‑CCTGCCAATATGGAGGATACAAACTATTTCCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCctc                                           <  1:1071928/237‑1 (MQ=255)
                                                                                                                                                                                               cAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCtga                                                                                                                                                                                                     >  2:935118/1‑80 (MQ=255)
                                                                                                                                                                                               cAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCtga                                                                                                                                                                                                     <  1:935118/80‑1 (MQ=255)
                                                                                                                                                                                                    tttAC‑CCGGCAAATATTGAGAAAAAAAACTTTTTCCTTCTCCCCCTAAGGGACGATGCCCGTTATGCCTCCTTGATGATGTATACCAAAGTAAATCCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTggg                                 <  1:578142/238‑1 (MQ=255)
                                                                                                                                                                                                             gCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACAt                                                                                                                          >  1:1301237/1‑142 (MQ=255)
                                                                                                                                                                                                             gCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCCTAACCAATAAAAACCCTATCCTCTCAACCCGGGGCAACGGa                         >  2:160198/1‑239 (MQ=255)
                                                                                                                                                                                                               cAATATTGGGGATAAAACTTATTTTCTTCTACCTCAAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATAAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATa                       <  1:449659/239‑1 (MQ=255)
                                                                                                                                                                                                                        ggATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCa                                                                                                                                                                                                                     >  2:560289/1‑40 (MQ=255)
                                                                                                                                                                                                                        ggATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCa                                                                                                                                                                                                                     <  1:560289/40‑1 (MQ=255)
                                                                                                                                                                                                                            aCAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCtga                                                                                                                                                                                                     >  2:1035672/1‑52 (MQ=255)
                                                                                                                                                                                                                            aCAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCtga                                                                                                                                                                                                     <  1:1035672/52‑1 (MQ=255)
                                                                                                                                                                                                                            aCAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGc                                              >  1:846077/1‑203 (MQ=255)
                                                                                                                                                                                                                            aCAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTTGcc          >  1:659264/1‑239 (MQ=255)
                                                                                                                                                                                                                                     tttCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAAt  >  1:105222/1‑238 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                
TTATCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAG‑AG‑AACG‑CGTCCCAGTCCATGCCCTCCGCCGAGAG‑CAACGACACGATC‑AAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGG‑GAAACAGCAATTAA‑CCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAAT  >  NC_000913/816517‑816976

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