Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,080,616 C→A 5.6% G205C (GGC→TGC)  fdhE ← formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,080,6160CA 5.6% 688.4 / 23.8 197G205C (GGC→TGC) fdhEformate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base C (89/97);  new base A (5/6);  total (94/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.95e-02

Only 100 of 199 total aligned reads displayed.

CGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCGCTGCTGACGGACGAGAAATCAGAGGCAAACAGCGCGCTGGCCATATCTTCCAGCTCCTGAGTCGATGCCTTCTCCAGATTCTCAATCACTGCCAGCGCCGGGCCGCTCATTTCAGGTTT  >  NC_000913/4080342‑4080886
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                              
  ggCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGCGTGGTGCCaa                                                                                                                                                                                                                                                                              >  2:79223/1‑275 (MQ=255)
         ccTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCa                                                                                                                                                                                                                                                                        <  2:273348/274‑1 (MQ=255)
         agttcttGCTCCCTTCTTGCGTCCAGTACCAGAGAGTCGAGGTCATCTGCTTCGGCTTCGATTTTTGTATCTTTTTCCTGATAGAGAATTTTAAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCGTTTCGTCATCCAGCTACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCa                                                                                                                                                                                                                                                                        <  1:112166/272‑1 (MQ=255)
                          tGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCa                                                                                                                                                                                                                                                      <  2:230114/275‑1 (MQ=255)
                                         agagGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTATGTGCCCTTCAGTTCCACACAGGTTGCAGTGCAGGTAACGTAACCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCg                                                                                                                                                                                                                                       >  2:166788/1‑275 (MQ=255)
                                           agGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCggg                                                                                                                                                                                                                                     <  1:121526/275‑1 (MQ=255)
                                             gCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCCGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCCCACTTCCTTCCATTCAGTTTCCCCCACGTTGCCGTTCAGGTAACGTAGCCCCTGCGTGGTTCCAATTTCCCCCCTTCTGGCCCCCGCCCTAGAAcccccaccccggc                                                                                                                                                                                                                                    >  4:388185/1‑263 (MQ=255)
                                             gCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGc                                                                                                                                                                                                                                    >  2:453126/1‑274 (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)
                                                                    gATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTAtt                                                                                                                                                                                                            <  2:289937/275‑1 (MQ=255)
                                                                                  tttCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGacac                                                                                                                                                                                                                                                           <  4:110611/214‑1 (MQ=255)
                                                                                              gaATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGa                                                                                                                                                                                  <  2:183523/275‑1 (MQ=255)
                                                                                                       ggTAAGTGGCGCAGTCATCGCAGCTTTCGGCTTTAAGCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCAt                                                                                                                                                                          <  2:388186/274‑1 (MQ=255)
                                                                                                        gTAAGTGTCGCATGCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCt                                                                                                                                                                        <  1:92641/275‑1 (MQ=255)
                                                                                                           aGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCCCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTTCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCCATATTGACGTTGCCCCCCGACTTCAGCGCGGGCTTTGCCGGGGATCAGGTTGGCCATCTgcc                                                                                                                                                                     >  4:402870/1‑274 (MQ=255)
                                                                                                             tgtCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCCGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTCACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCCGGGTTCAGATTGGCCATCTGCGcc                                                                                                                                                                   >  1:58737/1‑275 (MQ=255)
                                                                                                                        tCGCTGCTTTCGGCTTTAACCGCTGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGa                                                                                                                                                        <  4:185389/275‑1 (MQ=255)
                                                                                                                                ttCGGCTTTTATCGCGGTCGGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAg                                                                                                                                                <  1:45790/275‑1 (MQ=255)
                                                                                                                                ttCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCa                                                                                                                                                                                                                                                                               >  6:1119633/1‑148 (MQ=255)
                                                                                                                                 tCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCaa                                                                                                                                                                                                                                                                              <  6:500709/148‑1 (MQ=255)
                                                                                                                                     cTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTg                                                                                                                                                                                                                                                                          <  6:1039389/148‑1 (MQ=255)
                                                                                                                                     cTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTg                                                                                                                                                                                                                                                                          >  8:206294/1‑148 (MQ=255)
                                                                                                                                          ctcccgGCCTTTTCCTCTCCCAGCGTCCATTATGGCAGTTTGCCGCTCTTCTTGCATTTTCGTCACTTTACGCGCCCAACGTGCCATTCAGTTTCACCCAGGTTGCGTTGCAGGTAACGTAGCCCCTGGGTAATGACACATTGCACCATGCTCGACACCAGCATGGAACCACTAACCAGGCAATGTTGGCGTTGCACGCCGTATTCAGCGCGGGCTTTGCCAGGAATCAGATTGGGCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAgtgaag                                                                                                                                      <  4:432042/271‑6 (MQ=255)
                                                                                                                                          aTCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCa                                                                                                                                                                                                                                                                     <  5:1033112/148‑1 (MQ=255)
                                                                                                                                               ggCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCTCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCTGGCTTTGCCGGGGATCAGATTGGCCATTTGCGCCCAGTAGAGCGACAGTGCAGCCTAGATAAACGgc                                                                                                                                  >  2:159756/1‑274 (MQ=255)
                                                                                                                                                     ttCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACAcc                                                                                                                                                                                                                                                          >  7:488410/1‑148 (MQ=255)
                                                                                                                                                          cATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCAt                                                                                                                                                                                                                                                     <  6:398263/148‑1 (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:548733/1‑148 (MQ=255)
                                                                                                                                                               aGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAAc                                                                                                                                                                                                                                                >  7:1083421/1‑148 (MQ=255)
                                                                                                                                                                    ccAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACaa                                                                                                                                                                                                                                           <  8:1083421/148‑1 (MQ=255)
                                                                                                                                                                       gTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCCCACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAAcc                                                                                                                                                                                                                                        <  7:322395/148‑1 (MQ=255)
                                                                                                                                                                             cAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTgg                                                                                                                                                                                                                                                               >  8:336617/1‑119 (MQ=255)
                                                                                                                                                                              aGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAAtt                                                                                                                                                                                                                                                                            >  5:1066934/1‑105 (MQ=255)
                                                                                                                                                                               gTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCACTGGGGCCCAGTAGCGCGGCAGTGCAGCCCAGATAAACGGCGCTTTAACGCTGCTGACGGACAAGAAATAagag                                                                                                 >  4:8792/1‑275 (MQ=255)
                                                                                                                                                                                       ctctGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTg                                                                                                                                                                                                                             >  5:344369/1‑143 (MQ=255)
                                                                                                                                                                                        tctGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACca                                                                                                                                                                                                                                              <  5:579617/125‑1 (MQ=255)
                                                                                                                                                                                        tctGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACca                                                                                                                                                                                                                                              >  6:579617/1‑125 (MQ=255)
                                                                                                                                                                                            cTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCGCTGCTGACGGACGAGAAATCAGAGGCAAACAgcgcg                                                                                     <  4:207494/274‑1 (MQ=255)
                                                                                                                                                                                                        gcATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCa                                                                                                                                                                                                                                                                     >  6:863874/1‑86 (MQ=255)
                                                                                                                                                                                                        gcATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGaa                                                                                                                                                                                                                                                 >  7:791426/1‑106 (MQ=255)
                                                                                                                                                                                                        gcATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGaa                                                                                                                                                                                                                                                 <  8:791426/106‑1 (MQ=255)
                                                                                                                                                                                                          aTTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCaa                                                                                                                                                                                                                                                                              >  7:126264/1‑75 (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)
                                                                                                                                                                                                            tttACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCggg                                                                                                                                                                                                   >  5:1102641/1‑148 (MQ=255)
                                                                                                                                                                                                               aCGAGCACGAAGTGACATTCAGGTTCACACAGGTTTCAGTGCAGTTAACGTAACCCCTGAGTGGGGCCAATTTGCACCATGCTGGACACCGTCAGAGAATCACAATCCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCGCTGCTGACGGACGAGAAATCAGAGGCAAACa                                                                                          <  1:109600/250‑1 (MQ=255)
                                                                                                                                                                                                                cgcgCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTg                                                                                                                                                                                                                       >  6:441416/1‑124 (MQ=255)
                                                                                                                                                                                                                cgcgCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTg                                                                                                                                                                                                                       <  5:441416/124‑1 (MQ=255)
                                                                                                                                                                                                                       tACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTAt                                                                                                                                                                                                             >  7:653404/1‑127 (MQ=255)
                                                                                                                                                                                                                       tACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTAt                                                                                                                                                                                                             <  8:653404/127‑1 (MQ=255)
                                                                                                                                                                                                                            gCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAg                                                                                                                                                                                   >  6:746794/1‑148 (MQ=255)
                                                                                                                                                                                                                               aTTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGcc                                                                                                                                                                            >  4:168770/1‑152 (MQ=255)
                                                                                                                                                                                                                                ttCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgggg                                                                                                                                                                                        <  6:899230/139‑1 (MQ=255)
                                                                                                                                                                                                                                 tCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATa                                                                                                                                                                                                                                                    <  8:1113446/78‑1 (MQ=255)
                                                                                                                                                                                                                                  cAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTc                                                                                                                                                                                                           >  6:250293/1‑118 (MQ=255)
                                                                                                                                                                                                                                   aGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCt                                                                                                                                                                        <  2:168771/152‑1 (MQ=255)
                                                                                                                                                                                                                                    gTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACGCCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAg                                                                                                                                                                                   <  8:51027/140‑1 (MQ=255)
                                                                                                                                                                                                                                          cacaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTg                                                                                                                                                                               >  5:615504/1‑138 (MQ=255)
                                                                                                                                                                                                                                          cacaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTgcg                                                                                                                                                                     >  7:1073007/1‑148 (MQ=255)
                                                                                                                                                                                                                                          cacaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATAAGATTg                                                                                                                                                                               <  6:615504/138‑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)
                                                                                                                                                                                                                                                 tgcagtgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGt                                                                                                                                                                                                                         <  6:1085701/89‑1 (MQ=255)
                                                                                                                                                                                                                                                     gtgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGc                                                                                                                                                          <  8:1073007/148‑1 (MQ=255)
                                                                                                                                                                                                                                                      tgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGt                                                                                                                                                                                                               >  8:622019/1‑94 (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)
                                                                                                                                                                                                                                                            tAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAg                                                                                                                                                                                   >  6:130640/1‑116 (MQ=255)
                                                                                                                                                                                                                                                              aCGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTg                                                                                                                                                                                                                       <  5:8706/78‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGcc                                                                                                                                              >  7:278392/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                    ccccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAatat                                                                                                                                                                                                                               >  8:184722/1‑64 (MQ=255)
                                                                                                                                                                                                                                                                      ccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTg                                                                                                                                                                                                                             <  8:565244/64‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAg                                                                                                                                                                                   <  7:634680/106‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAg                                                                                                                                                                                   >  8:634680/1‑106 (MQ=255)
                                                                                                                                                                                                                                                                           gtggtgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACg                                                                                                                                    <  6:351780/148‑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)
                                                                                                                                                                                                                                                                             ggtgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAg                                                                                                                                                                                                         <  5:769762/77‑1 (MQ=255)
                                                                                                                                                                                                                                                                               tgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATACTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCTCTGCTGACGGACGAGAAATCAGCGGCAAACAGCGCGCTGGCCCACTCTTCCAGCCCCTTAGTCGATCCCTTCTCCACATTCCCAATCACTGCCCGCCCCGGGCCGCACATTTccgcctt  >  4:331385/1‑268 (MQ=255)
                                                                                                                                                                                                                                                                                gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGTCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTg                                                                                                                                                                       <  7:961541/108‑1 (MQ=255)
                                                                                                                                                                                                                                                                                gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGAt                                                                                                                                                                                 <  5:198506/98‑1 (MQ=255)
                                                                                                                                                                                                                                                                                gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTaga                                                                                                                                                            <  5:377903/119‑1 (MQ=255)
                                                                                                                                                                                                                                                                                gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCt                                                                                                                               >  7:1102884/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTGGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTg                                                                                                                                                                       >  8:961541/1‑108 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                              
CGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCGCTGCTGACGGACGAGAAATCAGAGGCAAACAGCGCGCTGGCCATATCTTCCAGCTCCTGAGTCGATGCCTTCTCCAGATTCTCAATCACTGCCAGCGCCGGGCCGCTCATTTCAGGTTT  >  NC_000913/4080342‑4080886

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