Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,334,705 C→T 5.7% intergenic (+24/‑204) ispB → / → sfsB octaprenyl diphosphate synthase/malPQ operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,334,7050CT 5.7% 634.2 / 16.4 193intergenic (+24/‑204)ispB/sfsBoctaprenyl diphosphate synthase/malPQ operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (82/98);  new base T (5/6);  total (88/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 202 total aligned reads displayed.

TGCGATGCATCATGGCACACCAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGG  >  NC_000913/3334435‑3334969
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                        
ggacggggaccaaggCACCCCAGAACAGGCCCATCTGTTCCGTACCGCCATCGAACAGGGTAACTGTCGCCATCTTCCGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcg                                                                                                                                                                                                                                                                       <  4:342541/261‑1 (MQ=255)
     aataacgaacacgacCGGAAGAGACACAGAGGGTCCGTAGCGCAATCGGACAGGGTAACGGTCGCCATCGTCGGGANCCGGTGCGGGAAGCAATGAACGCTTGTGGTTCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGGG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGa                                                                                                                                                                                                                                                                  <  3:163307/261‑1 (MQ=255)
                    cAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACATGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGACTCGCGCACATCGATGTTCAACGCGATCGTTAATGCCCTCCCCTCACCCCTCGCCGCCCGCGGGATGGTTCCAGTAAGCTc                                                                                                                                                                                                                                                  >  4:384669/1‑275 (MQ=255)
                          gggCCCAGATGACACGTACGGCCATCGAACAGGGTAACGGCCGCCATCTTCTGGAACCGGTGCTGGAAGCAAGGAAGGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGATAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCAGCCCGCGCAGCGCGCGGGATGGTTCCAGGAAGTTCCATaaa                                                                                                                                                                                                                                            <  4:483120/274‑1 (MQ=255)
                          aGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaa                                                                                                                                                                                                                                            <  2:622946/275‑1 (MQ=255)
                            gCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGACGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcgc                                                                                                                                                                                                                                                                      >  1:163308/1‑247 (MQ=255)
                                      tCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCa                                                                                                                                                                                                                                                >  2:483122/1‑259 (MQ=255)
                                       ccggaacGCAATCGAAAAGGGTAACGGGCGCCAACTTATGGAAGCGGTTCTGGAAGCAAGGAACGCTTGGGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAg                                                                                                                                                                                                                               <  4:226500/269‑1 (MQ=255)
                                                          ggTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACCCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGaaa                                                                                                                                                                                                              >  2:167922/1‑273 (MQ=255)
                                                                                  acgggtagggAAGAAATTGACGTGTGTGGGTGTCTTGAATGGACGCGTCTGCGTGGCGAGGAATAAGTAGGCAAAGGCATCGCAGCGTTAGCGGTGACTCCTGG‑CACGCCATGGCGGGGAGCAGTCATCGGACTCGCGCAAATCGCTGGTCAACGCGATAGTTAATCCCCTCCCCTAATCCCGCGGAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTatatat                                                                                                                                                                                     <  4:112053/266‑1 (MQ=255)
                                                                                    ggTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCa                                                                                                                                                                                                                                                >  2:640080/1‑213 (MQ=255)
                                                                                        cTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaa                                                                                                                                                                                                                                            <  4:640078/213‑1 (MQ=255)
                                                                                                    aaCGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGt                                                                                                                                                                                                                                                        <  4:168009/189‑1 (MQ=255)
                                                                                                                   ctTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      >  2:647980/1‑224 (MQ=255)
                                                                                                                       aaTGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  <  4:647978/224‑1 (MQ=255)
                                                                                                                       aaTGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTagag                                                                                                                                                <  2:384671/274‑1 (MQ=255)
                                                                                                                               actCAGCTGGCCAAGGAAAGAGCGGCCAAAGGCCTAGCCGCGTTACAGGTG‑CGCCCGGAAACCCCTGGGCGAGACGCACTCATCGGCCTGGCGCACATGCCGGGTCAACGCGATCGTTAATCCACACACCTCATCACACGCTGCGCGCGGGATGGTTCCAGTAATTTCCATAAAAACAGTTATTCAGCTCTGTTAAACCTGAAAAATAGGGATGTGCATCTTTATAAATTCTGAATATTCACACTCGTTAAAGGAACTTTTTAGAGCAATAGGcc                                                                                                                                       <  4:538669/273‑1 (MQ=255)
                                                                                                                                 tCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCTGCGCggg                                                                                                                                                                                                                                                                   >  8:626353/1‑149 (MQ=255)
                                                                                                                                  cAGCGTGCCGAGGAGGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCTGCGCggg                                                                                                                                                                                                                                                                   <  6:148968/148‑1 (MQ=255)
                                                                                                                                  cAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATa                                                                                                                                                                                                                                              >  4:599981/1‑169 (MQ=255)
                                                                                                                                   aGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCTGCGCGGGAt                                                                                                                                                                                                                                                                 <  7:626353/149‑1 (MQ=255)
                                                                                                                                     cGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCTGCGCGGGATga                                                                                                                                                                                                                                                               <  8:78118/149‑2 (MQ=255)
                                                                                                                                      gTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGt                                                                                                                                                                                                                                                              >  7:87641/1‑149 (MQ=255)
                                                                                                                                           gAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCa                                                                                                                                                                                                                                                          <  6:383705/148‑1 (MQ=255)
                                                                                                                                              gaagaaGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTaa                                                                                                                                                                                                                                                      <  8:170060/149‑1 (MQ=255)
                                                                                                                                               aagaagCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTaa                                                                                                                                                                                                                                                      <  6:91249/148‑1 (MQ=255)
                                                                                                                                                      aGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATa                                                                                                                                                                                                                                              >  7:581278/1‑149 (MQ=255)
                                                                                                                                                       gACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAt                                                                                                                                                                                                     >  2:588857/1‑189 (MQ=255)
                                                                                                                                                          aaaGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaaa                                                                                                                                                                                                                                           <  5:773637/148‑1 (MQ=255)
                                                                                                                                                           aaGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTcc                                                                                                                                                                                                                                                           <  5:429463/131‑1 (MQ=255)
                                                                                                                                                           aaGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaaaa                                                                                                                                                                                                                                          >  5:34227/1‑148 (MQ=255)
                                                                                                                                                            aGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAaca                                                                                                                                                                                                                                        <  7:903165/149‑1 (MQ=255)
                                                                                                                                                                aTCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTa                                                                                                                                                                                                                                    >  8:490233/1‑149 (MQ=255)
                                                                                                                                                                 tCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAt                                                                                                                                                                                                                                   >  8:440617/1‑149 (MQ=255)
                                                                                                                                                                  cGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAtt                                                                                                                                                                                                                                  <  7:568193/149‑1 (MQ=255)
                                                                                                                                                                  cGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCta                                                                                                                                                                                                                          >  1:33377/1‑157 (MQ=255)
                                                                                                                                                                     aGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAg                                                                                                                                                                                                                               <  7:274025/149‑1 (MQ=255)
                                                                                                                                                                      gCGTTACAGGTG‑CTCCCGNACACCCNTTGNCGAGAAGNACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAac                                                                                                                                                                                                                      <  3:33376/157‑1 (MQ=255)
                                                                                                                                                                       cGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGc                                                                                                                                                                                                                              >  6:490932/1‑148 (MQ=255)
                                                                                                                                                                        gTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcgc                                                                                                                                                                                                                                                                      >  2:656931/1‑107 (MQ=255)
                                                                                                                                                                         ttACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGctc                                                                                                                                                                                                                            <  6:64788/148‑1 (MQ=255)
                                                                                                                                                                          tACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAtt                                                                                                                                                                                                                                  <  2:84978/141‑1 (MQ=255)
                                                                                                                                                                            cAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGa                                                                                                                                                                                                                                                                  <  4:656929/107‑1 (MQ=255)
                                                                                                                                                                             aGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATaa                                                                                                                                                                                                                       >  7:195722/1‑149 (MQ=255)
                                                                                                                                                                              ggTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATaa                                                                                                                                                                                                                       <  6:34227/148‑1 (MQ=255)
                                                                                                                                                                                    tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCt                                                                                                                                                                                                                  >  5:557603/1‑148 (MQ=255)
                                                                                                                                                                                    tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTg                                                                                                                                                                                                                 >  8:42784/1‑149 (MQ=255)
                                                                                                                                                                                      ccGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATg                                                                                                                                                                                                                                                                <  6:676817/100‑1 (MQ=255)
                                                                                                                                                                                            acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTc                                                                                                                                                                                                                                                            <  8:582848/98‑1 (MQ=255)
                                                                                                                                                                                            acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                                                                                                                                                         >  7:798121/1‑149 (MQ=255)
                                                                                                                                                                                             ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCTGCGCGGGATg                                                                                                                                                                                                                                                                <  5:690784/93‑1 (MQ=255)
                                                                                                                                                                                             ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                                                                                                                                                         <  5:13961/148‑1 (MQ=255)
                                                                                                                                                                                               ccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      <  8:195722/149‑1 (MQ=255)
                                                                                                                                                                                                cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      <  6:557603/148‑1 (MQ=255)
                                                                                                                                                                                                cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAt                                                                                                                                                                                                     <  8:83828/149‑1 (MQ=255)
                                                                                                                                                                                                cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCt                                                                                                                                                                                             >  1:1501/1‑157 (MQ=255)
                                                                                                                                                                                                  tGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgt                                                                                                                                                                                                   <  7:832679/149‑1 (MQ=255)
                                                                                                                                                                                                    gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcg                                                                                                                                                                                                                                                                       <  2:118825/79‑1 (MQ=255)
                                                                                                                                                                                                    gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      >  4:667252/1‑144 (MQ=255)
                                                                                                                                                                                                        gaAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  <  2:667254/144‑1 (MQ=255)
                                                                                                                                                                                                          aGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCAACCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCtt                                                                                                                                                                                            <  6:651944/148‑1 (MQ=255)
                                                                                                                                                                                                              cTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      >  2:546492/1‑134 (MQ=255)
                                                                                                                                                                                                               tCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  >  6:618439/1‑137 (MQ=255)
                                                                                                                                                                                                               tCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtatata                                                                                                                                                                                      >  7:403580/1‑149 (MQ=255)
                                                                                                                                                                                                                cATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtatata                                                                                                                                                                                      >  5:295683/1‑148 (MQ=255)
                                                                                                                                                                                                                 aTCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                                                                                                                                                    >  8:531183/1‑149 (MQ=255)
                                                                                                                                                                                                                  tCGGCGTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                                                                                                                                                                   <  7:149089/149‑1 (MQ=255)
                                                                                                                                                                                                                  tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                                                                                                                                                    >  6:168541/1‑148 (MQ=255)
                                                                                                                                                                                                                  tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                                                                                                                                                                   <  8:798121/149‑1 (MQ=255)
                                                                                                                                                                                                                    ggCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTg                                                                                                                                                                                 <  8:87641/149‑1 (MQ=255)
                                                                                                                                                                                                                      ccTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGTCTGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTTCATCTTTTNATATTCTGAATATTCACAATCTTAACCGGAACTTTTTAGAGTAATAGGCCATCTGGAGAATAGTGATTCTCGAANGAAGCAGTGTCTTGAATGAAAGCAAATAAANNAAGAAGGAAGGAAAAAGGGtaaaacattt                                                 >  3:88229/1‑268 (MQ=255)
                                                                                                                                                                                                                      ccTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGaa                                                                                                                                                                               <  7:490233/149‑1 (MQ=255)
                                                                                                                                                                                                                       cTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGaa                                                                                                                                                                               <  5:889742/148‑1 (MQ=255)
                                                                                                                                                                                                                           cgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      >  4:743242/1‑121 (MQ=255)
                                                                                                                                                                                                                               caTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcg                                                                                                                                                                                                                                                                       <  4:29526/52‑1 (MQ=255)
                                                                                                                                                                                                                               caTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  <  2:743244/121‑1 (MQ=255)
                                                                                                                                                                                                                                 tCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAACTTCCATAAAAACACTTATTGAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCGCtctc                                                                                                                                                                    <  8:573657/149‑1 (MQ=255)
                                                                                                                                                                                                                                  cGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGAGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCa                                     <  1:88230/275‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGct                                                                                                                                                                                                                             >  4:447259/1‑89 (MQ=255)
                                                                                                                                                                                                                                        tCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCtat                                                                                                                                                                                                                         <  2:447261/89‑1 (MQ=255)
                                                                                                                                                                                                                                        tCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAgg                                                                                                                                         >  4:800805/1‑169 (MQ=255)
                                                                                                                                                                                                                                           aCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGaa                                                                                                                                                          >  7:730291/1‑149 (MQ=255)
                                                                                                                                                                                                                                           aCGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCTGCGCggg                                                                                                                                                                                                                                                                   <  6:284221/44‑1 (MQ=18)
                                                                                                                                                                                                                                            cgcgATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCa                                                                                                                                      >  2:270150/1‑168 (MQ=255)
                                                                                                                                                                                                                                             gcgATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTa                                                                                                                                                                >  5:739974/1‑141 (MQ=255)
                                                                                                                                                                                                                                                aTCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAg                                                                                                                                  <  4:270149/168‑1 (MQ=255)
                                                                                                                                                                                                                                                  cGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTa                                                                                                                                                   <  8:635756/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  cGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTa                                                                                                                                                   <  8:403580/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        tcccctccccTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCaa                                                                                                                                             >  8:279176/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          ccctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAAt                                                                                                                                            >  5:616401/1‑148 (MQ=255)
                                                                                                                                                                                                                                                           cctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATa                                                                                                                                           >  5:668523/1‑148 (MQ=255)
                                                                                                                                                                                                                                                            ctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTcac                                                                                                                                                                        <  5:649137/118‑1 (MQ=255)
                                                                                                                                                                                                                                                            ctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAg                                                                                                                                          >  6:820800/1‑148 (MQ=255)
                                                                                                                                                                                                                                                              cccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGcc                                                                                                                                       <  7:42784/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCaag      >  2:214760/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                   cATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCa                                                                                                                                   >  6:737285/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                    aTCCCGCGCGCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttt                                                                                                                                                    >  5:21542/1‑130 (MQ=255)
                                                                                                                                                                                                                                                                     tCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCAATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAAATTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCCTCCCGCAGACATCATTGCGGGttttggacgtaagg  >  4:582359/1‑264 (MQ=255)
                                                                                                                                                                                                                                                                         gcgcAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGAt                                                                                                                      >  2:787312/1‑155 (MQ=255)
                                                                                                                                                                                                                                                                          cgcAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGtat                                                                                                                            <  5:490932/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                        
TGCGATGCATCATGGCACACCAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGG  >  NC_000913/3334435‑3334969

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑