Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,339,833 G→A 5.4% S78L (TCG→TTG)  mlaF ← ABC transporter maintaining OM lipid asymmetry, ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,339,8330GA 5.4% 361.2 / 6.7 111S78L (TCG→TTG) mlaFABC transporter maintaining OM lipid asymmetry, ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (59/45);  new base A (2/4);  total (61/50)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 124 total aligned reads displayed.

CCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCT  >  NC_000913/3339576‑3340107
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                  
ccGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGCAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGa                                                                                                                                                                                                                                                                   >  1:114446/1‑275 (MQ=255)
          cgcAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCCCa                                                                                                                                                                                                                                                          >  2:166780/1‑274 (MQ=255)
                    tgggatggcgcTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGaaaa                                                                                                                                                                                                                                                >  2:300891/8‑274 (MQ=255)
                         cAGCGCTGCCAGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGAtct                                                                                                                                                                                                                                          <  2:460119/275‑1 (MQ=255)
                             gcTGCACGCCGCGCCATCCAACCGGAAAGTTCCAAAGGCATCAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTATGATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCAcc                                                                                                                                                                                                                                      <  4:300887/275‑1 (MQ=255)
                                                            cAGAAGGCATTATTTTAGCCGCACCCAGNGNGCCCACGGCCTCAAGCTTCATCATCCCCGTACTATGCAACAAGGGNGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACNCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGa                                                                                                                                                                                                       <  3:42799/275‑1 (MQ=255)
                                                                               cgacgcACGGAGCCCCACGCACCCTAGCTNNAACATCAACCGACGATGCAACAATGGCGCGNAGAGTTGGGTATGTTAGCGCAGTGGATAGGCGACGTTGTAAAATACGTTCATATCATTGAACAACGACCCGGACNGAAATAACCTGCTCATCCGGTTGCGCACTGTAGACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTAc                                                                                                                                                                                    <  3:60647/270‑1 (MQ=255)
                                                                                 cTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGCTCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTAcc                                                                                                                                                                                   >  2:12481/1‑274 (MQ=255)
                                                                                  tCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATCCCGGCGCGAACGAGACATCCCCGCAATATTCCCACCATCGAACAGGATCCCACCATGAACTGGTGCGCTTTGCCCGCCAATCAAGCGGACTCGCGTCGTTTTACCGa                                                                                                                                                                                 >  4:129949/1‑275 (MQ=255)
                                                                                                         gCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGcc                                                                                                                                                                                                              <  2:346721/223‑1 (MQ=255)
                                                                                                               tcatcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTGATCCGTTTGCGCACTGTATACAGGCGCGaa                                                                                                                                                                                                                                                                                  >  5:309038/1‑149 (MQ=255)
                                                                                                                catcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGCCTGTATACAGTGCGCAAAc                                                                                                                                                                                                                                                                                 <  8:125995/149‑1 (MQ=255)
                                                                                                                 atcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACg                                                                                                                                                                                                                                                                                >  8:212454/1‑149 (MQ=255)
                                                                                                                 atcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACg                                                                                                                                                                                                                                                                                >  8:332774/1‑149 (MQ=255)
                                                                                                                  tcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGCCTGTATACAGTGCGCAAACgg                                                                                                                                                                                                                                                                               <  5:34897/149‑2 (MQ=255)
                                                                                                                    aCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACgaga                                                                                                                                                                                                                                                                             <  6:289526/149‑1 (MQ=255)
                                                                                                                     ccGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCt                                                                                                                                              <  2:129950/275‑1 (MQ=255)
                                                                                                                     ccGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGCCTGTATACAGTGCGCAAACggat                                                                                                                                                                                                                                                                             >  6:72679/1‑145 (MQ=255)
                                                                                                                     ccGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGCCTGTATACAGTGCGCAAACggat                                                                                                                                                                                                                                                                             <  5:72679/148‑4 (MQ=255)
                                                                                                                             tGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGAACGCCGTGACCTTCCCTCGCGGCAc                                                                                                                                      >  2:510805/1‑275 (MQ=255)
                                                                                                                                aacaaTGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAat                                                                                                                                                                                                                                                                 >  8:50186/1‑149 (MQ=255)
                                                                                                                                aacaaTGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAat                                                                                                                                                                                                                                                                 >  5:135469/1‑149 (MQ=255)
                                                                                                                                aacaaTGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAat                                                                                                                                                                                                                                                                 >  6:186641/1‑149 (MQ=255)
                                                                                                                                aacaaTGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCCCgg                                                                                                                                    >  2:446910/1‑274 (MQ=255)
                                                                                                                                      ggCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATtctc                                                                                                                                                                                                                                                           <  8:416028/149‑1 (MQ=255)
                                                                                                                                       gcgcGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCa                                                                                                                                                                                                                                                          <  8:372736/149‑1 (MQ=255)
                                                                                                                                                     gTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATGACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGAt                                                                                                                                                                                                                                            >  8:60195/1‑149 (MQ=255)
                                                                                                                                                         gTTCGCGCAGTGGATAGTCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCNGAATATTNCCANCATCGANAAGGAt                                                                                                                                                                                                                                            >  3:24497/1‑145 (MQ=255)
                                                                                                                                                             gcgcAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCa                                                                                                                                                                                                                                        <  1:24497/145‑1 (MQ=255)
                                                                                                                                                                 aGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATc                                                                                                                                                                                                                                <  7:332774/149‑1 (MQ=255)
                                                                                                                                                                  gTGGATAGGCGACTTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAAGAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGcc                                                                                                  <  2:395219/274‑1 (MQ=255)
                                                                                                                                                                    ggATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTgg                                                                                                                                                                                                                             >  5:348281/1‑149 (MQ=255)
                                                                                                                                                                    ggATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGcc                                                                                                                                                                             >  4:162151/1‑197 (MQ=255)
                                                                                                                                                                      aTAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTg                                                                                                                                                                                                                           <  5:186641/149‑1 (MQ=255)
                                                                                                                                                                        aGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGAt                                                                                                                                                                         <  2:162153/197‑1 (MQ=255)
                                                                                                                                                                          gCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGAt                                                                                                                                                                                                                       >  7:33318/1‑149 (MQ=255)
                                                                                                                                                                          gCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGAt                                                                                                                                                                                                                       >  6:72181/1‑149 (MQ=255)
                                                                                                                                                                           cGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGAtt                                                                                                                                                                                                                      <  7:50186/149‑1 (MQ=255)
                                                                                                                                                                            gACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATtc                                                                                                                                                                                                                                                             <  6:30196/109‑1 (MQ=255)
                                                                                                                                                                               gTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGaa                                                                                                                                                                                                                                                                                  <  1:76239/85‑1 (MQ=255)
                                                                                                                                                                               gTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGAACGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAAGATTAGCCGAGAAGACACGATTGCCCCGCGTaaaccgt                                                                                     >  4:332175/1‑270 (MQ=255)
                                                                                                                                                                               gTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGCCTGTATACAGTGCGCAAACg                                                                                                                                                                                                                                                                                >  8:143844/1‑87 (MQ=37)
                                                                                                                                                                                 tGTCAAATACGTTCATATCAGTGAACAACGCCCAGGACTGAAAGAACATGCGCATCCTTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCt                                                                                                                                              <  2:322815/215‑1 (MQ=255)
                                                                                                                                                                                  gTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGc                                                                                                                                                                                                               <  7:212454/149‑1 (MQ=255)
                                                                                                                                                                                   tCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATCACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCCGCACGGTCAGGGAACTATTATCGACGATGCCGCGATTGCCACGTGTaaaaattaaatc                                                                                >  1:39390/1‑267 (MQ=255)
                                                                                                                                                                                   tCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACga                                                                                                                                                                                                                                                                               <  8:80836/84‑1 (MQ=255)
                                                                                                                                                                                   tCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACga                                                                                                                                                                                                                                                                               >  7:80836/1‑84 (MQ=255)
                                                                                                                                                                                     aaaTACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGt                                                                                                                                                                                                                            <  5:402135/133‑1 (MQ=255)
                                                                                                                                                                                     aaaTACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGt                                                                                                                                                                                                                            >  6:402135/1‑133 (MQ=255)
                                                                                                                                                                                       aTACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATc                                                                                                                                                                                                          >  5:158925/1‑149 (MQ=255)
                                                                                                                                                                                        tACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAg                                                                                                                                                                                                                                               >  7:379309/1‑111 (MQ=255)
                                                                                                                                                                                        tACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAg                                                                                                                                                                                                                                               <  8:379309/111‑1 (MQ=255)
                                                                                                                                                                                             tCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACTCGTAAAACTGACATCGCGCATATCg                                                                      >  2:273296/1‑275 (MQ=255)
                                                                                                                                                                                                tatCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGt                                                                                                                                                                                                 >  5:384742/1‑149 (MQ=255)
                                                                                                                                                                                                  tCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAg                                                                                                                                                                                               >  8:267472/1‑149 (MQ=255)
                                                                                                                                                                                                    cgataacaacGCCCCAGACTGAAATAACATGCTCATCCGTTTGCGCACGGTATACAGGCGCTAACGGGACATCGCCGGAATATTCTCACCATCGAAAAGGTGCTCACCAGGATCTGGTGCGATTGGCCCGCCAAGCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAAt                                                               <  4:273292/271‑1 (MQ=255)
                                                                                                                                                                                                    aGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTTACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcg                                                                                                                                                                                             >  6:429885/1‑149 (MQ=255)
                                                                                                                                                                                                     gTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgt                                                                                                                                                                                            >  8:273549/1‑149 (MQ=255)
                                                                                                                                                                                                      tGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgtc                                                                                                                                                                                           >  7:17042/1‑149 (MQ=255)
                                                                                                                                                                                                       gAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgtcg                                                                                                                                                                                          <  6:232479/149‑1 (MQ=255)
                                                                                                                                                                                                        aacaacGCCCCGGACCGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgtcgt                                                                                                                                                                                         >  8:125793/1‑149 (MQ=255)
                                                                                                                                                                                                           aacGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCa                                                                                                                                                                  >  1:83597/1‑169 (MQ=255)
                                                                                                                                                                                                             cGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCCAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATAACGACTAAATTCGccccc                                                       >  2:9835/1‑271 (MQ=255)
                                                                                                                                                                                                               ccccGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCAcc                                                                                                                                                                                                                                      <  1:14813/97‑1 (MQ=255)
                                                                                                                                                                                                               ccccGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGAt                                                                                                                                                              <  3:83597/169‑1 (MQ=255)
                                                                                                                                                                                                                cccGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGaaaa                                                                                                                                                                                                                                                >  5:380537/1‑86 (MQ=255)
                                                                                                                                                                                                                 ccGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGAt                                                                                                                                                                                <  7:108497/149‑1 (MQ=255)
                                                                                                                                                                                                                  cGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGaaaa                                                                                                                                                                                                                                                >  1:78891/1‑84 (MQ=255)
                                                                                                                                                                                                                   caaCTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGTCCCGCACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGACGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCt                                                <  2:514111/273‑1 (MQ=255)
                                                                                                                                                                                                                       tGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGc                                                                                                                                                                              <  7:229791/145‑1 (MQ=255)
                                                                                                                                                                                                                       tGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGc                                                                                                                                                                              >  8:229791/1‑145 (MQ=255)
                                                                                                                                                                                                                       tGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGc                                                                                                                                                                              <  5:453857/145‑1 (MQ=255)
                                                                                                                                                                                                                       tGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGc                                                                                                                                                                              >  6:453857/1‑145 (MQ=255)
                                                                                                                                                                                                                         aaaTAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATtc                                                                                                                                                                                                                                                             >  5:51822/1‑64 (MQ=255)
                                                                                                                                                                                                                         aaaTAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATtc                                                                                                                                                                                                                                                             <  6:51822/64‑1 (MQ=255)
                                                                                                                                                                                                                         aaaTAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGa                                                                                                                                                                                                   <  5:281077/122‑1 (MQ=255)
                                                                                                                                                                                                                         aaaTAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGa                                                                                                                                                                                                   >  6:281077/1‑122 (MQ=255)
                                                                                                                                                                                                                           aTAACATGCTCCTCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCTGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATAAtt                                        <  1:56522/275‑1 (MQ=255)
                                                                                                                                                                                                                                   cTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCAAGGTCAGGGAAATCTTATCGAAGATGCAGCGATTGCCAAGCGTAAAACTGACATCGCGATTATAGAATCAATTCTCCACACCCTGCCACACAATTCACCATTc                                >  4:274836/1‑275 (MQ=255)
                                                                                                                                                                                                                                      aTCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTg                                                                                                                                                                                                                    <  8:61186/92‑1 (MQ=255)
                                                                                                                                                                                                                                      aTCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTg                                                                                                                                                                                                                    >  7:61186/1‑92 (MQ=255)
                                                                                                                                                                                                                                      aTCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGc                                                                                                                                                           >  5:307493/1‑149 (MQ=255)
                                                                                                                                                                                                                                          gTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTTACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTg                                                                                                                                                       <  5:429885/149‑1 (MQ=255)
                                                                                                                                                                                                                                            ttGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGAt                                                                                                                                                     <  7:84713/149‑1 (MQ=255)
                                                                                                                                                                                                                                                gcACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGGCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTAAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATAATTCACCCTTCGTCTTGCGTTGAt                   <  2:53217/275‑1 (MQ=255)
                                                                                                                                                                                                                                                gcACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTc                                                                                                                                                 >  7:302990/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    tGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCg                                                                                                                                                                                                                <  8:169630/82‑1 (MQ=255)
                                                                                                                                                                                                                                                    tGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCg                                                                                                                                                                                                                >  7:169630/1‑82 (MQ=255)
                                                                                                                                                                                                                                                     gTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTcg                                                                                                                                            >  6:98535/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         aCAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGAt                                                                                                                                                                         <  2:45622/116‑1 (MQ=255)
                                                                                                                                                                                                                                                           aGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCAc                                                                                                                                      <  6:348281/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               cgcgAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGcc                                                                                                                                                                                                              >  3:130860/1‑74 (MQ=255)
                                                                                                                                                                                                                                                               cgcgAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCAc                                                                                                >  1:25377/1‑184 (MQ=255)
                                                                                                                                                                                                                                                               cgcgAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTTCTCCAATAATCACCCCTCGTTTTGCGTTGCTTTTTTAAGccaggc     >  4:189138/1‑270 (MQ=255)
                                                                                                                                                                                                                                                                gcgAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTc                                                                                                                                  >  8:102005/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 cgAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGNGGCACGGTCAGGGAAATCTTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAANATCGCCACAGCATGCCCCATAATTAACCCCNNCTCTTGCGTTGATTTTATAAATATGgctc   >  3:63952/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                 cgAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACt                                                                   >  2:213094/1‑211 (MQ=255)
                                                                                                                                                                                                                                                                 cgAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAg           >  2:74407/1‑267 (MQ=255)
                                                                                                                                                                                                                                                                 cgAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGgcgc   >  2:378497/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                 cgAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGAGCATATCCACTAAATTCGCCACAGACTGCTCCATAATTCACCCTTCGTCTTGCGTTGATTGTCTa             >  2:42424/1‑265 (MQ=255)
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                  
CCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCT  >  NC_000913/3339576‑3340107

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