Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 26,648 C→T H124H (CAC→CAT ispH → 4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091326,6480CT93.2% 204.4 / 3.9 74H124H (CAC→CATispH4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base T (36/32);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

ACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAGGAAG  >  NC_000913/26401‑26892
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
aCCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCg                                                                                                                                                                                                                                                   <  1:184093/251‑1 (MQ=255)
 ccGCTATGTGGTCGCTAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCgg                                                                                                                                                                                                                                                  >  1:390823/1‑251 (MQ=255)
      aTGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTg                                                                                                                                                                                                                                             <  1:346815/251‑1 (MQ=255)
       tgtgGTCGATAGCTTTCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCAATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTGCGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGCTGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGCAGAATCTATTCTCCTCGGTCACGCCGGGCACCCGGAAGTgg                                                                                                                                                                                                                                            >  1:221344/1‑251 (MQ=255)
           gTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGACGTGGAAgg                                                                                                                                                                                                                                        >  1:138882/1‑251 (MQ=255)
              gATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAggtac                                                                                                                                                                                                                                     <  1:284288/251‑4 (MQ=255)
               aTAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCAGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCGAGTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACa                                                                                                                                                                                                                                    <  1:143586/251‑1 (MQ=255)
               aTAGCTTGCGTGAGAGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAggtacg                                                                                                                                                                                                                                    >  1:110585/1‑247 (MQ=255)
                tAGCTTGCGTGAACGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGGTCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCCGACGTCGTGGCGAAGAGTCGGTTCTCATAGGCCACGCCGGGCCTCCGGGAGTGGCAggta                                                                                                                                                                                                                                      >  2:77843/1‑246 (MQ=255)
                     ggCGTGAGCGTGGGGCTTTCTTTATTGAGCAGATTAGCGAAGTCCCGGACGGCGCGATCCTGATTTTCTCCGCACACGGGGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGc                                                                                                                                                                                                                              <  2:186182/250‑1 (MQ=255)
                              gTGGGGCTATATTTATTGAGCATATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCGGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCcagtacagt                                                                                                                                                                                                                     <  1:179530/251‑1 (MQ=255)
                                        cTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGa                                                                                                                                                                                                             >  1:429125/1‑249 (MQ=255)
                                            aTTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAgggg                                                                                                                                                                                                        <  1:56036/250‑1 (MQ=255)
                                              tGAGCAGATGAGGGAGGTACCGGTCGGCGGGAGCCTTAGTGTCTCCGCACACGGTGTTTCTCAGGCGGGACGTAACGAAGCAAAAGGCCGCGATTTGTCGGTGTTTGATGCCACCTGTCCGCTGGGGACCAAAGTGCATATGGAAGTGGCCCGCGCCAGCGGGCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCGTCCGGAAGTGGGAGGTGCGATGGGCCAGTACAGTAACCCGGAAggggg                                                                                                                                                                                                       <  2:423598/249‑1 (MQ=255)
                                                        agcgaatgaCCGGACGGGTCGATCCGGATTTGCGCCGCACACGGGTTTTCTCTGGCGGTACGTAACGAAGCAAAAAGCCGCGAGTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAATTGCTTATGGAAGTCGCCCGCTCCAGCCGTCGTGGCGAAGAGTCGATTCTCAGCGGCCTCGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGTGGAAGGTATCTGGt                                                                                                                                                                                           <  2:207404/243‑1 (MQ=255)
                                                        aGCGAAGTACCGGACGGCCCGATCCTGATTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGt                                                                                                                                                                                           <  1:403460/250‑1 (MQ=255)
                                                         gCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTc                                                                                                                                                                                          <  1:192114/251‑1 (MQ=255)
                                                             attACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTGCTCAGGCGGTACGTAACAAAGCAAAAAGCCGCGATTTGACGTTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACATTAACCCGGAAGGGAGAATTTATCTGGTCGAAt                                                                                                                                                                                      <  2:211996/249‑1 (MQ=255)
                                                             aGTACCGGACGGCGCGATCCTGATTTTCTCCACACACGGTGTTTCTCAGGCGGAACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCACGCGCCAGCAGTCGTGGCGACGAGTCGGTTCTAATCGGCCACGCCGGGCATACGGAAGAGGAAGGTACGATGGGCCAGGACAGTAACCCGGAAGGGGGCATGTATCTGGTCGAAt                                                                                                                                                                                      >  1:400064/1‑251 (MQ=255)
                                                                     aCGGCGCGATCCTGATTTCCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCCCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGAc                                                                                                                                                                              <  2:290766/251‑1 (MQ=255)
                                                                          gcgatctggATTTTCGACGCACACGGTGTTTCGAAGGCGGTCCGTAACTAAGCAAAAAGCCGCGAGTTTACGGTGTTTGATGCCACCTGTCCGCGGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGCACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGAtgt                                                                                                                                                                         <  2:238443/243‑1 (MQ=255)
                                                                                         cccGCACACGGGGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTTCGCTGGTTACCAAAGGGCATATGGAAGTCGCCCGCGCCAGCCGGCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGt                                                                                                                                                          <  2:247639/250‑1 (MQ=255)
                                                                                          ccGCACACGGTGTTTCTCAGGCGGTACTTAACGAAGCAAAAAGCCGCGATTTGACGGTGGTTGATGCCTCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGTCGGACGATGTGTGGAAACTGACGGTc                                                                                                                                                         <  2:334994/251‑1 (MQ=255)
                                                                                              acacGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCCTCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGGGGAAACTGACGGTCaaa                                                                                                                                                      >  2:452763/1‑250 (MQ=255)
                                                                                                 cGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACg                                                                                                                                                  <  1:346886/251‑1 (MQ=255)
                                                                                                   ggttttCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCCCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGaa                                                                                                                                                <  2:349340/248‑1 (MQ=255)
                                                                                                         ctcAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGGCGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTAGGTTGGGCCAGTACAGTTACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGCTGTGTGGGAACTGACGGTCCAAAACGAAgaga                                                                                                                                            >  2:279977/1‑249 (MQ=255)
                                                                                                             ggCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCAGAGCATCCGGAAGTGGAAGGTACGATTGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGctcc                                                                                                                                      >  1:156564/1‑250 (MQ=255)
                                                                                                                ggTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGGGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGGGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCt                                                                                                                                   >  2:233290/1‑251 (MQ=255)
                                                                                                                  ggcGTAACGAAGCAAAAAGCCGGGATTGGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGCAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTATCCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCttt                                                                                                                                 <  2:449416/249‑1 (MQ=255)
                                                                                                                   aCGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCttt                                                                                                                                 <  1:418132/250‑1 (MQ=255)
                                                                                                                    cgaaaCGAAGCAAAAAGCCGCGATGTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGTATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATGATCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTTGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTAt                                                                                                                               <  2:156564/248‑1 (MQ=255)
                                                                                                                    cGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTAt                                                                                                                               <  2:152340/251‑1 (MQ=255)
                                                                                                                       aaCGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACt                                                                                                                                                                >  1:372476/1‑215 (MQ=255)
                                                                                                                       aaCGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACt                                                                                                                                                                <  2:372476/215‑1 (MQ=255)
                                                                                                                             gCAAAAATCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAACGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGGACGATGGGACAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCACTCCTTTATGACCCAGa                                                                                                                       >  2:422743/1‑250 (MQ=255)
                                                                                                                                 aaaGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGACGAGCAGCTCTCCTTTATGACCCAGAccccg                                                                                                                  >  2:138097/1‑248 (MQ=255)
                                                                                                                                    gccgcgATTTGACGGGGTTTGAGGCCCCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTg                                                                                                               <  2:340042/249‑1 (MQ=255)
                                                                                                                                        cgATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCgg                                                                                                           >  1:331912/1‑251 (MQ=255)
                                                                                                                                          aTTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTg                                                                                                         >  1:294900/1‑251 (MQ=255)
                                                                                                                                                  gtgtTTGATGCAACCTGTCCGCTGGTGACCCAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGGGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGGTGTGTGGAAACTTACTGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGagtc                                                                                                 >  2:2683/1‑248 (MQ=255)
                                                                                                                                                          tGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGAtg                                                                                         >  2:236040/1‑251 (MQ=255)
                                                                                                                                                              cccTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGAt                                                                                     <  2:221344/250‑1 (MQ=255)
                                                                                                                                                               ccTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGACGGGGGAATGTATCTGGTCGAATCGCCGGACGCTGTGTGGAAACTGACGGGCAAAAACGAAGCGAAGCTCTCCTTTATGACCCAGACCACGCAGTCGGCGGATGACACGACTGATGTGATc                                                                                    >  2:222308/1‑251 (MQ=255)
                                                                                                                                                               ccTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATc                                                                                    >  1:395063/1‑251 (MQ=255)
                                                                                                                                                                  gTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAgaga                                                                                                                                            >  1:457606/1‑192 (MQ=255)
                                                                                                                                                                         ggTGACCAAATTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCAAAGAGTCGATTCGCAGCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCt                                                                            <  1:232211/249‑1 (MQ=255)
                                                                                                                                                                           tGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGGATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAACAACGAAGAGAAGCTCTCCTTTATGGCCCAGACCAAGCTGTCGGTGGATGCCAAGTCTGATGTGATCGACTCGCTGCGt                                                                        >  2:119095/1‑251 (MQ=255)
                                                                                                                                                                                 aaGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACCGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGCCGGTGAGGGGAAACTGCCGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGCCCACGCTGCCGGGGGATGACCAGTCTTATGTGATCgcgtcgacgcgtaaaccc                                                                  >  2:275136/1‑234 (MQ=255)
                                                                                                                                                                                 aaGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGGACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAACACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCGGCTGTGTTCGACGCGCTGCGTAAACGc                                                                  >  2:73447/1‑251 (MQ=255)
                                                                                                                                                                                     gCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTcc                                                              <  1:236040/251‑1 (MQ=255)
                                                                                                                                                                                      cATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTAGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGAGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTccc                                                             >  1:184557/1‑251 (MQ=255)
                                                                                                                                                                                        tatGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCg                                                            >  1:261290/1‑250 (MQ=255)
                                                                                                                                                                                        tatGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCAGGGCATCCGGGAGTGGCAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACCCGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCg                                                            >  1:198622/1‑250 (MQ=255)
                                                                                                                                                                                             aaGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCCCGCTGTCGGTGGATGACACGTCTGATGTGATCGACTCGCTGCGTAAACGCTTCCCGAAAAtt                                                      <  1:138097/251‑1 (MQ=255)
                                                                                                                                                                                                  ccccGCGCCAGTCGTCGTGGGTAAGAGTAGATTCTCATCGGCCACGCCGGGCATCCGGAAGATGAAGGTACGATGGGCCTGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTTCTTTATGACCCCGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCgt                                                 <  2:33816/250‑2 (MQ=255)
                                                                                                                                                                                                   cccGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACCCGTCTGATGTGGTCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCggg                                                >  2:483825/1‑250 (MQ=255)
                                                                                                                                                                                                         ccagccgtcGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGCATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCgcgc                                          >  2:170254/9‑251 (MQ=255)
                                                                                                                                                                                                            gccgtcGTGGCGAAGAGGCGATTCGCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCaaa                                       <  2:11148/246‑1 (MQ=255)
                                                                                                                                                                                                              cgtcGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGACGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGGGGATGACACGTCTGAAGTGATCGACGCACTGCGTAAACGATTCCCGAAAATTGTCGGTCCGCGCacagg                                     >  1:325957/4‑247 (MQ=255)
                                                                                                                                                                                                                gcGTGGCGAAGAGGCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGGAACCCGGAAGGGGGAATGTATCTGGGCGAAGCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTTTCTGTGGATGACTCGTCTGATGTGATCGACGCGCTGCGGAAACGCTTCCCGTAAATTGTCGGTCCGCGCAAAGATg                                   <  2:235848/250‑1 (MQ=255)
                                                                                                                                                                                                                  gTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCCAAAATGAc                                 >  2:408050/1‑251 (MQ=255)
                                                                                                                                                                                                                  gTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGAc                                 >  1:16220/1‑251 (MQ=255)
                                                                                                                                                                                                                     gCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCTCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCGTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATc                              <  2:313940/251‑1 (MQ=255)
                                                                                                                                                                                                                      cGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGACAATTGTCGGTCCGCGCAAAGATGACATCt                             >  1:36467/1‑251 (MQ=255)
                                                                                                                                                                                                                       gaagagTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTg                            <  1:452763/251‑1 (MQ=255)
                                                                                                                                                                                                                       gaagagTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTg                            >  1:41687/1‑251 (MQ=255)
                                                                                                                                                                                                                       gaagagTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGACCGATCGGCCAGTACAGAAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGCGGAAACTGACGGTCAAAAACGAAGAGACGCTCTCCTATATGACCCAGACCACGCTGTCGGTGGATTACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCAGCGCAAAGATGACAtgtg                            >  1:198095/1‑248 (MQ=255)
                                                                                                                                                                                                                                aTTCTCATCGGCCACGCCGTGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCAc                   <  1:272539/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cATCGGCCACGCCGGGCCTCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTa              <  2:261290/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aTCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGCCCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACAGCCGCCACGCCACGACTa              >  2:10230/1‑250 (MQ=255)
                                                                                                                                                                                                                                        cgggcACGCCGGGCTGCCGGAAGTGGAAGGTACGATGGGCCAGTACAGGAACCCGGAAGGGGGAATGTATCTGGGCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGTGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAAcc           <  2:138882/247‑1 (MQ=255)
                                                                                                                                                                                                                                              cGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAgg     <  2:452882/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 cGGGCATCCGGAAGTGGAAGGTACGGTGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGCTGTGTGGAAACTGACGGGCAAAAACGAAGAGAAGCTCTCCTTTATGACCCCGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTAAGCCACGACTAACCGTccgggtg  >  2:247744/1‑245 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
ACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAGGAAG  >  NC_000913/26401‑26892

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: