Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,553,609 T→A L237I (TTA→ATA)  uxuB → D‑mannonate oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,553,6090TA95.8% 216.8 / NA 74L237I (TTA→ATA) uxuBD‑mannonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  major base A (39/29);  minor base C (3/0);  total (43/31)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGAAGTCCGCGATTGGTTACATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGT  >  NC_000913/4553371‑4553853
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
gcCGAAGTCCGCGATTGGTTACATCGTCGAAGCGTTGCGTCTGCGTCGGGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGAGCTGGCTCAGGAGCGTGACCCGCCGCTGGCGGCATGGATTGAAGAACATGTCACCTTCACGTGAACCCTGGTTAGCCGCATAGTTCCGGCGGCATAGCCAGCAACCCCACAGGAAAtgg                                                                                                                                                                                                                                          >  2:264868/1‑249 (MQ=255)
  cGAAGTCCGCGATTGGTTACATCGTCGAAGCGGTGCGTCTGCGTCGTGAAAAAGGGTTGAAATCGTTTACGGGGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGGCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGTCCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCt                                                                                                                                                                                                                                        <  2:301320/251‑1 (MQ=255)
  cGAAGTCCGCGATTGGTTACAGCGTCGAAGCGTTTCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGGCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCt                                                                                                                                                                                                                                        <  1:355070/251‑1 (MQ=255)
          gcgATTGGTTACATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCt                                                                                                                                                                                                                                <  1:9966/251‑1 (MQ=255)
                tgTTACATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGgtg                                                                                                                                                                                                                           <  1:357400/249‑1 (MQ=255)
                ggTTACATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGgtgt                                                                                                                                                                                                                          <  1:44855/251‑1 (MQ=255)
                ggTTACATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGgtg                                                                                                                                                                                                                           <  1:245689/250‑1 (MQ=255)
                    aCATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTAc                                                                                                                                                                                                                      <  1:146104/251‑1 (MQ=255)
                     cATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGCTAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACg                                                                                                                                                                                                                     >  1:498980/1‑251 (MQ=255)
                      aTCGTCGCAGCGTTGCGTCTGCGTCGTGAGAAAGGGTTGCTAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCCTGGTTGACCGCATCGTTCCGGCGGCGCCGCCAGAAACCATACAGGAAAGTGCTGACCAGCTGGGTGtttccgc                                                                                                                                                                                                                    >  2:46141/1‑247 (MQ=255)
                        cgtcgAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGAc                                                                                                                                                                                                                   <  1:339165/250‑1 (MQ=255)
                                  ttGCGTCTGCGGCGTGAAAAAGGGTTGAAAGCGTTTACGGGGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGGGGCGAAGGTCGCGGTACTGGGGCTGGCTAAGGCGCGTGACCCGCAGATGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGGGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGACACCATCCAGGAAATTGCTGACCAGCTGGGTGTTTAAGGCACGTGCGCCAt                                                                                                                                                                                                        >  2:145255/1‑251 (MQ=255)
                                       gctgcgtaGTGAAAAAGGGATGAAAGCGTTTACGGTGATGTCCTGAGATAACGTGCGTGAAAGCGGTCATGTGGCGACGGTAGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCGTGGATTGAAGAAAATGTCACATTCCCGTGCACCATGGTTGACCGCCTCGTTCCGGCGGCGACGCCAGCCACCCTCCCGGCAATTGCTGACCAGCTGGGGGATTCAGACCCGTGCGCCATTGCCt                                                                                                                                                                                                   >  2:71078/2‑251 (MQ=255)
                                         tgcgtcGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGGGTTTACGACCCGTGCGCCATTGCCTGc                                                                                                                                                                                                 >  1:419758/1‑251 (MQ=255)
                                            gtcgtGCAAAAGGGTTGTAAGCGTTAACGGTGCTGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCGGGCTCAGGAGCGTGACACGCAGCTGGCGGCATGGATTGAAGAAAGTGTCACCTTCCCGTGCACCATGGTTGACCGCAGCGTTCCGGCAGCGACGCCAGAAAACCTACAGGAAATAGCCGACCAGCTGGGTGTTTGAGCCCGGTGCGCCACTGCCCGCGa                                                                                                                                                                                               >  2:237792/1‑250 (MQ=255)
                                              cgtGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATAGCTGACCAGCTGGGCGTTTACGACCCGTGCGCCATTGCCTGCGAAcc                                                                                                                                                                                            >  1:43381/1‑251 (MQ=255)
                                                      aGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTc                                                                                                                                                                                    >  1:454041/1‑251 (MQ=255)
                                                       gggTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCAATTGCCTGCGAACCGTTCCGTc                                                                                                                                                                                    >  1:110966/1‑250 (MQ=255)
                                                        ggTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTccg                                                                                                                                                                                  >  2:491574/1‑249 (MQ=255)
                                                           ggAAAGCGTTTACGGTGATTTCCGTCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTgg                                                                                                                                                                               <  2:403142/250‑1 (MQ=255)
                                                             aaaGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTgggg                                                                                                                                                                             >  1:199072/1‑250 (MQ=255)
                                                              aaGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTg                                                                                                                                                                            >  1:269056/1‑251 (MQ=255)
                                                               aGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCCCCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACCGGAAATTGCTGACCAGCTGGGTGTTTACGCCCCGTGCGCCATTGCCTGCGACCCGTTCCGTCAGggggggg                                                                                                                                                                           >  2:474471/1‑248 (MQ=255)
                                                                gCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTg                                                                                                                                                                                                              >  1:375284/1‑215 (MQ=255)
                                                                gCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTg                                                                                                                                                                                                              <  2:375284/215‑1 (MQ=255)
                                                                   tttACGGAGATGTCCTGCGGTAACGTGCGTGAAAACGGTCCTGTGGCGCAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCAGCAGCTGGCGGCATGGATTGAAGAAAATGTCACATTCCCGTGCACCAGGGCTGGCCGAATCGTTCCGGCGGCGCCCCCAGAAACCGTACCGGGAATTGCTGCCCAGCTGGGTGTTTAAGACCCCGGCGCCATTACCGGCCAACCGTTcgtgccatggtggattgc                                                                                                                                                                       >  2:260454/1‑234 (MQ=255)
                                                                     tACGGTGGTGTCCTGCGATAACGTGCGTGAAAACGGTGATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCGGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTGGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGaa                                                                                                                                                                      <  2:1334/250‑1 (MQ=255)
                                                                         gTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGCAGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGCTTGAAGAAAATGTAACCTTACCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGCAACCATACAGGGAATTGCTGCCCAGCTGGGTGTTTAAGACCCGTGCGCCATTGCCTGCGAACCGTTCTGTCAGTTGGAGATTgtagcta                                                                                                                                                                  >  1:303127/1‑244 (MQ=255)
                                                                         gTGAGGTCCGGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATaa                                                                                                                                                                 <  2:84998/251‑1 (MQ=255)
                                                                                ccTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATCGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCCGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGtt                                                                                                                                                          >  1:61503/1‑251 (MQ=255)
                                                                                   gCGATAACGTGCGTGAAAACGGTCATGGGGCGAAGGTCGCGGTACTGTGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCTCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTaa                                                                                                                                                        <  2:454041/250‑1 (MQ=255)
                                                                                   gCGATAACGTGCGGGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGGGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTGCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTa                                                                                                                                                         <  2:425581/249‑1 (MQ=255)
                                                                                    cGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATg                                                                                                                                                      >  1:171844/1‑251 (MQ=255)
                                                                                    cGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATg                                                                                                                                                      <  2:171844/251‑1 (MQ=255)
                                                                                      aTAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGt                                                                                                                                                    >  1:313024/1‑251 (MQ=255)
                                                                                        aaCGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACCGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGa                                                                                                                                                                       >  2:62018/1‑230 (MQ=255)
                                                                                        aaCGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGa                                                                                                                                                                       <  1:62018/230‑1 (MQ=255)
                                                                                        aaCGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCg                                                                                                                                                  <  1:48192/251‑1 (MQ=255)
                                                                                         aCGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGGGATTGAAGATAACTTCGTTAATGGTCGc                                                                                                                                                 >  1:147014/1‑251 (MQ=255)
                                                                                         aCGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTAACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGCAACCATACAGGAAATTGCTGACAAGCTGGGTGTTTACGAACCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGCTAACTTCGGTTATGGTCGc                                                                                                                                                 >  2:241035/1‑251 (MQ=255)
                                                                                           gtgcgtgAAAACGGTCATGTGGCGAAGGTCGCGGTCCGGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGCTTGAAGAAAATGTCACCTTCCCGTGCACCATGGGTGACCGCATCGTTCCGGCGGCGACGCCAGAAACAATACAGGACATTGCTGACCAGCTGGGTGGTTACGCCCCGTGCGCCATTGCCTGCGCACCGTTCCGTCCGTGGGGGGTTGTAGCTAACTTGGTTAATggaggacc                                                                                                                                               >  2:262930/1‑245 (MQ=255)
                                                                                             gcgtgAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCTCCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCgg                                                                                                                                             <  2:355613/251‑1 (MQ=255)
                                                                                               gtgAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGCAAATGTCACCTTCCCGTGCACCATGGTTGCCCGCATCGTTCCGGCGGCGACGCCAGCAACCATACAGGAAATTGCTGGCCCGCTGGGGGTTTACGAACCGTGCGCCAATGCCTGCGAACCGTGCCGCCAGTGGGGGTATGACGATAACTTCGTTAATGGTCGCCCGGAt                                                                                                                                           >  2:132673/1‑251 (MQ=255)
                                                                                               gtgAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGAt                                                                                                                                           >  1:238003/1‑251 (MQ=255)
                                                                                                                    aGGTCGCGGTACTGGGGCTGGCTCAGGTGCGTGACCCGCAGCTGGCGGCAGGGATTGAAGAAAATGTCTCCCTCCCGTGCACCATGGTTGCCCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTGCGTTAATGGTCTCCCGGATTGGGATAAAGTGGGCGCACAg                                                                                                                      <  2:338061/251‑1 (MQ=255)
                                                                                                                                       ggctCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTtgtg                                                                                                    <  1:142290/250‑1 (MQ=255)
                                                                                                                                                gcgTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCCCCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGaaa                                                                                          <  1:274366/251‑1 (MQ=255)
                                                                                                                                                         gCAGCTGGCGGCCTGGATTGCAGAAAATGTCAACTTCCCGTGCCCCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATCCAGGAACATGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCCTTGCCTGCGCACCGGTCCGTCAGGGGGCGATTGCAGATAACTTCGTTAATGGTCGCCCGGTTTGGGATAATGTGGGCGCACAGTTCGTTGCAGACGttataccattagaacagatgatg                                                                                  >  2:309844/1‑229 (MQ=255)
                                                                                                                                                           aGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTg                                                                               >  2:153991/1‑251 (MQ=255)
                                                                                                                                                            gCTGGCGGCCTGGATTGAAGAAAATGTCACCTTCACGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGCCGCACGAACCCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCGGCTACCCTTTCCGTCTGTGGGTGTTTGAAGTTAACTTCTTTAATGGTCGCCGGGTTTGGGTTAACGTGGACGCACAGTTCGTTGAAGACTTTGTCCCGTTCGAAAGGATACAGCTGc                                                                              >  2:53575/1‑251 (MQ=255)
                                                                                                                                                               ggcggcATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCg                                                                             <  2:237850/249‑1 (MQ=255)
                                                                                                                                                                 gggcATGGATTGAAGAAAATGTCCCCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATg                                                                         <  2:54438/250‑1 (MQ=255)
                                                                                                                                                                        gATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCTCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTGCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCCCAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACg                                                                  <  2:211979/251‑1 (MQ=255)
                                                                                                                                                                                  aaTGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGCCCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCActc                                                        >  2:204044/1‑251 (MQ=255)
                                                                                                                                                                                  aaTGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCCGCGCCGCCAGACACCATACAGGGAATTGCTGACCAGCTGGGGGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGAGCCGTTCGCAATGATGAAGCTGCGTATGCTGAACGGCAGCCActc                                                        >  1:418178/1‑251 (MQ=255)
                                                                                                                                                                                     gTCACCTTCCCGAGCACCATGGTTGACCGCATCGTTCCGGCGGCGCCGCCAGCAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCttt                                                     >  1:406876/1‑251 (MQ=255)
                                                                                                                                                                                           cacccggtcACCAAGGTTGACCTCAATGTTCCGGCGGCGACGCCAGAAACCAAACAATAAATTGCTTACCAGCTTGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGTGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACATTTCGTTGCAGACGTTGTGCCGTTCGAAATGATTAATCGGCGTAGGCTGAACGGCAGCTACGCTTGTCTGGc                                               <  2:215075/243‑1 (MQ=255)
                                                                                                                                                                                            tCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGc                                               <  1:109017/250‑1 (MQ=255)
                                                                                                                                                                                                 ggCACCATGGGTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCt                                         >  2:70141/2‑251 (MQ=255)
                                                                                                                                                                                                   cACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCg                                       >  2:300646/1‑251 (MQ=255)
                                                                                                                                                                                                   cACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCAGCGAACCGTTCCGTCAGTGGGTGATTGAAGATCACTTCGTTAATGGTCGCCCGGGTTGGGATAAAGTGGGCCCACAGTTCGTTGCAGACGTTATGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGAAGCCACTCTTTTCTGGCGGACCTCg                                       >  2:303473/1‑251 (MQ=255)
                                                                                                                                                                                                   cACCATGGATGACCGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGAGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCg                                       >  1:268311/1‑251 (MQ=255)
                                                                                                                                                                                                       aTGGTTGACCGCATCGTTCCGGCGTCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGGTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGTCACTCTTTTCTGGCGTACCTCGGTTa                                   <  1:241035/251‑1 (MQ=255)
                                                                                                                                                                                                               cctcATCGGTCCGGCGGCGTCTCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTcggcg                           <  1:196522/248‑1 (MQ=255)
                                                                                                                                                                                                               ccGCATCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAg                                                              >  1:469985/1‑216 (MQ=255)
                                                                                                                                                                                                                   aTCGTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTACGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCt                        >  1:63323/1‑250 (MQ=255)
                                                                                                                                                                                                                      gTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATg                     >  2:457772/1‑250 (MQ=255)
                                                                                                                                                                                                                      gTTCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATg                     >  1:26950/1‑250 (MQ=255)
                                                                                                                                                                                                                        tCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCCGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGCAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGaa                   >  2:403377/1‑250 (MQ=255)
                                                                                                                                                                                                                        tCCGGCGGCGACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGaa                   <  1:505294/250‑1 (MQ=255)
                                                                                                                                                                                                                             cggcgACGCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCAtt             <  1:438210/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gCCAGAAACCATACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGac       >  1:245589/1‑250 (MQ=255)
                                                                                                                                                                                                                                       aGAAACCATACAGGCAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCCGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGCTTCGTATAACGCGGGCGCACAGTTCGTTGCAGACGTTTTGCATTTCGATATGCTGACTCTGCGTATGCTTAACGGCAGCCACTCTTTTCTGGCTTCCCACGATTACCTCGGCGTCTTTGAAGCCATTCCCgcccccc   >  2:312291/1‑245 (MQ=255)
                                                                                                                                                                                                                                        gAAACCATACAGGAAAGTGCTGACCAGCTGGGGGTTTACGACCCGGGCCCCAGTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAAAGGTCGCCCGGATTGGGATAAAGTGGGCGCGCAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTTCCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGt  <  2:199519/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
GCCGAAGTCCGCGATTGGTTACATCGTCGAAGCGTTGCGTCTGCGTCGTGAAAAAGGGTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGT  >  NC_000913/4553371‑4553853

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