Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,879,473 A→C 22.3% V38G (GTG→GGG)  yoaG ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,879,4730AC22.3% 161.9 / 13.2 81V38G (GTG→GGG) yoaGuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (48/12);  new base C (2/17);  total (51/30)
Fisher's exact test for biased strand distribution p-value = 7.39e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACAGAATGCCCATTGTTTCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGTATCGTCAGTCATGGCTTCAATGTTGTGCCA  >  NC_000913/1879241‑1879705
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                         
tGACAGAATGCCCATTGTTTCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCTGGGGGGCAAGACTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATAAACGATCGGTGATAATTACCAACCACTAACATCCAGTTCGTTTTCCGGGGCAAAAAAAAGAACGGAATTACACACAAATTTTATCACTTcg                                                                                                                                                                                                                                      >  2:319788/1‑237 (MQ=255)
 gACAGAATGCCCATTGTTTCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAAAACGTACGGAATTCACCAAACCTTTTATCACttaac                                                                                                                                                                                                                                     >  2:1337732/1‑235 (MQ=255)
             aTTGTTTCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATAAACGATCAGTGATATTTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGAATTCACCAGTGCTTTTGTAACTTCACCCGCCATTTCTg                                                                                                                                                                                                                        >  1:403486/1‑238 (MQ=255)
              ttGTTCCCTCCCTTATCTGTTTTGCATTCGGGGGTTTGTTGTTGAAAAAAAGCATAGCCCAACATCGGGAGGGCAAGATTTGTTGCGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCCCTATCAACGATCAGTGATAATTCCCAACCAAAAACATCATGTTTGTTTTCCGTGTAATAAGAACGTACGGTATTCACCAGATCTTTTATCATTTCAGCCGCCACTTCTGGc                                                                                                                                                                                                                      <  1:1100991/239‑1 (MQ=255)
                  ttCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCAcc                                                                                                                                                                                                                   >  1:268609/1‑238 (MQ=255)
                   tCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAAAGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCa                                                                                                                                                                                                                  >  1:472460/1‑238 (MQ=255)
                    ccTCCCTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAAATTTTTTATCACTTCAGCCGCCACTTCTGGCAcaaca                                                                                                                                                                                                                >  2:451580/1‑235 (MQ=255)
                            tttttttttGTGTTAGGGGGTTAGTCGATAAAAAAACGCATACCACAACATCGGGAGGGCAAGTTTTGGGACGAGCATCACGGAGGTTTTTTTTCCGATGGCACAGAATTGCCCCCAACACCAACCAGTAAAAATTCCAAACCACAACCTTCATGTTCTTTTTCCGTGTCATAAGAACGTCCGGTATTCACCAGATTTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCAt                                                                                                                                                                                                        <  2:852696/234‑1 (MQ=255)
                            aacTGTTTTGCATTAGGGGGTTAGTCGCTGAAAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCCCCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCAt                                                                                                                                                                                                        <  1:1339102/237‑1 (MQ=255)
                                gTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCCCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGCCATACGAACGGACGGTATTCACCAACTCTTTTATCCCTTCCGCCGCCAATTCTGGCACCCGCCAAACCACCGGc                                                                                                                                                                                                    >  2:230128/1‑239 (MQ=255)
                                  tttGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGCCATAAGACCGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGCATTCGGCg                                                                                                                                                                                                   >  1:825211/1‑238 (MQ=255)
                                    tGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGGAGAGCATCACGGAGGTTTTTTTTGCGAGGGCGCAGAAATTGAGCCATCAACGATCAGTGATATTTACCAAACACAAAACACTAGTTCGTTTTCCGGGTACTAAGAACGGACGGAATTGCCCAGATCTTTTATATCTTAAGCAGCAACTTCTGGCACCAGCAAAGTcctgggagtac                                                                                                                                                                                                >  2:1278985/1‑229 (MQ=255)
                                    tGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGTTATTCACCAGATCTTTTATCACTTCCGCCGCCACTTCTTGCCCCACCAAATTCATCGGCgcc                                                                                                                                                                                                 >  1:1005252/1‑236 (MQ=255)
                                    ggaaTTAGCGGGATAGTCGCTGAGAAAAAGCATAGCACAACATCGGGTGGGCAAGATTTGTGTCGAGCATCACGGAGGTTTTTTTTTCGATGGCGCAGAAATTGCGCCATCAACGCTCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGGTTCATAAGAACGTACCGCATTCACCAGATCTTTTATCCCATCAACCGCCACTTATGCGAACAGCAAAGACAtacgggcc                                                                                                                                                                                                 >  2:1278959/4‑231 (MQ=255)
                                    ggCATTAGGGGTTTAGTGGCTGATAAAAAGCAAAGCACAACATCGGGAGGGAAAGATTTGTGCCGAGCATCACGGAGTTTTTTTTTGCGATGGCGCAGAAATTGGTCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGtct                                                                                                                                                                                                <  1:1148202/238‑1 (MQ=255)
                                      cATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCCGCCAAGTCATCGGCGTctct                                                                                                                                                                                              >  1:1320173/1‑239 (MQ=255)
                                           gCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGGTAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAAATCTTTTATTACTTCAGCCGCCACTTCTGGCACCAGACAAGCAAACGGCGTCTCTTTTTc                                                                                                                                                                                         >  2:1224193/1‑239 (MQ=255)
                                           gCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTc                                                                                                                                                                                         >  1:1268647/1‑239 (MQ=255)
                                                         gATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGCCCCTGTTTCATAATCGCCAGaa                                                                                                                                                                            >  1:275204/1‑238 (MQ=255)
                                                                    tAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCCCCAGCAAAGTCATCGGCGTCTCTGTTTCATAAACGACAGAAACGCCCTTGcgg                                                                                                                                                                >  2:1264359/1‑237 (MQ=255)
                                                                     aGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAATCCATCGGCGTCTCTGTTtaataaccacccaaaaacccatggcggt                                                                                                                                                               >  2:502478/1‑217 (MQ=255)
                                                                       gacaacaTGGGGGGGGAAAGTTTTGGGAAGAGCACCGGGGAGTTTTTTTTGGCGATGGCGCAAAAATTGGACCATCAACGATAAGGGATAATTCCCAACCACAAACATCATGTTCTTTGTCCGTGTAAAAAGAACTGACGGTATTCAACAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTa                                                                                                                                                             <  2:357642/238‑1 (MQ=255)
                                                                                 ggAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCCCCAGCAAAGTCATCGGCGTCCCTGTTTCATAATCGACAGAAACGCCATTCCTGTTATTGGGGACgg                                                                                                                                                   >  2:58709/1‑239 (MQ=255)
                                                                                 ggAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGAATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGAACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCCGTTCTTGTTTACg                                                                                                                                                    >  1:1103628/1‑238 (MQ=255)
                                                                                  gggggAAAGTTTTGTGACGGGCATCAGGGGGGTTTTTTTTGCGATGGCGCAAAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTCCGGTGTCATAAGACCGTACGGTATTCACCAGATCTTATATCCCTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACgg                                                                                                                                                   <  1:177064/236‑1 (MQ=255)
                                                                                   gggtcAAGATTGGTAGCGAGCATCACGGAGGTTTGTTGTGCGATGGCGCAGAAATTGCCCCAACAAGGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCATTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGt                                                                                                                                                  <  1:1358/234‑1 (MQ=255)
                                                                                   aGGGCAAGTTTTGTGCCGACATACACGGAGGTTTTTTTTGCGATGGCGCAAAAATTGCCCCACAAACGATCAGTGATAATTACCAACCAAAAACATCATGTTCGTTTCCCGTGTCATAAGAACGTACGGTATTCCCAAGATTTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGt                                                                                                                                                  <  1:1264359/238‑1 (MQ=255)
                                                                                    gggCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTggggcaggtac                                                                                                                                                >  2:629070/1‑230 (MQ=255)
                                                                                    gggCAAAATTTGGGACGACCATCCCGGAGGTTTTTTTTGCGATGGCGAAAAAATTCCGCCACCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCCCTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTc                                                                                                                                                 <  1:16592/238‑1 (MQ=255)
                                                                                     ggCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACCGAAACGCCATTGCTGTTATTGGTGACGGTCAc                                                                                                                                               >  1:802846/1‑239 (MQ=255)
                                                                                      gCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACg                                                                                                                                              >  1:188702/1‑239 (MQ=255)
                                                                                              tttgggacagaatggcggggtttttttttGGGAGGGACGAGAATGTGGCCCATCAATGGTAAGGTATCATTAACAACCAAAAACTACGTGGTCGTTTTCCGGGTAAAAAGACCGTGCGGTATTCCCCAGATCTTTTATCCCTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACg                                                                                                                                       <  2:1080588/218‑1 (MQ=255)
                                                                                               gtgACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGtt                                                                                                                                     >  1:781022/1‑239 (MQ=255)
                                                                                                  aCGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCCAAGTCATCGGCGTCTCTGGTTCATAATCGCCAGAAACGCCCTTGCTGTTATTGGTGACGGTACCGGGATAAGTTGCt                                                                                                                                  >  1:1273364/1‑239 (MQ=255)
                                                                                                             gggggtttttttttGCGAGGGCAAAAAAATTGCCCCAACAACGAAAATTAAAAATCCCCACCCAAAACCTCACTGTTTGTTTTCCGTGTCAAAAGACCTCACGTTTTTCCCAAGATCTTTTATCACTGCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGa                                                                                                                        <  2:825211/233‑1 (MQ=255)
                                                                                                                     ttttttGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATa                                                                                                                                                                                      >  1:349021/1‑169 (MQ=255)
                                                                                                                     ttttttGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGGATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCCAAGTCATCGGAGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGGGGAGGCAACGGGATACGTTGCTTGGCCCagggttacttt                                                                                                                 >  2:806708/1‑228 (MQ=255)
                                                                                                                        tttGCGATGGCGACAAAATTGCCCCATCAACGATCATGAATAATTCCAAACCAAAAACCTCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCCCCAGATCTTTTATCCCTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCg                                                                                                             <  1:109373/239‑1 (MQ=255)
                                                                                                                        ttggCGATGGCGCAAAAATTGCGCCATAAACAATCGTTAATACTAACCAACCACAAACATTATGTTCTTTTCCCGTGTAAAAAGAACGTCCGGTATTCACAAGATCTTTTATCCCTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCg                                                                                                             <  1:601235/236‑1 (MQ=255)
                                                                                                                            cGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCAACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACCGAAACGCCATGGCTGTTATTGGTGACGGTCACGGTATAAGTTGCTTTGCCCATTATTAATTTCCCGTTat                                                                                                         >  2:220064/1‑239 (MQ=255)
                                                                                                                                         attgcccccacaaagtacAGGGAAAATTCCCAACCAAAAACATCATGTTCTTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATTTTTTATCCCTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCg                                                                                            <  1:58709/222‑1 (MQ=255)
                                                                                                                                            gcgcCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCgttg                                                                                         >  1:451902/1‑239 (MQ=255)
                                                                                                                                            gcgcCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCgttg                                                                                         >  1:1305524/1‑239 (MQ=255)
                                                                                                                                                 ctcacattttagtgtaaaTTCCCACCCAAAAACCTCTTGTTCGTTTCCCGTGTCATAAGAACGTACGTTATTCACCAGATCTTTTACCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTgcg                                                                                    <  2:231333/223‑1 (MQ=255)
                                                                                                                                                 aTCAATGATCATTAATACTTACCAACCCAAAACTTCTTTTTCTTTTTCCGTCTCATAAGAACGTACGGTATTCACCAGTTCTTTTATCCCTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTgcg                                                                                    <  2:451902/239‑1 (MQ=255)
                                                                                                                                                          cAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTTTTGCGCACCTTCCa                                                                           >  1:842722/1‑239 (MQ=255)
                                                                                                                                                           aGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATTGGCGTCTCTGTTTCATAATCGACAGAAACACCATTGCTGTTATTGGTGCCGGTCACGGTATACGTTGCTTTGCCCCTGATTCATTTCCCGTAATGAAGGACTTTCCGTTTTTGGGCACCTTCCaa                                                                          >  2:687678/1‑238 (MQ=255)
                                                                                                                                                            gtattAATTACCAACCCAAACCTTTTTTTTCTTTTTCCGTGTCATAAAAACTTCCGGTTTTCACCAGATTTTTTACCCCTTCAGCCGCCCCTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATc                                                                         <  2:1268647/235‑1 (MQ=255)
                                                                                                                                                              taTAATTCCCAACCAAAACCTTCATGTTCGTTTCCCGTGTTAAAAAACCGTCCGGTATTCCCCAGTTTTTTTACACCTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAg                                                                       <  1:1337732/238‑1 (MQ=255)
                                                                                                                                                                       caaccaCAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAgg                                                              >  1:780148/1‑239 (MQ=255)
                                                                                                                                                                       caaccaCAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAgg                                                              >  1:362772/1‑239 (MQ=255)
                                                                                                                                                                         accaCAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTTCGCACCTTCCATCAGGACTTAAGGa                                                             >  2:61212/1‑238 (MQ=255)
                                                                                                                                                                           cacaAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGcc                                                          >  2:337773/1‑239 (MQ=255)
                                                                                                                                                                           caaaaaCCTGATGTTCGTTTCCCGTGTCATAAGACCGTCGGGTTTTCCCCAGATCTTTTAACACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGcc                                                          <  1:181023/236‑1 (MQ=255)
                                                                                                                                                                                acaatctggttctgttctggtctcatcaGAACTTACTGTTTTCCCCAGATCTTTTATCTCTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACga                                                      <  1:1278959/217‑1 (MQ=255)
                                                                                                                                                                                 actctcgtttCCTTTTCCGTGTAAAAAGGCCGTACGGTATGCCCCAGTTTTTTTACCCCTTCAGCCGCCACTTAGGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACgaa                                                     <  1:687678/231‑1 (MQ=255)
                                                                                                                                                                                     atGTTCGTTCCCCGTGTAAAAAGAGCGTACGTTTTTCACCAGATTTTTTATCCGTTCGGCCGCCACTTCCGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGt                                                <  1:230128/239‑1 (MQ=255)
                                                                                                                                                                                         ttCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCCGGAGCCACGAAGTAGCCaa                                             >  1:560716/1‑239 (MQ=255)
                                                                                                                                                                                               ttCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGa                                       >  1:896406/1‑239 (MQ=255)
                                                                                                                                                                                                   gtgtCATACGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCCCCTTCCATCCGGACTTCCGGAGCACCGAAACAATCAACGTTGAAATa                                   >  2:120782/1‑239 (MQ=255)
                                                                                                                                                                                                    tttCATAAGAACGTACGGTTTTCCCCATTTTTTTTACCCATTGAGCCGCCACTTCTGGCCCCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATaa                                  <  2:780148/237‑1 (MQ=255)
                                                                                                                                                                                                      ttataagcagattcGGTATCCCCCAGTTTTTTTATACTTCGAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGt                                <  2:268609/226‑1 (MQ=255)
                                                                                                                                                                                                           aGACCGTTCGGTTTTCACAAGTTTTTTTACACCTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGTATCgt                           <  1:1224193/239‑1 (MQ=255)
                                                                                                                                                                                                                 tcgggTATTCCCCAGATTTTTTACCCCTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTCGAAATAAGTATCgtcagtca                     <  1:502478/236‑1 (MQ=255)
                                                                                                                                                                                                                 tACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATg                                           >  1:1329561/1‑217 (MQ=255)
                                                                                                                                                                                                                 tACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGAATCgtcagcaa                     >  1:294487/1‑236 (MQ=255)
                                                                                                                                                                                                                 tACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGACGTCTCTGCTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATGCATTTCCCGTTATGAATGACTTTCCGTTGTTCTGCACCTTCGACACGGCCTTCAGGATGCCCCAAGAATTTCATGTTGAAACaaattgcgctag                        >  2:1191412/1‑226 (MQ=255)
                                                                                                                                                                                                                  aCGGTTTTCACCAGATTTTTTAGCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGTATCGTCAGTCAt                    <  1:451580/239‑1 (MQ=255)
                                                                                                                                                                                                                   cGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGTATCGAGAGTCATg                   >  2:976124/1‑239 (MQ=255)
                                                                                                                                                                                                                        ttCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCCCGAAGAAGTCCATGTTGAAATAAGATTCGTCAGTCATGGCTTc              >  1:559165/1‑239 (MQ=255)
                                                                                                                                                                                                                           aCCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAGCAAGTCAATGTTGACATAAGTATGGTCAGTCATGGCTTAAAt           >  1:708962/1‑239 (MQ=255)
                                                                                                                                                                                                                                aTCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGCCAGAAACGCCATTGCTGTTTTTGGGGACGGTCACGGTATACGTTGCTTTGACCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCAACGCAACAGTAAATGGTGAAATAAGTAACGGCAGTCACAGTTTCAAttttgt      >  2:787876/1‑239 (MQ=255)
                                                                                                                                                                                                                                aTCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGCCAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTTCGCCCCTTCACTCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGTAATAAGTATCGGAAGTCATGGCTTCatattttg      >  2:1126437/1‑232 (MQ=255)
                                                                                                                                                                                                                                aTCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGTACTTCAGGAGCCACGAAGCAGTCAATGTTGCAATAAGTATCGTCAATCCTGGCTTCAAtgttgt      >  2:1037644/1‑239 (MQ=255)
                                                                                                                                                                                                                                aTCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGTATCGTCAGTCATGGCTTCAAtgttgt      >  1:174468/1‑239 (MQ=255)
                                                                                                                                                                                                                                aTCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGTATCGTCAGTCATGGCTTCAAtgttgt      >  2:1296743/1‑239 (MQ=255)
                                                                                                                                                                                                                                 acTTTTATACCTTGAGCCGCCATTGCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCAGAATCGACAGATACGCCTTTGCTTTTATTGGTGACGTTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGCCTTTCCGTTGTTGCGCTCCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGTATCGTCAGTCATGGCTTCAAtgttg       <  1:787876/236‑1 (MQ=255)
                                                                                                                                                                                                                                  tttttATCCCTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCa                                                                                                                                                <  1:1141084/98‑1 (MQ=255)
                                                                                                                                                                                                                                  cTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCa                                                                                                                                                >  2:1141084/1‑99 (MQ=255)
                                                                                                                                                                                                                                  cTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGCAGAAGTCAATGTTGAAATAAGTATCGCCAGTCATGGCTTCAAtg          >  1:881221/1‑233 (MQ=255)
                                                                                                                                                                                                                                    tttATCCCTTCCCCCGCCCCTTCCGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGTATCGTCAGTCATGGCTTCAATGTTGTGCCa  <  2:1305524/239‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                         
TGACAGAATGCCCATTGTTTCCTCACTTATCTGTTTTGCATTAGCGGGTTAGTCGCTGATAAAAAGCATAGCACAACATCGGGAGGGCAAGATTTGTGACGAGCATCACGGAGGTTTTTTTTGCGATGGCGCAGAAATTGCGCCATCAACGATCAGTGATAATTACCAACCACAAACATCATGTTCGTTTTCCGTGTCATAAGAACGTACGGTATTCACCAGATCTTTTATCACTTCAGCCGCCACTTCTGGCACCAGCAAAGTCATCGGCGTCTCTGTTTCATAATCGACAGAAACGCCATTGCTGTTATTGGTGACGGTCACGGTATACGTTGCTTTGCCCATGATTCATTTCCCGTTATGAATGACTTTCCGTTGTTGCGCACCTTCCATCAGGACTTCAGGAGCCACGAAGAAGTCAATGTTGAAATAAGTATCGTCAGTCATGGCTTCAATGTTGTGCCA  >  NC_000913/1879241‑1879705

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