Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 26,606 T→C S110S (AGT→AGC 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,6060TC91.8% 179.2 / 1.8 74S110S (AGT→AGCispH4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase
Reads supporting (aligned to +/- strand):  ref base T (1/5);  new base C (36/31);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGT  >  NC_000913/26364‑26846
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
ggCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCATCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTgg                                                                                                                                                                                                                                          >  1:192375/1‑251 (MQ=255)
  cgcACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGCAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGACGTCGCCCGCGCCAGCCGTCGTGGCg                                                                                                                                                                                                                                        >  1:451091/1‑251 (MQ=255)
     aCCGATATATGTCCGTCACGAAGTGGTACATACCCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGCAGCAAAAAGCCGCGCTTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGCAGTCGCCCGCGCCAGCCGGCGGGGCga                                                                                                                                                                                                                                       >  2:416333/1‑249 (MQ=255)
         atatatGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCgcagaggc                                                                                                                                                                                                                                 >  1:38164/1‑248 (MQ=255)
           atatGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCgaagagtcg                                                                                                                                                                                                                                >  1:170975/1‑246 (MQ=255)
              tGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGAACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCAACCTGTCCGCTGGTGACCAAAGTGAATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCgacagatcgcttt                                                                                                                                                                                                                            >  2:174070/1‑240 (MQ=255)
               gTCCGGCACGAAGTGGTACGTAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATtc                                                                                                                                                                                                                            <  1:416971/250‑1 (MQ=255)
                 ccGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCa                                                                                                                                                                                                                         <  1:479364/251‑1 (MQ=255)
                  cGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCAt                                                                                                                                                                                                                        <  1:121061/251‑1 (MQ=255)
                   gTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATc                                                                                                                                                                                                                       >  1:235848/1‑251 (MQ=255)
                       cGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCggcc                                                                                                                                                                                                                   <  1:188297/251‑3 (MQ=255)
                            tGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGCTTCTCATCggccccgcc                                                                                                                                                                                                              >  2:346815/1‑244 (MQ=255)
                             ggTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGcc                                                                                                                                                                                                              <  1:82045/250‑1 (MQ=255)
                             ggTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCGGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGcc                                                                                                                                                                                                              <  2:479631/250‑1 (MQ=255)
                              gTACAGAACCGCTATGTGGTCTATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCCCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCTCGCCgg                                                                                                                                                                                                            <  2:142910/251‑1 (MQ=255)
                                   tAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGCAGCGTCGATTCTCCTCGGCCACGCCGGGCatc                                                                                                                                                                                                       >  2:192114/1‑249 (MQ=255)
                                     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)
                                                      aGCTTGCGTGAACGTGGGGCTATCTTTATTGAGCAGCTTAGCGAAGTAACGGACGGCGCGCTCCTGATTTTCTCTGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCACAAAGCCGCGCTTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGAGCATATGGATGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGACGATTCTCATCAGCCACTCCGGGcctcctgacgtggcaggaccgct                                                                                                                                                                                    >  2:166229/1‑229 (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                                                                                                                  <  2:372476/215‑1 (MQ=255)
                                                                                                                                                            aaCGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACt                                                                                                                  >  1:372476/1‑215 (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)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
GGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGT  >  NC_000913/26364‑26846

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: