Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,424,497 C→G 5.9% H32Q (CAC→CAG PP_4764 → histone deacetylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,424,4970CG 5.9% 216.4 / 4.1 68H32Q (CAC→CAGPP_4764histone deacetylase family protein
Reads supporting (aligned to +/- strand):  ref base C (33/31);  new base G (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

CAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGCC  >  NC_002947/5424351‑5424645
                                                                                                                                                      |                                                                                                                                                      
cAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCt                                                                                                                                                      >  7:203384/1‑149 (MQ=255)
cAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCt                                                                                                                                                      >  7:176903/1‑149 (MQ=255)
      cccGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCgg                                                                                                                                               <  8:176903/149‑1 (MQ=255)
      cccGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCgg                                                                                                                                               <  8:203384/149‑1 (MQ=255)
      cccGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCgg                                                                                                                                               <  4:192756/149‑1 (MQ=255)
             tGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCggg                                                                                                                                        >  2:33969/1‑149 (MQ=255)
               gCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCAGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCt                                                                                                                                      >  6:319198/1‑149 (MQ=255)
                 aaGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGa                                                                                                                                    >  5:153944/1‑149 (MQ=255)
                 aaGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGa                                                                                                                                    >  5:148437/1‑149 (MQ=255)
                 aaGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCCGACCAGGTGGTCgtgc                                                                                                                                              >  7:61397/1‑136 (MQ=255)
                 aaGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCCGACCAGGTGGTCgtgc                                                                                                                                              <  8:61397/139‑4 (MQ=255)
                  aGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGac                                                                                                                                   >  7:28700/1‑149 (MQ=255)
                        gCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGa                                                                                                                             >  2:178136/1‑149 (MQ=255)
                           cTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGTTCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGAc                                                                                                                            <  1:90684/147‑1 (MQ=255)
                           cTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGTTCTACCACGAAGATTACAGCCCGGAGTTACCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGAc                                                                                                                            >  2:90684/1‑147 (MQ=255)
                           cTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCa                                                                                                                          <  8:28700/149‑1 (MQ=255)
                              gCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGc                                                                                                                       >  6:334697/1‑149 (MQ=255)
                              gCCCTGATTGCGACCCCATCCATCCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGACTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGc                                                                                                                       <  4:311036/149‑1 (MQ=255)
                                 cttaTTGCTACCCCATCCATGCCGTTGCCGCTCATCTACCACGAACATTCCCGCCCGTCCGTTCCCCGCTGACTCCAACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAgg                                                                                                                        <  1:112127/146‑1 (MQ=255)
                                    aTTGCGACCCCATCCATGCCGTTGCAGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGACGGCAGCGCGCTGACCACCGCCCAGGCCctgct                                                                                                                 >  4:196308/1‑149 (MQ=255)
                                     ttGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCctgct                                                                                                                 <  3:196308/148‑1 (MQ=255)
                                       gCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCtgctgc                                                                                                               <  5:162245/148‑1 (MQ=255)
                                          aCCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGccc                                                                                                           <  5:319198/149‑1 (MQ=255)
                                           ccccATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACAAGGCCCTGCTCCGCCCg                                                                                                          >  5:302023/1‑149 (MQ=255)
                                             ccatccatGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGcc                                                                                                            >  5:352734/1‑145 (MQ=255)
                                             ccatccatGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGa                                                                                                        >  4:255149/1‑149 (MQ=255)
                                             ccatccatGCCGTTGCCGATGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGcc                                                                                                            <  6:352734/145‑1 (MQ=255)
                                                tccatGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACa                                                                                                      >  6:245810/1‑148 (MQ=255)
                                                    tGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGc                                                                                                <  5:334697/149‑1 (MQ=255)
                                                    tGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGc                                                                                                <  1:33969/149‑1 (MQ=255)
                                                         ttGCCGCTGATCTACCACGAAGATTACAGTCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGaccacc                                                                                                                               >  1:210474/1‑114 (MQ=255)
                                                         ttGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGaccacc                                                                                                                               <  2:210474/114‑1 (MQ=255)
                                                            ccGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACg                                                                                         <  5:245810/148‑1 (MQ=255)
                                                                    cTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCg                                                                                >  7:323200/1‑149 (MQ=255)
                                                                       ccACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGccc                                                                             >  7:215530/1‑149 (MQ=255)
                                                                       ccACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGccc                                                                             >  2:89654/1‑149 (MQ=255)
                                                                         aCGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGcc                                                                                               <  8:102807/129‑1 (MQ=255)
                                                                         aCGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGcc                                                                                               >  7:102807/1‑129 (MQ=255)
                                                                         aCGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGc                                                                               >  5:334432/1‑145 (MQ=255)
                                                                         aCGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGc                                                                               <  6:334432/145‑1 (MQ=255)
                                                                             aGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGccc                                                                       <  2:298543/149‑1 (MQ=255)
                                                                              gATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCa                                                                      >  8:242317/1‑149 (MQ=255)
                                                                               aTTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCAt                                                                     <  8:323200/149‑1 (MQ=255)
                                                                                   cAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGAcc                                                                 >  4:353898/1‑149 (MQ=255)
                                                                                                ccccGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGAc                                                                                       >  8:10414/1‑115 (MQ=255)
                                                                                                ccccGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGAc                                                                                       <  7:10414/115‑1 (MQ=255)
                                                                                                ccccGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATc                                                      >  2:159960/1‑148 (MQ=255)
                                                                                                 cccGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCCGCTACATCgc                                                    >  3:183192/1‑148 (MQ=255)
                                                                                                   cGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAgc                                                  <  7:242317/149‑1 (MQ=255)
                                                                                                   cGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAgc                                                  <  8:215530/149‑1 (MQ=255)
                                                                                                    gCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGGTGCGCCCGGACATTCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAgc                                                  <  1:89654/149‑1 (MQ=255)
                                                                                                              acCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACAt                                           >  3:223276/1‑148 (MQ=255)
                                                                                                                cGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTgcgc                                                                                                             <  5:272290/81‑1 (MQ=255)
                                                                                                                cGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTgcgc                                                                                                             >  6:272290/1‑81 (MQ=255)
                                                                                                                      ccGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGa                                  <  1:178136/149‑1 (MQ=255)
                                                                                                                             aCAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcg                           >  8:50790/1‑149 (MQ=255)
                                                                                                                              cAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgc                          <  1:159960/149‑1 (MQ=255)
                                                                                                                               aaGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCAt                                                                     <  6:138347/105‑1 (MQ=255)
                                                                                                                               aaGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgcg                         >  3:80192/1‑149 (MQ=255)
                                                                                                                                      gCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGAc                  <  4:80192/149‑1 (MQ=255)
                                                                                                                                           gctgccGACCAGGTGGTCGTGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGa                                  <  7:7613/124‑1 (MQ=255)
                                                                                                                                           gctgccGACCAGGTGGTCGTGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGa                                  >  8:7613/6‑129 (MQ=255)
                                                                                                                                           ctgcACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTTCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAgc              <  6:153944/148‑1 (MQ=255)
                                                                                                                                           ctgcACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAgc              <  6:148437/148‑1 (MQ=255)
                                                                                                                                              cACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACgg                                     <  6:185845/122‑1 (MQ=255)
                                                                                                                                              cACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACgg                                     >  5:185845/1‑122 (MQ=255)
                                                                                                                                                cGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCgccgcc        >  2:2601/1‑149 (MQ=255)
                                                                                                                                                cGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCgccgcc        >  1:182944/1‑149 (MQ=255)
                                                                                                                                                      ccTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGcc  >  1:304681/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGA‑‑ACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGCC  >  NC_002947/5424351‑5424645

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