Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,080,621 T→G 5.9% Q203P (CAA→CCA)  fdhE ← formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,080,6210TG 5.9% 653.2 / 29.3 188Q203P (CAA→CCA) fdhEformate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base T (85/91);  new base G (5/6);  total (90/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-02

Only 100 of 190 total aligned reads displayed.

CCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCGCTGCTGACGGACGAGAAATCAGAGGCAAACAGCGCGCTGGCCATATCTTCCAGCTCCTGAGTCGATGCCTTCTCCAGATTCTCAATCACTGCCAGCGCCGGGCCGCTCATTTCAGGTTTCAGCTCA  >  NC_000913/4080351‑4080893
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                
ccTTCTTGCTCGATTCGTGCGTCGAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTAATCGCAGCTTTCGGCTTTAATCGCGGCCTGGTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTgaa                                                                                                                                                                                                                                                                               <  4:181819/274‑3 (MQ=255)
ccTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCa                                                                                                                                                                                                                                                                               <  2:248626/274‑1 (MQ=255)
ccTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCa                                                                                                                                                                                                                                                                               <  2:273348/274‑1 (MQ=255)
agttcttGCTCCCTTCTTGCGTCCAGTACCAGAGAGTCGAGGTCATCTGCTTCGGCTTCGATTTTTGTATCTTTTTCCTGATAGAGAATTTTAAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCGTTTCGTCATCCAGCTACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCa                                                                                                                                                                                                                                                                               <  1:112166/272‑1 (MQ=255)
 cttcttGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGAAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCAcc                                                                                                                                                                                                                                                                             <  2:174936/275‑1 (MQ=255)
                       cAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTTCCAATTTGCACCATGCTGGACACCGGCATAGGAc                                                                                                                                                                                                                                                       >  2:160106/1‑275 (MQ=255)
                                  agGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCggg                                                                                                                                                                                                                                            <  1:121526/275‑1 (MQ=255)
                                    gCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCCGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCCCACTTCCTTCCATTCAGTTTCCCCCACGTTGCCGTTCAGGTAACGTAGCCCCTGCGTGGTTCCAATTTCCCCCCTTCTGGCCCCCGCCCTAGAAcccccaccccggc                                                                                                                                                                                                                                           >  4:388185/1‑263 (MQ=255)
                                            aTCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAAcc                                                                                                                                                                                                                                                      >  4:397284/1‑255 (MQ=255)
                                               tGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCGGCATTTTACGCGGACTACGTGCCATTCAGTTTCACACAGGTTGCCGGGCAGGGAACGTAGCCCCTGCGTGGTACCAAGTTTCACCATCCTGGCCACCGCACCAGAACCACACACCCCGCAatactaacatg                                                                                                                                                                                                                               >  4:28660/1‑267 (MQ=255)
                                                gCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACaa                                                                                                                                                                                                                                                  <  2:397285/255‑1 (MQ=255)
                                                  ttaggcttttttttttGTACCTTGTTCCNNNTAGGGTATTTTCCGGTAAGTGTCGCAGTCATCGCAGTTTTCCGCTTTAATCTCCGCCTGTTCGTCATCCCGCGACCAGTAATGCCGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGc                                                                                                                                                                                                                             <  3:148468/263‑1 (MQ=255)
                                                  tACGGCTTGTATTTTTGGATCTGTTTCCTGATAGCGAATTTTCAGGTAAGTGTCGCAGTCATCGTCGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGAACAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGc                                                                                                                                                                                                                             <  4:244133/274‑1 (MQ=255)
                                                      gCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGc                                                                                                                                                                                                                         <  4:194222/274‑1 (MQ=255)
                                                           gATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTAtt                                                                                                                                                                                                                   <  2:289937/275‑1 (MQ=255)
                                                                     tCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTgg                                                                                                                                                                                                                                                                      >  2:110612/1‑214 (MQ=255)
                                                                         tttCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGacac                                                                                                                                                                                                                                                                  <  4:110611/214‑1 (MQ=255)
                                                                           tcccggtaggcaaTTTTTAGGTAAGTGTCGCAGGCATCGCCGCTTTCGGCTTTAACCGGGGCCTTTTCGGCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTGGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGcc                                                                                                                                                                                                   <  4:453125/264‑1 (MQ=255)
                                                                                              ggTAAGTGGCGCAGTCATCGCAGCTTTCGGCTTTAAGCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCAt                                                                                                                                                                                 <  2:388186/274‑1 (MQ=255)
                                                                                                    tgtCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCCGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTCACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCCGGGTTCAGATTGGCCATCTGCGcc                                                                                                                                                                          >  1:58737/1‑275 (MQ=255)
                                                                                                                       ttCGGCTTTTATCGCGGTCGGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAg                                                                                                                                                       <  1:45790/275‑1 (MQ=255)
                                                                                                                        tCGGCTTTAATCGCGGCCTGTTCTTCATCCGGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAg                                                                                                                                                       <  4:111928/274‑1 (MQ=255)
                                                                                                                            cTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTg                                                                                                                                                                                                                                                                                 >  8:206294/1‑148 (MQ=255)
                                                                                                                                 aTCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCa                                                                                                                                                                                                                                                                            <  5:1033112/148‑1 (MQ=255)
                                                                                                                                      ggCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGgcg                                                                                                                                        >  2:532205/1‑275 (MQ=255)
                                                                                                                                            ttCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACAcc                                                                                                                                                                                                                                                                 >  7:488410/1‑148 (MQ=255)
                                                                                                                                                    ccAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCa                                                                                                                                                                                                                                                                            >  7:931166/1‑129 (MQ=255)
                                                                                                                                                    ccAGCGACCAGTAATGCAGTTTGCCGCACTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCa                                                                                                                                                                                                                                                                            <  8:931166/129‑1 (MQ=255)
                                                                                                                                                     cAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGaa                                                                                                                                                                                                                                                        <  8:447630/148‑1 (MQ=255)
                                                                                                                                                      aGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAAc                                                                                                                                                                                                                                                       >  7:1083421/1‑148 (MQ=255)
                                                                                                                                                      aGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAAc                                                                                                                                                                                                                                                       >  7:548733/1‑148 (MQ=255)
                                                                                                                                                              gTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCCCACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAAcc                                                                                                                                                                                                                                               <  7:322395/148‑1 (MQ=255)
                                                                                                                                                                   gCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCa                                                                                                                                                                                                                                          >  8:1094839/1‑148 (MQ=255)
                                                                                                                                                                    cAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTgg                                                                                                                                                                                                                                                                      >  8:336617/1‑119 (MQ=255)
                                                                                                                                                                         tGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCGCTGCTGACGGACGAGAAATCAGAGGc                                                                                                      >  2:207495/1‑274 (MQ=255)
                                                                                                                                                                              ctctGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTg                                                                                                                                                                                                                                    <  6:344369/143‑1 (MQ=255)
                                                                                                                                                                              ctctGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTg                                                                                                                                                                                                                                    >  5:344369/1‑143 (MQ=255)
                                                                                                                                                                               tctGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACca                                                                                                                                                                                                                                                     >  6:579617/1‑125 (MQ=255)
                                                                                                                                                                                               gcATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCa                                                                                                                                                                                                                                                                            >  6:863874/1‑86 (MQ=255)
                                                                                                                                                                                                  ttttACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCGCTGCTGACGGACGAGAAATCAGAGGCAAACAGCGCGCTGGCCATATCTTcc                                                                             <  2:167415/274‑1 (MQ=255)
                                                                                                                                                                                                  ttttACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCGCTGCTGACGGACGAGAAATCAGAGGCAAACAGCGCGCTGGCCATATCTTcc                                                                             <  2:281602/274‑1 (MQ=255)
                                                                                                                                                                                                  ttttACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTANTCAGCGCGGGCTTTGCCGGGTATCCGGTTGGCCCTCTGCGCCCAGTGGCGCGACAGTGCACCCCAGAAAAACGGCGCTTTCTCTCTGCGGGCGGAAGAGAAATc                                                                                                            >  3:109600/1‑243 (MQ=255)
                                                                                                                                                                                                   tttACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTg                                                                                                                                                                                                                                    >  5:1125490/1‑122 (MQ=255)
                                                                                                                                                                                                   tttACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTg                                                                                                                                                                                                                                    <  8:555038/122‑1 (MQ=255)
                                                                                                                                                                                                   tttACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTg                                                                                                                                                                                                                                    <  6:1125490/122‑1 (MQ=255)
                                                                                                                                                                                                     tACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCactca                                                                                                                                                                                                                                                                        >  7:628882/1‑80 (MQ=39)
                                                                                                                                                                                                     tACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCactca                                                                                                                                                                                                                                                                        <  8:628882/84‑5 (MQ=39)
                                                                                                                                                                                                      aCGAGCACGAAGTGACATTCAGGTTCACACAGGTTTCAGTGCAGTTAACGTAACCCCTGAGTGGGGCCAATTTGCACCATGCTGGACACCGTCAGAGAATCACAATCCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCGCTGCTGACGGACGAGAAATCAGAGGCAAACa                                                                                                 <  1:109600/250‑1 (MQ=255)
                                                                                                                                                                                                       cgcgCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTg                                                                                                                                                                                                                              >  6:441416/1‑124 (MQ=255)
                                                                                                                                                                                                           cACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCg                                                                                                                                                                                                  <  7:1094839/148‑1 (MQ=255)
                                                                                                                                                                                                                   gCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAg                                                                                                                                                                                          >  6:746794/1‑148 (MQ=255)
                                                                                                                                                                                                                       ttCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgggg                                                                                                                                                                                               <  6:899230/139‑1 (MQ=255)
                                                                                                                                                                                                                        tCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATa                                                                                                                                                                                                                                                           >  7:1113446/1‑78 (MQ=255)
                                                                                                                                                                                                                           gTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACGCCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAg                                                                                                                                                                                          <  8:51027/140‑1 (MQ=255)
                                                                                                                                                                                                                                 cacaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTgcg                                                                                                                                                                            >  7:1073007/1‑148 (MQ=255)
                                                                                                                                                                                                                                 cacaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATAAGATTg                                                                                                                                                                                      <  6:615504/138‑1 (MQ=255)
                                                                                                                                                                                                                                   caGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGa                                                                                                                                                                                              >  8:645512/1‑128 (MQ=255)
                                                                                                                                                                                                                                   caGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGa                                                                                                                                                                                              <  7:645512/128‑1 (MQ=255)
                                                                                                                                                                                                                                   caGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGcc                                                                                                                                                                          >  7:556558/1‑148 (MQ=255)
                                                                                                                                                                                                                                     ggTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCgg                                                                                                                                                                                                                                             <  7:186389/79‑1 (MQ=255)
                                                                                                                                                                                                                                       ttGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCt                                                                                                                                                                                                                                                                        >  7:824452/1‑50 (MQ=255)
                                                                                                                                                                                                                                       ttGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCt                                                                                                                                                                                                                                                                        <  8:824452/50‑1 (MQ=255)
                                                                                                                                                                                                                                       ttGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGt                                                                                                                                                                      >  6:653527/1‑148 (MQ=255)
                                                                                                                                                                                                                                       ttGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGt                                                                                                                                                                      <  8:680249/148‑1 (MQ=255)
                                                                                                                                                                                                                                       ttGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGt                                                                                                                                                                      >  7:680249/1‑148 (MQ=255)
                                                                                                                                                                                                                                        tgcagtgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGt                                                                                                                                                                                                                                <  6:1085701/89‑1 (MQ=255)
                                                                                                                                                                                                                                        tgcagtgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGt                                                                                                                                                                                                                                >  5:1085701/1‑89 (MQ=255)
                                                                                                                                                                                                                                        tgcagtgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTa                                                                                                                                                                     >  8:617476/1‑148 (MQ=255)
                                                                                                                                                                                                                                            gtgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGc                                                                                                                                                                 <  8:1073007/148‑1 (MQ=255)
                                                                                                                                                                                                                                             tgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGt                                                                                                                                                                                                                      <  7:622019/94‑1 (MQ=255)
                                                                                                                                                                                                                                             tgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgggg                                                                                                                                                                                               >  5:988421/1‑117 (MQ=255)
                                                                                                                                                                                                                                               cagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGAc                                                                                                                                                              <  5:653527/148‑1 (MQ=255)
                                                                                                                                                                                                                                                   tAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGa                                                                                                                                                                                                                                   >  7:42978/1‑75 (MQ=255)
                                                                                                                                                                                                                                                   tAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCAGCGCCCAGTAGAGCGAc                                                                                                                                                              >  8:1074852/1‑144 (MQ=255)
                                                                                                                                                                                                                                                     aCGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTg                                                                                                                                                                                                                              <  5:8706/78‑1 (MQ=255)
                                                                                                                                                                                                                                                        tAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGcc                                                                                                                                                     >  7:278392/1‑148 (MQ=255)
                                                                                                                                                                                                                                                         aGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGccc                                                                                                                                                    >  6:450367/1‑148 (MQ=255)
                                                                                                                                                                                                                                                           ccccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAatat                                                                                                                                                                                                                                      >  8:184722/1‑64 (MQ=255)
                                                                                                                                                                                                                                                             ccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTg                                                                                                                                                                                                                                    <  8:565244/64‑1 (MQ=255)
                                                                                                                                                                                                                                                               tGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTgcg                                                                                                                                                                            >  5:58359/1‑118 (MQ=255)
                                                                                                                                                                                                                                                                  atggtgCAAATTGGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTg                                                                                                                                                                              <  5:648085/112‑1 (MQ=255)
                                                                                                                                                                                                                                                                  atggtgCAAATTGGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTg                                                                                                                                                                              >  8:310124/2‑113 (MQ=255)
                                                                                                                                                                                                                                                                   tggtgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTag                                                                                                                                                                    <  7:9790/122‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tggtgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTag                                                                                                                                                                    >  8:9790/1‑122 (MQ=255)
                                                                                                                                                                                                                                                                   tggtgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAg                                                                                                                                                            >  8:1003364/1‑130 (MQ=255)
                                                                                                                                                                                                                                                                   tggtgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAg                                                                                                                                                            <  7:1003364/130‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tggtgCAAATTGGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCg                                                                                                                                                                                                  <  5:544661/92‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ggtgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAg                                                                                                                                                                                                                <  5:769762/77‑1 (MQ=255)
                                                                                                                                                                                                                                                                      tgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATACTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCTCTGCTGACGGACGAGAAATCAGCGGCAAACAGCGCGCTGGCCCACTCTTCCAGCCCCTTAGTCGATCCCTTCTCCACATTCCCAATCACTGCCCGCCCCGGGCCGCACATTTccgcctt         >  4:331385/1‑268 (MQ=255)
                                                                                                                                                                                                                                                                       gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGTCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTg                                                                                                                                                                              <  7:961541/108‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGAt                                                                                                                                                                                        >  6:198506/1‑98 (MQ=255)
                                                                                                                                                                                                                                                                       gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGAt                                                                                                                                                                                        <  5:198506/98‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTaga                                                                                                                                                                   <  5:377903/119‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTaga                                                                                                                                                                   >  6:377903/1‑119 (MQ=255)
                                                                                                                                                                                                                                                                       gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTGGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTg                                                                                                                                                                              >  8:961541/1‑108 (MQ=255)
                                                                                                                                                                                                                                                                        ccAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCGCTGCTGACGGACgag                                                                                                                 >  4:406331/1‑168 (MQ=255)
                                                                                                                                                                                                                                                                        caaaTTTCCACCAGGCGTGAACCCGGCATAAAACCAAAAACCTGGAAATATTGAAGTTGTTCGCCTTATTCAGGGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCGCTGCTGACGGACGAGAAATCAGAGGCAAACAGCGCGCTGGCCATATCTTCCAGCTCCTGAGTCGATGCCTTCTCCAGATTCTCAATCACTGCCAGCGCCGGGCCGCTCATTTCAGGTTtcag      <  4:532204/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aaTTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgggg                                                                                                                                                                                               >  8:803941/1‑88 (MQ=255)
                                                                                                                                                                                                                                                                          aaTTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgggg                                                                                                                                                                                               <  7:803941/88‑1 (MQ=255)
                                                                                                                                                                                                                                                                            tttGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCGCTGCTGACGGACGAGAAAt                                                                                                             <  2:406332/168‑1 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                
CCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCGCTGCTGACGGACGAGAAATCAGAGGCAAACAGCGCGCTGGCCATATCTTCCAGCTCCTGAGTCGATGCCTTCTCCAGATTCTCAATCACTGCCAGCGCCGGGCCGCTCATTTCAGGTTTCAGCTCA  >  NC_000913/4080351‑4080893

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