Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,890,401 T→C 7.0% E106E (GAA→GAG)  hpd ← 4‑hydroxyphenylpyruvate dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,890,4010TC 7.0% 164.8 / 2.5 57E106E (GAA→GAG) hpd4‑hydroxyphenylpyruvate dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (28/25);  new base C (2/2);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CGGCCCCTACCGGGTTGCGATCCACGCCTTCGATGAAGTTGAAGTCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCAGCGGCGAAATACGAAGCGATGCTCTTGGGTTCATTGTTCA  >  NC_002947/3890253‑3890543
                                                                                                                                                     |                                                                                                                                              
cggccCCTACCGGGTTGCGATCCACGCCTTCGATGAAGTTGAAGTCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGtt                                                                                                                                                >  7:121635/1‑149 (MQ=255)
 ggccCCTACCGGGTTGCGATCCACGCCTTCGATGAAGTTGAAGTCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGttt                                                                                                                                               <  7:192987/149‑1 (MQ=255)
    cccTACCGGGTTGCGATCCACGCCTTCGATGAAGTTGAAGTCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGa                                                                                                                                            >  5:172497/1‑149 (MQ=255)
     ccTACCGGGTTGCGATCCACGCCTTCGATGAAGTTGAAGTCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGAt                                                                                                                                           <  1:249283/149‑1 (MQ=255)
        aCCGGGTTGCGATCCACGCCTTCGATGAAGTTGAAGTCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTc                                                                                                                                        >  8:66746/1‑149 (MQ=255)
         ccGGGTTGCGATCCACGCCTTCGATGAAGTTGAAGTCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTc                                                                                                                                        >  2:308947/1‑148 (MQ=255)
                 cGATCCACGCCTTCGATGAAGTTGTAGTCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCt                                                                                                                                <  5:289428/148‑1 (MQ=255)
                 cGATCCACGCCTTCGATGAAGTTGAAGTCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTCTCGAt                                                                                                                                           >  4:16850/1‑137 (MQ=255)
                 cGATCCACGCCTTCGATGAAGTTGAAGTCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTCTCGAt                                                                                                                                           <  3:16850/137‑1 (MQ=255)
                    tCCACGCCTTCGATGAAGTTGAAGTCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAgc                                                                                                                            <  5:224780/149‑1 (MQ=255)
                       aCGCCTTCGATGAAGTTGAAGTCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGcc                                                                                                                         >  7:160665/1‑149 (MQ=255)
                       aCGCCTTCGATGAAGTTGAAGTCGATGTCGTAGACAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGa                                                                                                                              >  7:267789/1‑144 (MQ=255)
                       aCGCCTTCGATGAAGTTGAAGTCGATGTCGTAGACAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGa                                                                                                                              <  8:267789/144‑1 (MQ=255)
                            ttCGATGAAGTTGAAGTCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGtt                                                                                                                    >  3:303961/1‑149 (MQ=255)
                            ttCGATGAAGTTGAAGTCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGtt                                                                                                                    >  8:304710/1‑149 (MQ=255)
                                 tgaagttgaagtCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAgg                                                                                                               <  8:121635/149‑1 (MQ=255)
                                 tgaagttgaagtCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAgg                                                                                                               <  4:303961/149‑1 (MQ=255)
                                    agttgaagtCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCGACAGGCTGAGCGCCCATTTCCAGGgcg                                                                                                            >  8:221622/1‑149 (MQ=255)
                                           gtCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCa                                                                                                      >  1:200066/1‑148 (MQ=255)
                                            tCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGc                                                                                                                                 >  6:207639/1‑120 (MQ=255)
                                            tCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGc                                                                                                                                 <  5:207639/120‑1 (MQ=255)
                                              gATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATaa                                                                                                   <  2:123581/148‑1 (MQ=255)
                                                     tagatagaGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATg                                                                                           >  7:140326/1‑149 (MQ=255)
                                                       gatagaGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATgg                                                                                          <  6:172497/148‑1 (MQ=255)
                                                              cTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATaa                                                                                                   >  3:99504/1‑132 (MQ=255)
                                                              cTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATaa                                                                                                   <  4:99504/132‑1 (MQ=255)
                                                                  cTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCAcg                                                                              >  1:14459/1‑149 (MQ=255)
                                                                   ttCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCt                                                                                                                                <  4:43013/98‑1 (MQ=255)
                                                                   ttCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCt                                                                                                                                >  3:43013/1‑98 (MQ=255)
                                                                   ttCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCAcgc                                                                             >  8:34463/1‑149 (MQ=255)
                                                                       tCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCCCGCGaa                                                                          >  4:217259/1‑148 (MQ=255)
                                                                         gAAGCGGTCGATCAGATAGAGCGGTGCCCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGaa                                                                          >  1:203058/1‑147 (MQ=255)
                                                                         gAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAg                                                                       <  3:68627/149‑1 (MQ=255)
                                                                                 cGATCAGATAGCGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATTCca                                                               <  7:304710/149‑1 (MQ=255)
                                                                                    tCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCcaca                                                             >  6:164731/1‑148 (MQ=255)
                                                                                      aGATAGAGCGGTGCCCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATa                                                           <  2:203058/149‑1 (MQ=255)
                                                                                       gATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATAcc                                                         <  5:164731/149‑1 (MQ=255)
                                                                                            agCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACgg                                                    <  3:217259/149‑1 (MQ=255)
                                                                                            agCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACgg                                                    >  1:73944/1‑149 (MQ=255)
                                                                                                    cccccccGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCgtg                                              <  8:160665/148‑1 (MQ=255)
                                                                                                     ccccccGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCgtgt                                             >  4:238535/1‑148 (MQ=255)
                                                                                                      cccccGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGtt                                            >  7:241951/1‑148 (MQ=255)
                                                                                                      cccccGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGtt                                            >  5:123378/1‑148 (MQ=255)
                                                                                                             tACACTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCAgcggcg                                    <  8:140326/149‑1 (MQ=255)
                                                                                                                      tCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCAGCGGCGAAATACg                             >  8:254907/1‑147 (MQ=255)
                                                                                                                            gCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCAGCGGCGAAATACGAAGCGATg                     >  6:136337/1‑149 (MQ=255)
                                                                                                                              aGGCGCAATTCCATCGGGCCGGTCTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCAGCGGCGAAATACGAAGCGATg                     <  5:117638/147‑1 (MQ=255)
                                                                                                                              aGGCGCAATTCCATCGGGCCGGTCTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCAGCGGCGAAATACGAAGCGATg                     >  6:117638/1‑147 (MQ=255)
                                                                                                                                gcgcAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGa                                                                           <  5:125151/91‑1 (MQ=255)
                                                                                                                                gcgcAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGa                                                                           >  6:125151/1‑91 (MQ=255)
                                                                                                                                 cgcAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGTCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCAGCGGCGAAATACGAAGCGATGCTCtt                <  7:221622/149‑1 (MQ=255)
                                                                                                                                    aaTTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCAGCGGCGAAATACGAAGCGATGCTCTTggg             <  1:308947/149‑1 (MQ=255)
                                                                                                                                     aTTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCAGCGGCGAAATACGAAGCGATg                     <  4:256620/140‑1 (MQ=255)
                                                                                                                                     aTTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCAGCGGCGAAATACGAAGCGATg                     >  3:256620/1‑140 (MQ=255)
                                                                                                                                         cATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCAGCGGCGAAATACGAAGCGATGCTCTTGGGTTCAt        <  7:34463/149‑1 (MQ=255)
                                                                                                                                          aTCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCAGCGGCGAAATACGAAGCGATGCTCTCGGGTTCAtt       <  5:136337/149‑1 (MQ=255)
                                                                                                                                              ggCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCAGCGGCGAAATACGAAGCGATGCTCTTGGGTTCATTGtt    >  1:4731/1‑148 (MQ=255)
                                                                                                                                               gCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCAGCGGCGAAATACGAAGCGATGCTCTTGGGTTCATTGtt    >  1:318723/1‑147 (MQ=255)
                                                                                                                                               gCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCAGCGGCGAAATACGAAGCGATGCTCTTGGGTTCATTGTTCa  <  2:73944/149‑1 (MQ=255)
                                                                                                                                                    tttCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCa                                          <  6:107885/104‑1 (MQ=255)
                                                                                                                                                    tttCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCa                                          >  5:107885/1‑104 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
CGGCCCCTACCGGGTTGCGATCCACGCCTTCGATGAAGTTGAAGTCGATGTCGTAGATAGAGCTGCCTTCCTCGAAGCGGTCGATCAGATAGAGCGGTG‑CCCCCCCGATACCCTTGATCGCCGGCAGGCGCAATTCCATCGGGCCGGTTTCGATTTCTACCGGCTGAGCGCCCAGTTCCAGGGCGCGGGCATAAGCCTCATGGGCATTGCGCACGCGAAAAGCCATGCCACATACCGACGGGCCGTGTTCAGCGGCGAAATACGAAGCGATGCTCTTGGGTTCATTGTTCA  >  NC_002947/3890253‑3890543

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