Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 784,012 T→A 6.0% Q271L (CAG→CTG)  zitB ← zinc efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913784,0120TA 6.0% 534.6 / 17.4 151Q271L (CAG→CTG) zitBzinc efflux system
Reads supporting (aligned to +/- strand):  ref base T (73/68);  new base A (6/3);  total (79/72)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-01

Only 100 of 154 total aligned reads displayed.

ACGTGATTAGCGTCGTTATTTTTGCACTTCAGGGCGTTTACGCCATGTCTCTGGAGTACATCATCCTGACCTTTATTGCGCTCAACGGCAGCCGGATGTGGATCAACAGCGCACGTGAAAGAGGCTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATATCCAGCGATACCGGTGCACCTTCAAGTAATTCATTCACGCTATCTTTCAACAGTCGCCACGCGCTACGCAGAACCAGAAGCGACACCAGTATCGAGAGAATGGGGTCAGCAGGTGTCCAGCC  >  NC_000913/783740‑784283
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                               
aCGTGATTAGCGTCGTTATTTTTGCACTTCAGGGCGTTTACGCCATGTCTCTGGAGTACATCATCCTGACCTTTATTGCGCTCAACGGCAGCCGGATGTGGATCAACAGCGCACGTGAAAGAGGCTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACACCTTCATTCAGGTGGCAGTCTGGCCCAGGACAAGGTTGATATTCCATCTGAATGGTGGCGCGCTCAATCTNNNAGTGATACATCAGGTAGTGTTGGTTCtta                                                                                                                                                                                                                                                                               >  3:142239/1‑273 (MQ=255)
            tcggTAGTTTTGCACTTCAGGGCGTTTACGCCATGTCTCTGGAGTACATCATCCGGACCTTTATTGCGCTCAACGGCAGCCGGATGTGGATCAACAGCGCACGTGAAAGAGGCTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCa                                                                                                                                                                                                                                                                   <  2:114263/275‑1 (MQ=255)
                                                   tGGAGTACATCATCCTGACCTTTATTGCGCTCAACGGCAGCCGGATGTGGATCAACAGCGCACGTGAAAGAGGCTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGGCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCTGGTAGTGTTGGATCTGATCCAACAAGGACTCGTGAGCGTGTGGGGGGATTACCTGCACATGCacggtc                                                                                                                                                                                                                            >  4:437646/1‑270 (MQ=255)
                                                                     ccTTTATTGCGCTCAACGGCAGCCGGATGTGGATCAACAGCGCACGTGAAAGAGGCTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGGGTGGCGGGATCACCTGACCATGCAACGTCATCACCGGCTTCCCCCCt                                                                                                                                                                                                          >  4:179498/1‑275 (MQ=255)
                                                                       tttATTGCGCTCAACGGCAGCCGGATGTGGATCAACAGCGCACGTGAAAGAGGCTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCt                                                                                                                                                                                                          <  2:41577/273‑1 (MQ=255)
                                                                             gggcTCAACGGCAGCGGGATGTTGATAAACCGCGCACGTGAAATAGGCGAAACCCGGCTGTACAATTAAGGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTATGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatcca                                                                                                                                                                                                  <  4:154277/273‑1 (MQ=255)
                                                                             gcgcTCAACGGCAGCCGGATGTGGATCAACAGCGCACGTGAAAGAGGCTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatcca                                                                                                                                                                                                  <  2:476870/275‑1 (MQ=255)
                                                                               gcTCAACGGCAGCCGGATGTGGATCAACAGCGCACGTGAAAGAGGCTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAta                                                                                                                                                                                                <  2:12220/275‑1 (MQ=255)
                                                                                    aCGGCAGCCGGATGTGGATCAACAGCGCACGTGAAAGAGGCTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgt                                                                                                                                                                                           <  2:179500/275‑1 (MQ=255)
                                                                                                      tCAACAGCGCACGTGAAAGAGGCTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCg                                                                                                                                                                          <  2:49966/274‑1 (MQ=255)
                                                                                                        aaCAGCGCACGTGAAAGAGGCTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGcc                                                                                                                                                                                                           >  2:235020/1‑239 (MQ=255)
                                                                                                              gcACGTGAAAGAGGCTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTgg                                                                                                                                                                                                                                                     >  3:103704/1‑191 (MQ=255)
                                                                                                               cACGTGAAAGAGGCTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGTGAGCCCCTGCACATGCAACTTCATCCCCTGCTTCTCACCTACCATCCATCCATGTACATGCTTCACATTGCGGACTCCCgtt                                                                                                                                                                >  4:145420/1‑273 (MQ=255)
                                                                                                                  gTGAAAGAGGCTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCggg                                                                                                                                                                                                                                                 <  1:103704/191‑1 (MQ=255)
                                                                                                                             tCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTg                                                                                                                                                                                                                                                                                >  7:303368/1‑149 (MQ=255)
                                                                                                                                  ggcgcGGTGCCATTAAAGGTGAGGAGGTGAATGGCCATACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGAGATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATacga                                                                                                                                              <  4:239629/273‑1 (MQ=255)
                                                                                                                                   gcgcTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCcaa                                                                                                                                                                                                                                                                          <  8:53643/149‑1 (MQ=255)
                                                                                                                                       tGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAg                                                                                                                                                                                                                                                                      >  5:384019/1‑149 (MQ=255)
                                                                                                                                            cATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatc                                                                                                                                                                                                                                                                 >  5:306651/1‑149 (MQ=255)
                                                                                                                                            cATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatc                                                                                                                                                                                                                                                                 <  8:387968/149‑1 (MQ=255)
                                                                                                                                                   ggTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtgat                                                                                                                                                                                                                                                            >  8:171993/1‑147 (MQ=255)
                                                                                                                                                   ggTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtgat                                                                                                                                                                                                                                                            <  7:171993/147‑1 (MQ=255)
                                                                                                                                                      gatgatGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCgtgt                                                                                                                                                                                                                                                       >  7:448579/1‑149 (MQ=255)
                                                                                                                                                       atgatgTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACAAGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGat                                                                                                                             >  1:41715/1‑270 (MQ=255)
                                                                                                                                                        tgatgTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatc                                                                                                                                                                                                    >  4:338497/1‑198 (MQ=255)
                                                                                                                                                           ttggAATGCCCAGACACCCCTGTNTNCGGGTGGAATTTTGGCCCATGACAAGGTTGATATTCCATCTGAATNGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGNGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACAAGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATATcc                                                                                                                         <  3:41715/267‑1 (MQ=255)
                                                                                                                                                            gtgAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAta                                                                                                                                                                                                <  2:338501/198‑1 (MQ=255)
                                                                                                                                                            gtgAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATATCCAGCGAt                                                                                                                   >  2:79295/1‑275 (MQ=255)
                                                                                                                                                                  ggCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCAcc                                                                                                                                                                                                                                           <  8:31804/149‑1 (MQ=255)
                                                                                                                                                                   gCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCt                                                                                                                                                                                                                                          <  8:448579/149‑1 (MQ=255)
                                                                                                                                                                    ccAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGCTATCCAGCGactgtcgct                                                                                                            >  2:481383/1‑266 (MQ=255)
                                                                                                                                                                    ccAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTCATCCATCAGGTAGTGTTGGGTCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCAc                                                                                                                                                                                                                                            <  7:276320/146‑1 (MQ=255)
                                                                                                                                                                    ccAGAAACTCCTTAATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCAATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATATCCAGCGATACCGGTg                                                                                                            <  4:79294/274‑1 (MQ=255)
                                                                                                                                                                     cAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGc                                                                                                                                                                                                                                        >  8:256021/1‑149 (MQ=255)
                                                                                                                                                                         cacCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACAt                                                                                                                                                                                                                                    <  7:222811/149‑1 (MQ=255)
                                                                                                                                                                          acCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATg                                                                                                                                                                                                                                   >  6:145737/1‑149 (MQ=255)
                                                                                                                                                                            cccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatc                                                                                                                                                                                                                                                                 >  6:414694/1‑117 (MQ=255)
                                                                                                                                                                            cccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatc                                                                                                                                                                                                                                                                 <  5:414694/117‑1 (MQ=255)
                                                                                                                                                                            cccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCa                                                                                                                                                                                                                                 >  5:364331/1‑149 (MQ=255)
                                                                                                                                                                               ttcattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACg                                                                                                                                                                                                                              >  5:333905/1‑149 (MQ=255)
                                                                                                                                                                                 cattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtc                                                                                                                                                                                                                            <  7:302760/149‑1 (MQ=255)
                                                                                                                                                                                 cattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtc                                                                                                                                                                                                                            <  6:333905/149‑1 (MQ=255)
                                                                                                                                                                                    tcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCg                                                                                                                                                                                                                     >  3:126789/1‑153 (MQ=255)
                                                                                                                                                                                      aGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTgg                                                                                                                                                                                                                                                     <  5:102810/119‑1 (MQ=255)
                                                                                                                                                                                      aGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTgg                                                                                                                                                                                                                                                     >  6:102810/1‑119 (MQ=255)
                                                                                                                                                                                      aGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCAc                                                                                                                                                                                                                       >  5:116901/1‑149 (MQ=255)
                                                                                                                                                                                      aGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCAc                                                                                                                                                                                                                       <  8:291449/149‑1 (MQ=255)
                                                                                                                                                                                      aGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCAc                                                                                                                                                                                                                       <  7:274024/149‑1 (MQ=255)
                                                                                                                                                                                        gTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCtt                                                                                                                                                                                                                 <  1:126789/153‑1 (MQ=255)
                                                                                                                                                                                           gCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCt                                                                                                                                                                                                                  >  8:200389/1‑149 (MQ=255)
                                                                                                                                                                                           gCAGTCTGGCCAATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATATCCAGCGATACCGGTGCACCTTCAAGTAATTCATTCACGc                                                                                    <  1:142239/275‑1 (MQ=255)
                                                                                                                                                                                            gatgtgtataaaggACAGGGTTGATATTACATCTGAATGGTGTCTGGCGCAATCTGATAGTGATCCATCAGGTAGTGTTGGAGCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCAACTGCACATGCAACGTCAGCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATATCCAGCGATACCGGTGCACCTTCAAGTAATTCATTCACGc                                                                                    <  4:206126/261‑1 (MQ=255)
                                                                                                                                                                                                 tGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGc                                                                                                                                                                                                            <  8:29807/149‑1 (MQ=255)
                                                                                                                                                                                                  ggCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCcaacaa                                                                                                                                                                                                                                                                       <  6:435745/89‑1 (MQ=255)
                                                                                                                                                                                                    cccATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCt                                                                                                                                                                                                                                          >  5:463367/1‑116 (MQ=255)
                                                                                                                                                                                                    cccATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCt                                                                                                                                                                                                                                          <  6:463367/116‑1 (MQ=255)
                                                                                                                                                                                                     ccATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAc                                                                                                                                                                                                        >  5:68836/1‑149 (MQ=255)
                                                                                                                                                                                                        tGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTa                                                                                                                                                                                                         <  4:521801/145‑1 (MQ=255)
                                                                                                                                                                                                          aCAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatcc                                                                                                                                                                                                   <  6:306651/149‑1 (MQ=255)
                                                                                                                                                                                                               gTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtaca                                                                                                                                                                                              >  7:8632/1‑149 (MQ=255)
                                                                                                                                                                                                               gTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtaca                                                                                                                                                                                              >  6:398986/1‑149 (MQ=255)
                                                                                                                                                                                                               gTTGATATTCCATCTGAATGGTGGCGTACTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtaca                                                                                                                                                                                              >  7:318168/1‑149 (MQ=255)
                                                                                                                                                                                                                    tatTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtac                                                                                                                                                                                         <  6:122219/149‑1 (MQ=255)
                                                                                                                                                                                                                      ttCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtacat                                                                                                                                                                                       >  5:165337/1‑149 (MQ=255)
                                                                                                                                                                                                                          aTCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGCCGCTTCAGTTCAGCGATATCCAGCGATACCGGTGAACCTTCAAGTAATTCATTAACCCTATCTTTCAACAGTCGCCACCCGCTACGCa                                                      >  2:157253/1‑274 (MQ=255)
                                                                                                                                                                                                                              gAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatc                                                                                                                                                                                                    <  2:527939/128‑1 (MQ=255)
                                                                                                                                                                                                                                    tggCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCgg                                                                                                                                                                         <  8:353142/149‑1 (MQ=255)
                                                                                                                                                                                                                                     ggCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGa                                                                                                                                                                        <  7:200389/149‑1 (MQ=255)
                                                                                                                                                                                                                                        gTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACtt                                                                                                                                                                     >  6:386988/1‑149 (MQ=255)
                                                                                                                                                                                                                                           cTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtac                                                                                                                                                                                               >  6:229817/1‑120 (MQ=255)
                                                                                                                                                                                                                                            tCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatcca                                                                                                                                                                                                  >  5:433221/1‑116 (MQ=255)
                                                                                                                                                                                                                                              aaTCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccctcc                                                                                                                                                                                                   <  7:358289/113‑1 (MQ=255)
                                                                                                                                                                                                                                                     tAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCg                                                                                                                                                        <  5:386988/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          aTCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATATCCAGCg                                                                                                                     >  4:29956/1‑179 (MQ=255)
                                                                                                                                                                                                                                                          aTCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATATCCAGCGATACCGGTGCAcc                                                                                                        >  2:484781/1‑192 (MQ=255)
                                                                                                                                                                                                                                                          aTCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATATCCAGCGATACCGGTGCACCTTCAAGTAATTCATTCACGCTATCTTTCAACAGTCGCCACGCGCTACGCAGAACCAGAAGCGACCGCAGTATCGATAGAATg                      >  2:194377/1‑274 (MQ=255)
                                                                                                                                                                                                                                                              aTCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATacg                                                                                                                                               <  7:256021/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              aTCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATATCCAGCGATAc                                                                                                                 <  2:29956/179‑1 (MQ=255)
                                                                                                                                                                                                                                                              aTCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATATCCAGCGATACCGGTGCACCTTCa                                                                                                    <  4:484776/192‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCg                                                                                                                                                        >  8:116299/1‑135 (MQ=255)
                                                                                                                                                                                                                                                                   gTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCg                                                                                                                                                        <  7:116299/135‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttc                                                                                                                                        >  7:344783/3‑149 (MQ=255)
                                                                                                                                                                                                                                                                      ttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTtctc                                                                                                                                                                                                              <  6:164497/77‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTtctc                                                                                                                                                                                                              >  5:164497/2‑78 (MQ=255)
                                                                                                                                                                                                                                                                      ttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGa                                                                                                                                                                        <  7:107698/115‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttca                                                                                                                                       >  5:300966/2‑149 (MQ=255)
                                                                                                                                                                                                                                                                      ttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGAGGCATCACCTGCACATGGAACGTCATCACCGGCTTCTCGCCTACAATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttca                                                                                                                                       >  5:72342/2‑149 (MQ=255)
                                                                                                                                                                                                                                                                       tgtTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATATCCAGCGATACCGGTGc                                                                                                           >  2:200752/1‑176 (MQ=255)
                                                                                                                                                                                                                                                                         ttGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatc                                                                                                                                                                                                    >  3:89622/1‑85 (MQ=255)
                                                                                                                                                                                                                                                                         ttGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCgg                                                                                                                                                                 <  8:101245/120‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ttGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCgg                                                                                                                                                                 <  8:146022/120‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ttGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCgg                                                                                                                                                                 >  7:146022/1‑120 (MQ=255)
                                                                                                                                                                                                                                                                         ttGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCgg                                                                                                                                                       <  7:453116/130‑1 (MQ=255)
                                                                                                                                                                                                                                                                            gATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATATCCAGCGATACCGGTGCACCTTCAAGTAATTCATTCACGCTATCTTTCAACAGt                                                                      >  2:280232/1‑208 (MQ=255)
                                                                                                                                                                                                                                                                             aTCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTtctc                                                                                                                                                                                                              >  6:391149/1‑71 (MQ=255)
                                                                                                                                                                                                                                                                             aTCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAta                                                                                                                                                                                                <  1:89622/85‑1 (MQ=255)
                                                                                                                                                                                                                                                                              tCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATATACAGCGATACCGGTGCACCTGCAAGTAATTCATTCACGCTATCTTTCAACAGTCGCCACGCGCTACGCAGAACCAGAAGCGACCACAGTATGGAGAGAATGGGGTAAGAAGGTGTacagac  >  1:14083/1‑272 (MQ=255)
                                                                                                                                                                                                                                                                              tCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTAAGCg                                                                                                                               >  5:473791/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                               cTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGa                                                                                                                              >  5:391900/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                tGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGat                                                                                                                             >  7:368636/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                               
ACGTGATTAGCGTCGTTATTTTTGCACTTCAGGGCGTTTACGCCATGTCTCTGGAGTACATCATCCTGACCTTTATTGCGCTCAACGGCAGCCGGATGTGGATCAACAGCGCACGTGAAAGAGGCTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATATCCAGCGATACCGGTGCACCTTCAAGTAATTCATTCACGCTATCTTTCAACAGTCGCCACGCGCTACGCAGAACCAGAAGCGACACCAGTATCGAGAGAATGGGGTCAGCAGGTGTCCAGCC  >  NC_000913/783740‑784283

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