Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 328,377 C→G 8.0% L119L (CTG→CTC) ‡ betB ← betaine aldehyde dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913328,3770CG 8.0% 431.2 / 22.8 126L119L (CTG→CTC) ‡betBbetaine aldehyde dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base C (61/54);  new base G (6/4);  total (67/59)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

Only 100 of 130 total aligned reads displayed.

CACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGG  >  NC_000913/328110‑328649
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                
cACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCgg                                                                                                                                                                                                                                                                           <  2:127046/275‑1 (MQ=255)
    accGGCAACCCGGTAAAAACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAGCGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCgta                                                                                                                                                                                                                                                                       <  4:394671/274‑1 (MQ=255)
                         ccGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGCCGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCAc                                                                                                                                                                                                                                                   >  2:153958/1‑274 (MQ=255)
                           gTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTCTACCCAAACGCCCTTTCCCGCACAGGGATCTGGCCGCCTTCCAGCCCC‑GGGATCAGCCCCGCCTCATCCTCCCGCCCCTCCGCACCgg                                                                                                                                                                                                                                                >  4:334385/1‑275 (MQ=255)
                              cGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCg                                                                                                                                                                                                                                                                         <  2:219535/247‑1 (MQ=255)
                               ggacggCCCGCTGCGCTGTAAATTTCCNCCGGCTTTGCGGCGGTAAGCGGGGTNACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCACCGCAGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTANTTACATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTTCAGCACGTCCGCACCGGTa                                                                                                                                                                                                                                              <  3:86340/269‑1 (MQ=255)
                                                        cagcagccTTTAACGCGGGAACGGGGGTACATGCGCTCGGTCTGAAAATCATTGCGTTGCCTGCCGCCAGCGCCGGAGCGGATTTCCACAGGGACACCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GTGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGa                                                                                                                                                                                                                    <  4:153957/267‑1 (MQ=255)
                                                                                   ttgtttggCTGGGTTTGAATTNNNNTTCGNGGGCGGCCGCCAGCGCCGGGGCGTGTGTCCACAGGGCACTCTGGATCGGGTAGTGCCATGCACCAAGCCCGGCCACTACGCCCAGCGGTTCGCGGCGTGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTcc                                                                                                                                                                                        <  3:589/268‑1 (MQ=255)
                                                                                                          aTTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCCCGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCGGGTGGCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTggc                                                                                                                                                                 >  4:200399/1‑273 (MQ=255)
                                                                                                          aTTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGc                                                                                                                                                                  >  2:183231/1‑274 (MQ=255)
                                                                                                              cGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGGCATTGCGTTCa                                                                                                                                                             >  2:461533/1‑275 (MQ=255)
                                                                                                                ttGCCGGCCGCCAGCGCCGGTGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACg                                                                                                                                                           <  2:533684/275‑1 (MQ=255)
                                                                                                                                  ggggCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCgt                                                                                                                                                                                                                                                                        <  7:717104/148‑1 (MQ=255)
                                                                                                                                           ttCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCgg                                                                                                                                                                                                                                                                           >  5:1016240/1‑136 (MQ=255)
                                                                                                                                               acaGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGcc                                                                                                                                                                                                                                                                              >  8:14161/1‑129 (MQ=255)
                                                                                                                                                   ggCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTCCCAGTTTTGCGAGTTCGTAATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAg                                                                                                                         >  2:470220/1‑274 (MQ=255)
                                                                                                                                                   ggCAAACTGGAGCGGGTAGTCCCATGGGCCAACCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAg                                                                                                                         <  4:183230/274‑1 (MQ=255)
                                                                                                                                                             aTCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATc                                                                                                                                                                                                                                                                                  <  7:136541/111‑1 (MQ=255)
                                                                                                                                                             aTCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCgtagta                                                                                                                                                                                                                                                                    <  7:769571/125‑1 (MQ=255)
                                                                                                                                                             aTCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCgtagta                                                                                                                                                                                                                                                                    >  8:769571/1‑125 (MQ=255)
                                                                                                                                                             aTCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTCTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCTTTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATCCGCGAGCGCTCCATGGc                                                                                                              >  1:159392/1‑275 (MQ=255)
                                                                                                                                                                   tAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCGGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGatat                                                                                                                                                                                                                                       >  8:731707/1‑148 (MQ=255)
                                                                                                                                                                      ttCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTcc                                                                                                                                                                                                                                                       >  2:348311/1‑129 (MQ=255)
                                                                                                                                                                         cATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGCTGATCCCgg                                                                                                                                                                                                                                                                           >  6:393730/1‑106 (MQ=255)
                                                                                                                                                                         cATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGCTGATCCCgg                                                                                                                                                                                                                                                                           >  8:695611/1‑106 (MQ=255)
                                                                                                                                                                          aTGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCAc                                                                                                                                                                                                                                                   <  4:348310/129‑1 (MQ=255)
                                                                                                                                                                          aTGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGt                                                                                                                                                                                                                                >  8:1004309/1‑148 (MQ=255)
                                                                                                                                                                           tGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCCGAATATCAACGGCCCGACGCACAAGACGCGAGCGCTCAATGGCGGTCATCGACTcc                                                                                                 >  2:189437/1‑274 (MQ=255)
                                                                                                                                                                           tGCGACAATCCCTGCCCCGACGCCCCGCGGGTGGCGGCGGGGTTACCCCAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGGAGTACTCCAGCACGTCCGGACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTGTCCGGTGTCGAGGGTTTCCAGTGTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGcc                                                                                                 <  4:288745/274‑1 (MQ=255)
                                                                                                                                                                            gcgcCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGccc                                                                                                >  2:323361/1‑274 (MQ=255)
                                                                                                                                                                                  aTCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGttt                                                                                                                                                                                                                        <  7:1004309/148‑1 (MQ=255)
                                                                                                                                                                                       tGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGtt                                                                                                                                                               >  4:95972/1‑200 (MQ=255)
                                                                                                                                                                                           ccaccGCCCGGCGGGTCGTGGCGGGTTAACACAAAGGCCTTTTCACGCAACGGGATCTGGCTGCCTTCCCGCGCC‑GGGATCAGCCCGGCGTAGGACTCCAGCACGTCCGCACCGGTAACGATATCTACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGccc                                                                                   <  4:239951/268‑1 (MQ=255)
                                                                                                                                                                                           aCTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGGTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACg                                                                                                                                                           <  2:95973/200‑1 (MQ=255)
                                                                                                                                                                                            cTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGc                                                                                                                                                                                                              <  7:731707/148‑1 (MQ=255)
                                                                                                                                                                                               cGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTc                                                                                                                                                                                                                       >  6:105023/1‑136 (MQ=255)
                                                                                                                                                                                               cGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTc                                                                                                                                                                                                                       <  5:105023/136‑1 (MQ=255)
                                                                                                                                                                                                    aGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTaa                                                                                                                                                                                                                                             <  6:762797/109‑1 (MQ=255)
                                                                                                                                                                                                     gCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTGGTACTCCAGCACGTCCGCACCGGTAACGata                                                                                                                                                                                                                                        <  8:188907/113‑1 (MQ=255)
                                                                                                                                                                                                        gTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACgg                                                                                                                                                                                                                                 >  8:20657/1‑117 (MQ=255)
                                                                                                                                                                                                        gTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACgg                                                                                                                                                                                                                                 <  7:20657/117‑1 (MQ=255)
                                                                                                                                                                                                          tCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGatat                                                                                                                                                                                                                                       <  8:151739/109‑1 (MQ=255)
                                                                                                                                                                                                          tCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGatat                                                                                                                                                                                                                                       >  7:151739/1‑109 (MQ=255)
                                                                                                                                                                                                             cggcggGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGat                                                                                                                                                                                                                                         >  8:321041/1‑104 (MQ=255)
                                                                                                                                                                                                              gacgggtaaaacaaaaGGCCGTTTCACGCACCGGGACCTGGCTGCCTTCCAGCGCC‑TGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTACGGATATCGACGGTTGAGTTTTCCGAATATGCTTTTCCGGGGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGgcg                                                             <  4:239440/262‑1 (MQ=255)
                                                                                                                                                                                                                  ggTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCCCCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAACGCTCCATGGCGGGACACGCCGCCACGATTTTTTGCCCCTCTTGGGGCCGTTTCCCgggccga                                                          >  4:457557/1‑269 (MQ=255)
                                                                                                                                                                                                                     ataCACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGat                                                                                                                                                                                                                                         <  8:1075338/96‑1 (MQ=255)
                                                                                                                                                                                                                     ataCACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGt                                                                                                                                                                                     >  7:901711/1‑148 (MQ=255)
                                                                                                                                                                                                                      taCACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGCATATCAACGGCCCGACGCAGAATACGCACGCGCTCCATGGCGGTCATCAACGCCCAACTTTTTTGCCCATTCTGGGCGATTTTCACGGCGCGATCGa                                                      >  2:237228/1‑274 (MQ=255)
                                                                                                                                                                                                                       acacAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTAGTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAg                                                                                                                                                                                               <  8:229712/136‑1 (MQ=255)
                                                                                                                                                                                                                       acacAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTAGTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAg                                                                                                                                                                                               >  7:229712/1‑136 (MQ=255)
                                                                                                                                                                                                                        cacaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGa                                                                                                                                                                                                                                    >  8:166144/1‑98 (MQ=255)
                                                                                                                                                                                                                        cacaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGa                                                                                                                                                                                                                                    <  7:166144/98‑1 (MQ=255)
                                                                                                                                                                                                                         acaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTg                                                                                                                                                                                 <  8:901711/148‑1 (MQ=255)
                                                                                                                                                                                                                          caAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACt                                                                                                                                                                                                                                                                  >  6:576271/1‑66 (MQ=255)
                                                                                                                                                                                                                          caAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACt                                                                                                                                                                                                                                                                  <  5:576271/66‑1 (MQ=255)
                                                                                                                                                                                                                          caAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAg                                                                                                                                                                                                                                                              >  5:545243/1‑70 (MQ=255)
                                                                                                                                                                                                                          caAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAg                                                                                                                                                                                                                                                              <  6:545243/70‑1 (MQ=255)
                                                                                                                                                                                                                          caAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACg                                                                                                                                                                                                                                  <  7:174319/98‑1 (MQ=255)
                                                                                                                                                                                                                          caAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAgg                                                                                                                                                                                              >  6:554398/1‑134 (MQ=255)
                                                                                                                                                                                                                          caAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTc                                                                                                                                                                                         <  6:1095342/139‑1 (MQ=255)
                                                                                                                                                                                                                            aaGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGCTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTAATCGCCGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCTGGATCGCAATCCtg                                               >  1:110755/1‑274 (MQ=255)
                                                                                                                                                                                                                             aGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCACCGGCCCGACGCAGAATACGCGAACGCTCCATGGCGGTCATCGCCGCCCAGATTTTTTGCCCCTGCTGGGCCCTTTTCACGGCCCTACCCACATc                                                  >  4:155786/1‑272 (MQ=255)
                                                                                                                                                                                                                             aGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGATGGCTTCATGGCGGTGATCGACGCCCAGGTTTTTTGCCCCCGATGGTCGCTTGTCACGTCGTGATGGGCATCCTc                                               >  2:89078/1‑274 (MQ=255)
                                                                                                                                                                                                                             aGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCt                                                <  4:509344/273‑1 (MQ=255)
                                                                                                                                                                                                                             aGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCGCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGACCGCCATcccc                                               >  4:465194/1‑272 (MQ=255)
                                                                                                                                                                                                                                 cGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAatat                                                                                                                                                                                                                >  6:704659/1‑109 (MQ=255)
                                                                                                                                                                                                                                 cGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAatat                                                                                                                                                                                                                <  5:704659/109‑1 (MQ=255)
                                                                                                                                                                                                                                                   cTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTg                                                                                                                                                                                                                              <  5:91387/77‑1 (MQ=255)
                                                                                                                                                                                                                                                   cTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTg                                                                                                                                                                                                                              >  6:91387/1‑77 (MQ=255)
                                                                                                                                                                                                                                                     ggCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTNACGCAGAATATCAAc                                                                                                                                               >  3:26986/1‑154 (MQ=255)
                                                                                                                                                                                                                                                       cTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCtt                                                                                                                                                                                                            >  6:604845/1‑91 (MQ=255)
                                                                                                                                                                                                                                                       cTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCtt                                                                                                                                                                                                            <  7:618219/91‑1 (MQ=255)
                                                                                                                                                                                                                                                       cTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCtt                                                                                                                                                                                                            >  8:618219/1‑91 (MQ=255)
                                                                                                                                                                                                                                                       cTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCtt                                                                                                                                                                                                            <  5:604845/91‑1 (MQ=255)
                                                                                                                                                                                                                                                         gCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGcc                                                                                                                                           <  1:26986/154‑1 (MQ=255)
                                                                                                                                                                                                                                                            ttCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTGTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACg                <  2:201762/274‑1 (MQ=255)
                                                                                                                                                                                                                                                             tCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGt                                                                                                                                                                                                                          <  7:281882/71‑1 (MQ=255)
                                                                                                                                                                                                                                                             tCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTcc                                                                                                                                                                                                                      >  7:671129/1‑75 (MQ=255)
                                                                                                                                                                                                                                                             tCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGt                                                                                                                                                                                            >  7:124038/1‑101 (MQ=255)
                                                                                                                                                                                                                                                             tCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGt                                                                                                                                                                                            <  8:124038/101‑1 (MQ=255)
                                                                                                                                                                                                                                                             tCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACgg                                                                                                                                             >  8:830062/1‑148 (MQ=255)
                                                                                                                                                                                                                                                              ccAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACg                                                                                                                                              >  7:616173/1‑146 (MQ=255)
                                                                                                                                                                                                                                                              ccAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACg                                                                                                                                              <  8:616173/146‑1 (MQ=255)
                                                                                                                                                                                                                                                              ccAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGc                                                                                                                                            <  7:830062/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                 acgcC‑GGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAgg                                                                                                                                                                                              <  7:399940/94‑1 (MQ=255)
                                                                                                                                                                                                                                                                 acgcC‑GGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAgg                                                                                                                                                                                              >  8:399940/2‑95 (MQ=255)
                                                                                                                                                                                                                                                                 acgcC‑GGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTg                                                                                                                                                                                 <  8:876817/107‑1 (MQ=255)
                                                                                                                                                                                                                                                                 acgcC‑GGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTg                                                                                                                                                                                 >  7:876817/2‑108 (MQ=255)
                                                                                                                                                                                                                                                                 acgcC‑GGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCAt                                                                                                                                                                     >  5:162621/2‑120 (MQ=255)
                                                                                                                                                                                                                                                                 acgcC‑GGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCAt                                                                                                                                                                     <  6:162621/119‑1 (MQ=255)
                                                                                                                                                                                                                                                                  cgcC‑GGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGttt                                                                                                                                                                                          >  7:116179/1‑98 (MQ=255)
                                                                                                                                                                                                                                                                    cc‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAat                                                                                                                                                                                                                  >  6:79495/1‑72 (MQ=255)
                                                                                                                                                                                                                                                                    cc‑GGCATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAat                                                                                                                                                                                                                  <  5:79495/72‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gggATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGa                                                                                                                                                                                                >  5:794005/1‑88 (MQ=255)
                                                                                                                                                                                                                                                                         gATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCgg                                                                                                                                                                                                      >  5:293060/1‑80 (MQ=255)
                                                                                                                                                                                                                                                                         gATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCAt                                                                                                                                                                     <  5:602856/113‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aTCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTCTCCGAATATGCTTTTCCGGTGTCCAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACTCAGAATCCGCGAGCGCTCCATGGCGGTCATCGCCGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGGTCCACCACCTCCGCCCCGTCCGCCGGCCCGGGCCCCAACACTTTCCCCGTGGcccg  >  4:522444/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                            cAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTg                                                                                                                                                                   >  7:1139329/1‑112 (MQ=255)
                                                                                                                                                                                                                                                                            cAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTg                                                                                                                                                                   <  8:1139329/112‑1 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                
CACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCC‑GGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGG  >  NC_000913/328110‑328649

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