Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 401,929 C→T 5.5% S199S (TCG→TCA)  PP_0335 ← glutathione S‑transferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947401,9290CT 5.5% 226.5 / 2.3 73S199S (TCG→TCA) PP_0335glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base C (38/31);  new base T (2/2);  total (40/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 = 9.25e-01

TGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGC  >  NC_002947/401783‑402070
                                                                                                                                                  |                                                                                                                                              
tGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGa                                                                                                                                              >  2:249823/1‑149 (MQ=255)
  gggTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGAt                                                                                                                                             <  2:62316/148‑1 (MQ=255)
       tCAGGCACCGGCTTCGCCGGTTTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGTACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGc                                                                                                                                       <  8:330399/149‑1 (MQ=255)
         gaacACCGGCGAAGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATg                                                                                                                                            >  5:113228/4‑142 (MQ=255)
         gaacACCGGCGAAGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATg                                                                                                                                            <  6:113228/139‑1 (MQ=255)
           gCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGa                                                                                                                                   >  4:315018/1‑149 (MQ=255)
           gCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGa                                                                                                                                   >  8:218783/1‑149 (MQ=255)
            cACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGa‑c                                                                                                                                 <  3:124181/149‑1 (MQ=255)
                   ttCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCg                                                                                                                          >  3:18278/1‑149 (MQ=255)
                      gCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCa                                                                                                                       >  7:79565/1‑149 (MQ=255)
                       ccGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCa                                                                                                                       >  8:185235/1‑148 (MQ=255)
                        cGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTc                                                                                                                           >  5:217022/1‑143 (MQ=255)
                        cGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTc                                                                                                                           <  6:217022/143‑1 (MQ=255)
                         ggTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATAc                                                                                                                    <  8:79565/149‑1 (MQ=255)
                            gtTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGCCCATCGATGGCCGCTTGA‑CCACCTCGGCATACCAc                                                                                                                 >  2:141494/1‑149 (MQ=255)
                             ggagcgGGCGTGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCAcg                                                                                                                >  4:65491/4‑149 (MQ=255)
                                           cGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGc                                                                                                  <  7:137105/149‑1 (MQ=255)
                                            gCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCt                                                                                                 >  5:32388/1‑149 (MQ=255)
                                              tCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTgg                                                                                               >  3:75789/1‑149 (MQ=255)
                                               cccAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGgc                                                                                              >  8:190494/1‑149 (MQ=255)
                                               cccAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGgc                                                                                              >  4:279486/1‑149 (MQ=255)
                                                ccAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGgcg                                                                                             <  3:202463/149‑1 (MQ=255)
                                                 cAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCCACCTCGGCATACCACGCCGCCAACCCAGGCTGgcg                                                                                             <  1:249823/149‑1 (MQ=255)
                                                  aaaaGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCt                                                                                           >  8:54763/1‑149 (MQ=255)
                                                  aaaaGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAAACCAGGCTGGCGCt                                                                                           >  1:286075/1‑149 (MQ=255)
                                                     aGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCAc                                                                                        <  2:20474/149‑1 (MQ=255)
                                                           ccacAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCTCTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACgccagc                                                                                  >  8:95198/1‑149 (MQ=255)
                                                             acaGATAGCGAACACTGTGGAACACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGTGAGCgccagcct                                                                                <  7:164784/149‑2 (MQ=255)
                                                              caGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCaa                                                                               <  7:190494/149‑1 (MQ=255)
                                                              caGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCaa                                                                               <  3:279486/149‑1 (MQ=255)
                                                                gATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGTGAGCgccagcctggg                                                                             <  7:185235/149‑5 (MQ=255)
                                                                    gCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCCAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCgg                                                                         <  5:323008/149‑1 (MQ=255)
                                                                     cGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCTGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATTGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCggg                                                                        >  5:25542/1‑149 (MQ=255)
                                                                     cGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCTGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATTGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCggg                                                                        >  5:30999/1‑149 (MQ=255)
                                                                     cGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCggg                                                                        <  4:75789/149‑1 (MQ=255)
                                                                     cGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCggg                                                                        >  4:96182/1‑149 (MQ=255)
                                                                     cGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCggg                                                                        <  6:294806/149‑1 (MQ=255)
                                                                      gAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCACCCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGc                                                                       >  5:342619/1‑149 (MQ=255)
                                                                              tgGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGCCGCa                                                                >  7:7359/1‑148 (MQ=255)
                                                                                gAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGa‑c                                                                                                                                 <  4:120211/81‑1 (MQ=255)
                                                                                gAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGa‑c                                                                                                                                 >  3:120211/1‑81 (MQ=255)
                                                                                gAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCaaa                                                              <  6:32388/148‑1 (MQ=255)
                                                                                gAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAAt                                                             <  4:18278/149‑1 (MQ=255)
                                                                                gAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAAt                                                             >  5:376982/1‑149 (MQ=255)
                                                                                  aaaCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATc                                                            >  7:351512/1‑148 (MQ=255)
                                                                                     ccTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGa                                                        <  1:238811/149‑1 (MQ=255)
                                                                                          tAGGAGCGGGTTTACCTGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATTGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGcc                                                   <  6:25542/149‑1 (MQ=255)
                                                                                          tAGGAGCGGGTTTACCTGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATTGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGcc                                                   <  6:30999/149‑1 (MQ=255)
                                                                                                     ttACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCTGGCTGACGCACATCCAGATAGCCCAGCGCACAgg                                        >  3:295458/1‑149 (MQ=255)
                                                                                                      tACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAgg                                        <  3:65491/148‑1 (MQ=255)
                                                                                                       aCCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGcc                                      >  8:252934/1‑149 (MQ=255)
                                                                                                         ccGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CAACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCAc                                    >  6:344039/1‑149 (MQ=255)
                                                                                                         ccGCGAAGCGGCCCTCAGGCCACTGCCGCTGTAGCGACCATCGATGGCCGCATGA‑CCACCTCGGCATACCACGCCGCAAACCCAGGCTAGCGCTCCCGCCAGCAAAACCCGCGCTCACGCAAATCCAGATAGCCAAGCGCACAGGCCAc                                    >  4:60769/1‑149 (MQ=255)
                                                                                                            cgAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGc                                  >  3:275511/1‑148 (MQ=255)
                                                                                                             gAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCa                                <  3:315018/149‑1 (MQ=255)
                                                                                                                ggggCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGC‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATc                             >  6:256036/1‑149 (MQ=255)
                                                                                                                ggggCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATc                             <  7:218783/149‑1 (MQ=255)
                                                                                                                   gCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATcgcc                          >  4:121627/1‑149 (MQ=255)
                                                                                                                     ccTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAgg                                        >  7:88500/1‑133 (MQ=255)
                                                                                                                     ccTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAgg                                        <  8:88500/133‑1 (MQ=255)
                                                                                                                      cTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATcgccgc                        >  1:52520/1‑148 (MQ=255)
                                                                                                                      cTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATcgccgc                        >  1:173069/1‑148 (MQ=255)
                                                                                                                      cTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCgccgcc                       >  8:213093/1‑149 (MQ=255)
                                                                                                                         aGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAgg                    <  3:60769/149‑1 (MQ=255)
                                                                                                                              aCTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGaa               <  2:166425/149‑1 (MQ=255)
                                                                                                                                tGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGc             >  7:188143/1‑149 (MQ=255)
                                                                                                                                tGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGc             >  4:77111/1‑149 (MQ=255)
                                                                                                                                   cgcTGTAGCGACCATCGCTGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCgg          <  4:39996/149‑1 (MQ=255)
                                                                                                                                   cgcTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCgg          <  7:252934/149‑1 (MQ=255)
                                                                                                                                    gcTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGga         <  2:173069/149‑1 (MQ=255)
                                                                                                                                        tAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGAt      >  7:16063/1‑148 (MQ=255)
                                                                                                                                          gCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAg   <  2:52520/149‑1 (MQ=255)
                                                                                                                                           cGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGc  >  2:50387/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
TGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA‑CCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGC  >  NC_002947/401783‑402070

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: ‑