Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,133 T→G S146S (TCT→TCG)  creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,637,1330TG98.7% 237.5 / ‑4.3 74S146S (TCT→TCG) creCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (38/35);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

GGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTC  >  NC_000913/4636888‑4637377
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
gggCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGat                                                                                                                                                                                                                                                  <  1:223812/251‑1 (MQ=255)
 ggCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGatg                                                                                                                                                                                                                                                 >  1:529102/1‑251 (MQ=255)
            ccggCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGAGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGTCGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCATGCCGTTCCCGAAAGTTCGGTGATGTATGTTGCCGc                                                                                                                                                                                                                                      >  2:4276/3‑251 (MQ=255)
                    tAATCCGCTACAACATATCCCGTTTCGCGCCAATATCGGTGGTATTAACAAAGTTCGCATTGA‑TATTCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATCACGTCTGCCTAACGTTGCGTGGTCAGTATGTTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGCCCGGTGATGTATGTTGCCGCTCCGAtt                                                                                                                                                                                                                               <  2:97660/250‑1 (MQ=255)
                    tAATCAGCTACAACATCGCCCGTTTCGCGCCACTATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAACGCAGTTGGGCAGGATTATTCGCGCTGGAACGACGTCTGGCTAACGTTGCGTGGTCAGTTTGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCTGAAAGTTCGGTGATGTATGGTGCCGCGCCggtt                                                                                                                                                                                                                               >  2:270227/1‑247 (MQ=255)
                        cAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGa                                                                                                                                                                                                                          <  1:134192/251‑1 (MQ=255)
                              caacaTCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTc                                                                                                                                                                                                                    <  1:304694/251‑1 (MQ=255)
                                 caTCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTTCCGCGCCGATTATGGACGGCTcgcg                                                                                                                                                                                                                 >  2:260222/1‑251 (MQ=255)
                                 caTCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTcgcg                                                                                                                                                                                                                 >  1:508434/1‑251 (MQ=255)
                                      cgcttttCGCGCCATTTGCGGTTGCATTAACAAAGTGCGCAATGAATAT‑CATGTCAATATGACCGATTCGCAGGTTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGGGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTa                                                                                                                                                                                                            <  2:656271/247‑1 (MQ=255)
                                      ccctgttCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAAGAAAGACGTTGGGCAGGATTTTTCGCGCTGGAATGTCGGCTGGTTAACGTTGCGTGGTCCGTATGGTGCGCGCCGCACGTTGGAAAATCCTGCCGATCCCGAAAGTGCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTa                                                                                                                                                                                                            <  1:594315/246‑1 (MQ=255)
                                          tttCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTgg                                                                                                                                                                                                        >  1:335265/1‑251 (MQ=255)
                                               cgcCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGttt                                                                                                                                                                                                   >  1:14871/1‑251 (MQ=255)
                                                  cAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTg                                                                                                                                                                                                 >  1:229176/1‑250 (MQ=255)
                                                      atCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGatg                                                                                                                                                                                                                                                 >  1:52283/1‑198 (MQ=255)
                                                      atCGGTGGCATTAACAAAGTGCGCAAGGAATAT‑CATGGCTATATGACTGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATGATGCGCGCTGGAATGACGTCTGGCTACCGTTGCGTGGTCAGTAGGGTGCGCGAAGCACGGTGCAAAACCCTGCCGATCCCGAAAGTTCGGTGagg                                                                                                                                                                                                                                                 <  2:52283/198‑3 (MQ=255)
                                                                    aaaaGTGCGAAAGGAATGT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCGGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTGGCAAACCCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAAcg                                                                                                                                                                              <  2:590219/250‑1 (MQ=255)
                                                                     aaaGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCACACCGAAcgc                                                                                                                                                                             >  2:121264/1‑251 (MQ=255)
                                                                            gctaTGAATAT‑CATGTCTATATGACCGATTCGGAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTGGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATg                                                                                                                                                                      <  2:255465/248‑1 (MQ=255)
                                                                              aaTGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑c                                                                                                                                                                   >  1:600317/1‑251 (MQ=255)
                                                                              aaTGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑c                                                                                                                                                                   >  1:196753/1‑251 (MQ=255)
                                                                                     aatcctgtctatctgctcgtttcgcAGGGTAAAGTATTGCTCGATTCGGAAAAGAATGCCGTTGGGAAGGTGTTTTCGCGCTGGAATGATGACTGTCTAACGGTGCGTGGTGAGTAGGGTGCTCGCAGCATGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGGAGGTTGCCGCGCCGATTATGGACGGGGCGCGGCTGATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGT‑CTCCGGt                                                                                                                                                             <  2:602867/229‑1 (MQ=255)
                                                                                        cATGTCTATATGACTGATTCGCAGGTTAAAGTTTTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGTCGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCa                                                                                                                                                           <  2:14871/250‑1 (MQ=255)
                                                                                                      cGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTagcgagc                                                                                                                                            >  1:183838/1‑251 (MQ=255)
                                                                                                        attcgcAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGCTGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGc                                                                                                                                                                                                       >  2:82391/4‑191 (MQ=255)
                                                                                                        attcgcAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGc                                                                                                                                                                                                       <  1:82391/188‑1 (MQ=255)
                                                                                                        agtcgcAGGGTAACGTATTGTGCGATGCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCTCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAAGCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGGAGCGAgcgg                                                                                                                                          <  2:126663/248‑1 (MQ=255)
                                                                                                            gcAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAAGCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGaa                                                                                                                                      <  1:620334/251‑1 (MQ=255)
                                                                                                                   aaaGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATg                                                                                                                               <  1:580329/251‑1 (MQ=255)
                                                                                                                    aaGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATgg                                                                                                                              <  1:260222/251‑1 (MQ=255)
                                                                                                                      gTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTATGGCTAACGTTGAGTGGACAGTATGGTGCGAGCAGCACGTAGCAAAATCCTGCCGATCCCGAAAGTGCGGAGATGTATGTTGACGCGCAGATTATGGACGGCTCGCGGCTTATTGGAGTTTTGAGCGTAGGCAAACCGAAAGCGGCG‑TGGACTCCGGGCCTTAAGCGTCGCGAGAGGCGGATTTTATgggg                                                                                                                            >  2:203695/1‑250 (MQ=255)
                                                                                                                       tATTGTTCGATGCGGCAAATAAAACCGTTGGGCAGGATTATTCGCGCTGGACTGACGTCTGGATAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGACAGTTCGGTGATGTATGTTGCCGAGACGATTATGGACGGCCCGCGGCTTAATGGCGTTTTGATCGTAGGCTACACGCACGCGGAGCAGG‑CTACGGTCATGACGCGTAGCGAGCCGCGAATTTTATTGGcc                                                                                                                           >  2:126654/1‑251 (MQ=255)
                                                                                                                        aTTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGGTCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCTTTAAGCGTAGCGAGCGGCGAATTTTAGGGGcc                                                                                                                           >  2:337048/1‑250 (MQ=255)
                                                                                                                          tgttcctttCGGAAAAAAAAGCCGTGGGGCAGGATTATTCGCGCTGGAATAACTTCTGTCTAACGTTGCTTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAg                                                                                                                         <  2:508434/243‑1 (MQ=255)
                                                                                                                              cGATTCGGCAAATAAAGACGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGAAAACTCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCACGGCTTATTGGCGTTTTTAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAgcaccc                                                                                                                    >  2:585194/1‑247 (MQ=255)
                                                                                                                                  tCGGCAAATAAAGCCGTTGGGCAGGATTACTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTGGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCAttt                                                                                                                 <  1:120465/250‑1 (MQ=255)
                                                                                                                                   cGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTttg                                                                                                               <  1:229325/251‑1 (MQ=255)
                                                                                                                                   cGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTttg                                                                                                               >  1:335950/1‑251 (MQ=255)
                                                                                                                                     gCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTttgct                                                                                                             >  1:310793/1‑249 (MQ=255)
                                                                                                                                          tAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTgggg                                                                                                         >  1:673389/1‑250 (MQ=255)
                                                                                                                                          tAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGAATATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGAAGGCCAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTAGTGGCAAGCGCCATTTTGCTggggg                                                                                                        >  2:68293/1‑250 (MQ=255)
                                                                                                                                             ggCCGTTGGGCAGGATTATTCGGGTTGGAATGAGGTCTGGCTAGCGTTGCGTGGTTAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTg                                                                                                     <  2:183773/250‑1 (MQ=255)
                                                                                                                                             aGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGAATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGAGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTg                                                                                                     >  1:246817/1‑251 (MQ=255)
                                                                                                                                               ccGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGc                                                                                                    >  1:116988/1‑250 (MQ=255)
                                                                                                                                                    gggcAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGt                                                                                              >  1:676809/4‑251 (MQ=255)
                                                                                                                                                         ggATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTg                                                                                         >  1:283024/1‑251 (MQ=255)
                                                                                                                                                            ttattCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGGGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTTTGGGCCAGCGCCATTTTGCTGGGGTTTGCACTGGTGATTgtcg                                                                                      >  2:340328/1‑248 (MQ=255)
                                                                                                                                                                tCGCGCTGGAATGACGTCTGGCTCACGTTGCGTGGTCAGTATGGTGCGCGAAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCACGCGGATTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCCTTAAGCGTAGCGAGCGTCGATTTTTATGGGCCAGCGCCCTTTTGCAGGGGATTGACCGGGTGATTGGCGCAgg                                                                                  >  2:631409/1‑251 (MQ=255)
                                                                                                                                                                      tGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGAACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTAGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGt                                                                            >  1:76117/1‑251 (MQ=255)
                                                                                                                                                                       ggAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATgg                                                                             <  1:181025/249‑1 (MQ=255)
                                                                                                                                                                        gaatgactttgtgcataCGTTGCGGGGGCAGTTTGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGGGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGttt                                                                          <  2:489152/235‑1 (MQ=255)
                                                                                                                                                                        gAATGACGTCTGTCTAACGTTGCGTGGGCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGGTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGttt                                                                          <  1:631409/251‑1 (MQ=255)
                                                                                                                                                                         aaTGACGGCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTtg                                                                         <  1:572861/251‑1 (MQ=255)
                                                                                                                                                                            gACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTAATTGGCCCAGGCATGGTTtggtg                                                                      >  1:237566/1‑251 (MQ=255)
                                                                                                                                                                                    gttAACGTGGAGTGGTCAGTATGGTGCGCGAGGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAAcc                                                              <  2:673389/249‑1 (MQ=255)
                                                                                                                                                                                         cgtggcgtgtcaaGTATGGGGCGCGCAGCCCGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGGTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATGTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGctcc                                                         <  2:310793/239‑2 (MQ=255)
                                                                                                                                                                                               gTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGcc                                                   <  1:121264/251‑1 (MQ=255)
                                                                                                                                                                                                tgggcAGTATGGTGCGCGCAGCACGTTGCCAAATCCTGCAGCTCCCGCACGTTCGGGGATGTATGTTGCCGCGCCGATTCTGGACGGCGCGCGGCTTCTTGGCGTTTTGAGCGGCGGCCAACCGAACGCGGCGATGG‑CTCCGGTCCTTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCAATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGGTTGGTGGATCAACCGCGCCATTGccc                                                  >  2:130575/5‑251 (MQ=255)
                                                                                                                                                                                                    cAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGTTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGAt                                                                                                                                                                       <  1:399393/131‑1 (MQ=255)
                                                                                                                                                                                                    cAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCt                                              <  1:278620/251‑1 (MQ=255)
                                                                                                                                                                                                        aTGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACt                                          >  1:151128/1‑251 (MQ=255)
                                                                                                                                                                                                        aTGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACt                                          >  1:226275/1‑251 (MQ=255)
                                                                                                                                                                                                        aTGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACt                                          >  1:426605/1‑251 (MQ=255)
                                                                                                                                                                                                              cgcgacgcaagTTGAAAAAGCCTGCCGATCCCGAAAGTTCGGGGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCTGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTAt                                    <  2:233709/241‑1 (MQ=255)
                                                                                                                                                                                                                    gcaCGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGAt                              >  1:429031/1‑251 (MQ=255)
                                                                                                                                                                                                                     caCGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGCCGTTTTGATCGGAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCTGCGCCGTTTTGCTGGGGATTGCACTGGTGATTGGCGCAGTCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGAtt                             <  1:295852/251‑1 (MQ=255)
                                                                                                                                                                                                                           gCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTc                        <  1:30617/250‑1 (MQ=255)
                                                                                                                                                                                                                                     gCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCCGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATa              >  1:77052/1‑250 (MQ=255)
                                                                                                                                                                                                                                       cGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGGTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTAGCTATGCCGATTCCGTCACTGACACTAAGc           >  2:196803/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aTCCCGAAAGTTCGGTGCTGTATGTTGCCTCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGGGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGGGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGc           <  2:116988/249‑1 (MQ=255)
                                                                                                                                                                                                                                         aTCCCGAAAGTTCGGTGATGTCTGTTGCCGAGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTACGTCACTGACAATAAGccc         >  1:67139/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aTCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAAc                                                               <  1:261896/197‑1 (MQ=255)
                                                                                                                                                                                                                                               aaaGTTCGGTGTTGTATGTTGCCGCGCCGATTATGGACGGCTCGAGGCTTATTGGCGTTTTGAGCGGAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCct   <  2:529102/251‑1 (MQ=255)
                                                                                                                                                                                                                                                aaGTTCGGAGATGTATGTTGCCGCGCCGATTATGGACGACTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCctc  >  1:174373/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
GGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTC  >  NC_000913/4636888‑4637377

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