Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,113,548 T→C 9.9% intergenic (+142/‑31) opgH → / → yceK OPG biosynthetic ACP‑dependent transmembrane UDP‑glucose beta‑1,2 glycosyltransferase; nutrient‑dependent cell size regulator, FtsZ assembly antagonist/outer membrane integrity lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,113,5480TC 9.9% 239.6 / 14.8 81intergenic (+142/‑31)opgH/yceKOPG biosynthetic ACP‑dependent transmembrane UDP‑glucose beta‑1,2 glycosyltransferase; nutrient‑dependent cell size regulator, FtsZ assembly antagonist/outer membrane integrity lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (37/35);  new base C (4/4);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.66e-03

TCTGCATTTCCGTGTCTGGAATTCCCCGGAGAGATATTCTTCATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATTAACGTT  >  NC_000913/1113292‑1113814
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                          
tCTGCATTTCCGTGTCTGGAATTCCCCGGAGAGATATTCTTCATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTTACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGa                                                                                                                                                                                                                                                          >  4:69581/1‑275 (MQ=255)
 ttGCATTTCTGGTGCTGGAATTCCCCGGTGAGTTATACTTCATGGGTGAGTTATGACGCAGGGATAAGGCTCAATCCAGTGGCTTTGCGGAAACGGGGTGCGGCTTCGCAATAAAAACGTAGTTTCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAAGTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGa                                                                                                                                                                                                                                                          <  4:408002/273‑1 (MQ=38)
       ttCCGTGTCTGGAATTCCCCGGAGAGATATTCTTCATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGt                                                                                                                                                                                                                                                   >  2:217521/1‑275 (MQ=255)
              tCTGGAATTCCCCGGAGAGATATTCTTCATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGa                                                                                                                                                                                                                                                          >  4:161404/1‑261 (MQ=255)
              tCTGGAATTCCCCGGAGAGATATTCTTCATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAAt                                                                                                                                                                                                                                            >  2:32091/1‑275 (MQ=255)
               cTGGAATTCCCCGGAGGGTTATTCTTAATGGTTGAGGTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAAt                                                                                                                                                                                                                                            <  4:391468/274‑1 (MQ=255)
                  gAATTCCCCGGAGAGATATTCTTCATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAaga                                                                                                                                                                                                                                                      <  2:161405/261‑1 (MQ=255)
                        cccGGAGAGATATTCTTCATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTCCGCTTATCAGGCCTACATCGTCCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTAAGAAAGAGGTATCAAATGCGTTTAAtt                                                                                                                                                                                                                                  >  2:155643/1‑275 (MQ=255)
                                  tatTCTTCATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGCTAAGGCGTTCACTCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTAGGCCGGTTTTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGcc                                                                                                                                                                                                                        >  4:49493/1‑274 (MQ=255)
                                        tcATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCTTTTAATTGTGGTGAGCATCATgg                                                                                                                                                                                                                  <  2:42046/275‑1 (MQ=255)
                                        acATGGGTTGTTTTTTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCTCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATgg                                                                                                                                                                                                                  <  4:405603/274‑1 (MQ=255)
                                         cATGGGGGAGTTATTACGAAGGGATAAAACTCAATCCACTGGGATTTCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGATACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATgg                                                                                                                                                                                                                  <  4:217520/274‑1 (MQ=255)
                                              gTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGAcc                                                                                                                                                                                                             <  4:367355/274‑1 (MQ=255)
                                              gTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCt                                                                                                                                                                                                            >  2:25319/1‑275 (MQ=255)
                                                       taCGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTg                                                                                                                                                                                                   >  2:278799/1‑275 (MQ=255)
                                                       agccaaGGGTTAAACCCCAACCCAATGGGATTGCGTAAACCGGATGCGCCTTCGCAATAAAAACTTAGTTGCCTGTTGCGCTACGCTTATCAGGCCGACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGt                                                                                                                                                                                                    <  4:380305/270‑1 (MQ=39)
                                                           aaGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCtg                                                                                                                                                                                               >  2:521219/1‑275 (MQ=255)
                                                                        cAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCACATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgcagctgttct                                                                                                                                                                                  >  2:380306/1‑269 (MQ=255)
                                                                        cAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCattatt                                                                                                                                                                                  >  2:314505/1‑275 (MQ=255)
                                                                         aaTCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCAttatta                                                                                                                                                                                 >  4:54290/1‑275 (MQ=255)
                                                                            ccgcgaggatgtcGTAAACCGTATGCGGCTTCGCAATAAAACCGTAGTTGCCTGATGCTCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTAGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAAGTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAg                                                                                                                                                                                <  4:522509/261‑1 (MQ=255)
                                                                            ccACTGGCATTGCGTAAACCGGATGCGGCTTCGAAATAAAAACGTAGTTGCCTGTTGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTCTTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAg                                                                                                                                                                                <  2:353864/273‑1 (MQ=255)
                                                                                     ttGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTATTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCACTTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAACGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCCTGGTGACCTTGCTGAGTGGCTGTGGCTGCATTCTTAGCCTCACTATAc                                                                                                                                                                     >  4:4093/1‑275 (MQ=255)
                                                                                     ttGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCCATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCGGTGGCAGCATTATTAGTCGCACtaacc                                                                                                                                                                     >  4:301104/1‑272 (MQ=255)
                                                                                        cGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGa                                                                                                                                                                                                                                                          >  2:556222/1‑187 (MQ=255)
                                                                                            aaCCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAaga                                                                                                                                                                                                                                                      <  4:556221/187‑1 (MQ=255)
                                                                                                  aTGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtggt                                                                                                                                                                                                                             >  2:75052/1‑206 (MQ=255)
                                                                                                  aTGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCggggcaggggca                                                                                                                                                        >  2:480917/1‑275 (MQ=255)
                                                                                                  aTGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGCGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGGAGGGATGTGGCAGCATTATTAGTCGAACTATACCggggcaggggcg                                                                                                                                                        >  4:448686/1‑274 (MQ=255)
                                                                                                      ggCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGc                                                                                                                                                                                                                         <  4:75051/206‑1 (MQ=255)
                                                                                                         ttCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAAc                                                                                                                                                 >  4:441850/1‑275 (MQ=255)
                                                                                                          tCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAAcc                                                                                                                                                <  4:278798/275‑1 (MQ=255)
                                                                                                           cGCAATAAAAACGTAGTTGCCTGATGCGCTACGCCTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATACGGCAAAAACAACGAGACAATCAAAATACCGGCGTTATGCCGGTATTTTTTTACGTACGAGGTATCAAATGCGTTTAATTTTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCCGCATTATTCGTCGCATTATACCGGGGCCGGGGAATGGCAACCa                                                                                                                                               >  2:109038/1‑275 (MQ=255)
                                                                                                                        tAGTTGCCTGATGCGCTACGCTTATCAGGCCTCCATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGgt                                                                                                                                  <  2:406037/275‑1 (MQ=255)
                                                                                                                        tAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGgt                                                                                                                                  <  4:314504/275‑1 (MQ=255)
                                                                                                                          gTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGgtgt                                                                                                                                <  2:448687/275‑1 (MQ=255)
                                                                                                                          gTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGCCAGGGGCATGGCAACCAATATTAACCCGgtgt                                                                                                                                >  2:484317/1‑275 (MQ=255)
                                                                                                                           ttGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATAAGTCGCACTATACCGGGGCAGGGGCATGCCAACCACTATTATCCCGgtgtc                                                                                                                               >  4:341779/1‑274 (MQ=255)
                                                                                                                             gCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCa                                                                                                                             <  2:441851/275‑1 (MQ=255)
                                                                                                                              ccTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTGTGAAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGAACTATACCGGGGCAGGGGCATGGCAACCAATATTGTCCCGGTGTGCaa                                                                                                                            <  2:389288/275‑1 (MQ=255)
                                                                                                                                  ttgggttACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTAGCAAGACTTTGTGGGTCGGAGAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATggg                                                                                                                        <  2:239648/269‑1 (MQ=255)
                                                                                                                                          cGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGt                                                                                                                >  4:530150/1‑275 (MQ=255)
                                                                                                                                               aTCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATc                                                                                                                                                                                                                                                >  4:93168/1‑142 (MQ=255)
                                                                                                                                               aTCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATTTTCGTGACTc                                                                                                           >  2:253479/1‑275 (MQ=255)
                                                                                                                                                   ggCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAAt                                                                                                                                                                                                                                            <  2:93169/142‑1 (MQ=255)
                                                                                                                                                   ggCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTa                                                                                                                                                                                                                                                                     >  7:722400/1‑117 (MQ=37)
                                                                                                                                                   ggCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTa                                                                                                                                                                                                                                                                     <  8:722400/117‑1 (MQ=255)
                                                                                                                                                           tCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATg                                                                                                                                                                                                                                                                           >  5:390119/1‑103 (MQ=255)
                                                                                                                                                           tCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATg                                                                                                                                                                                                                                                                           <  6:390119/103‑1 (MQ=255)
                                                                                                                                                              ttCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTAt                                                                                              >  2:55061/1‑273 (MQ=255)
                                                                                                                                                              ttCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGt                                                                                            >  4:307180/1‑275 (MQ=255)
                                                                                                                                                              ttCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATCGTGGTGAGCATCATGGTGACCTTGCTGAg                                                                                                                                                                                                     >  2:562615/1‑170 (MQ=255)
                                                                                                                                                                  tGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATCGTGGTGAGCATCATGGTGACCTTGCTGAGTGGc                                                                                                                                                                                                 <  4:562614/170‑1 (MQ=255)
                                                                                                                                                                   gCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGAAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCa                                                                                          <  1:72933/272‑1 (MQ=255)
                                                                                                                                                                         tATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCtt                                                                                  >  2:421065/1‑274 (MQ=255)
                                                                                                                                                                         tATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCAGGATCCTTg                                                                                 <  2:307181/275‑1 (MQ=255)
                                                                                                                                                                           ttGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACTAGCCAATAAAAAGACCGGCGTTATGCCGGTATTGTTTTACGAAAGGGGTATCAAATGCGTTTAATTGTGGTGTGCATCATGGTGACCTTGCTGAGTGACTGTGTCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGAt                                                                               <  2:356305/275‑1 (MQ=255)
                                                                                                                                                                                    cAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATtc                                                                      >  2:88138/1‑275 (MQ=255)
                                                                                                                                                                                                   cGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGTTCCTTGATCTGCCATTCTCTCTGGTTTTTTAt                                                       >  2:333850/1‑275 (MQ=255)
                                                                                                                                                                                                          ggCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATg                                                                                                                                                                                                                   >  3:145422/1‑112 (MQ=255)
                                                                                                                                                                                                            cGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCg                                                                                                                                                                                                                                         >  6:1024999/1‑88 (MQ=255)
                                                                                                                                                                                                            cGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCg                                                                                                                                                                                                                                         <  5:1024999/88‑1 (MQ=255)
                                                                                                                                                                                                              ttCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGa                                                                                                                                                                                                               <  1:145422/112‑1 (MQ=255)
                                                                                                                                                                                                                cACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATg                                                                                                                          >  1:136222/1‑195 (MQ=255)
                                                                                                                                                                                                                    ccgcATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGAt                                                                                                                      <  3:136222/195‑1 (MQ=255)
                                                                                                                                                                                                                                aaaaCAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATtcatca                          <  2:54291/275‑1 (MQ=255)
                                                                                                                                                                                                                                aaaaCAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATtcatca                          <  2:530151/275‑1 (MQ=255)
                                                                                                                                                                                                                                    caacGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGc                      <  4:88137/275‑1 (MQ=255)
                                                                                                                                                                                                                                          gCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTAt                <  4:55060/275‑1 (MQ=255)
                                                                                                                                                                                                                                                  aaaTACCGGGGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGAtt        <  4:480916/275‑1 (MQ=255)
                                                                                                                                                                                                                                                    aTACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATTAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAg                                                                                                                                                                                                     >  7:741308/1‑84 (MQ=38)
                                                                                                                                                                                                                                                    aTACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAg                                                                                                                                                                                                     <  8:741308/84‑1 (MQ=38)
                                                                                                                                                                                                                                                     tACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGt                                                                                                                                                                                                    <  6:531328/84‑1 (MQ=39)
                                                                                                                                                                                                                                                     tACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCgg                                                                                                                                                                  >  7:51334/1‑118 (MQ=255)
                                                                                                                                                                                                                                                     tACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCgg                                                                                                                                                                  <  8:51334/118‑1 (MQ=255)
                                                                                                                                                                                                                                                     caCCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGt                                                                                                                                                                                                    >  5:531328/2‑84 (MQ=39)
                                                                                                                                                                                                                                                         ggCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGAc                                                                                                                                                                                                              >  7:488535/1‑70 (MQ=255)
                                                                                                                                                                                                                                                         ggCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGAc                                                                                                                                                                                                              <  8:488535/70‑1 (MQ=255)
                                                                                                                                                                                                                                                          gCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCatt                                                                                                                                                                                     >  6:20378/1‑94 (MQ=255)
                                                                                                                                                                                                                                                          gCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCatt                                                                                                                                                                                     <  5:20378/94‑1 (MQ=255)
                                                                                                                                                                                                                                                          gCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATTAACGtt  >  2:357228/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                          
TCTGCATTTCCGTGTCTGGAATTCCCCGGAGAGATATTCTTCATGGGTGAGTTATTACGAAGGGATAAAGCTCAATCCACTGGCATTGCGTAAACCGGATGCGGCTTCGCAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATTAACGTT  >  NC_000913/1113292‑1113814

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: