Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 86,240 G→T 7.5% L160M (CTG→ATG)  PP_0075 ← choline sulfate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294786,2400GT 7.5% 210.8 / 6.3 67L160M (CTG→ATG) PP_0075choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base G (27/35);  new base T (2/3);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.12e-01

AACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCG  >  NC_002947/86092‑86385
                                                                                                                                                    |                                                                                                                                                 
aaCAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAg                                                                                                                                                   >  7:24770/1‑149 (MQ=255)
 aCAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAgg                                                                                                                                                  >  3:31722/1‑149 (MQ=255)
     gCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCa                                                                                                                                              >  4:160210/1‑149 (MQ=255)
             aGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCATGGCCCTGCCcagc                                                                                                                                      <  7:3176/149‑1 (MQ=255)
             aGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCATGGCCCTGCCcagc                                                                                                                                      <  7:7760/149‑1 (MQ=255)
             aGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCcagc                                                                                                                                      <  1:138171/149‑1 (MQ=255)
              gCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCATGGCCCTGCCcagcc                                                                                                                                     <  2:393/149‑2 (MQ=255)
                gggggCGCCAGCGTTTGACGGCCGCACCGCTCACGACGGGGGCCCGGCCCAGCCGCAAGGCAGGGAGGTCGACCTAGCCTGCGTGGCGGGCCAGGTCCTGCACGGGCATCAGCGCGGGGGCCTGGCGGGGCAGGGCCATGCCCAGCAgg                                                                                                                                   <  2:116103/147‑1 (MQ=255)
                   ggCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGtt                                                                                                                                 <  3:261197/148‑1 (MQ=255)
                   ggCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTg                                                                                                                                <  3:113144/149‑1 (MQ=255)
                     cgcCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTggg                                                                                                                              <  5:285524/149‑1 (MQ=255)
                      gcCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTggg                                                                                                                              >  7:58488/1‑148 (MQ=255)
                      gcCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGt                                                                                                                             <  8:26622/149‑1 (MQ=255)
                      gcCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCAAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGCTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGATgggc                                                                                                                             >  7:39128/1‑148 (MQ=255)
                       ccAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTggg                                                                                                                              >  4:180685/1‑147 (MQ=255)
                                gAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGccc                                                                                                                    <  3:160210/148‑1 (MQ=255)
                                   ggCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAgg                                                                                                                <  8:89229/149‑1 (MQ=255)
                                    gCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAgg                                                                                                                >  7:200058/1‑148 (MQ=255)
                                       gCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAAt                                                                                                            <  8:200058/149‑1 (MQ=255)
                                        cACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAAt                                                                                                            >  6:162914/1‑148 (MQ=255)
                                        cACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATg                                                                                                           <  3:180685/149‑1 (MQ=255)
                                           cGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGAc                                                                                                         >  7:239431/1‑148 (MQ=255)
                                                 cgacgGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATggcg                                                                                                  <  8:91828/149‑1 (MQ=255)
                                                        ggCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCAcc                                                                                             >  3:241792/1‑147 (MQ=255)
                                                           cAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAgg                                                                                                                >  8:69070/1‑125 (MQ=255)
                                                               cccAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGaa                                                                                    <  4:31722/149‑1 (MQ=255)
                                                                 cagcagCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTTATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGc                                                                                  <  6:26438/149‑1 (MQ=255)
                                                                  agcagcAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGcc                                                                                 >  8:173211/1‑149 (MQ=255)
                                                                    cagcaAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGa                                                                               <  6:31474/149‑1 (MQ=255)
                                                                    cagcaAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCACGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGa                                                                               <  5:264170/149‑1 (MQ=255)
                                                                            aaggTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACg                                                                       <  1:181619/149‑1 (MQ=255)
                                                                               cgtgcaggACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGaa                                                                    >  7:286636/8‑149 (MQ=255)
                                                                                 aggtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTGGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATg                                                                  >  6:299293/1‑149 (MQ=255)
                                                                                   gtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATgg                                                                 <  8:239431/148‑1 (MQ=255)
                                                                                      gACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGcaccac                                                                        <  8:58322/138‑1 (MQ=255)
                                                                                      gACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGcaccac                                                                        >  7:58322/1‑138 (MQ=255)
                                                                                      gACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGgaga                                                             <  5:162914/149‑1 (MQ=255)
                                                                                       aCCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGAc                                                            >  5:55954/1‑149 (MQ=255)
                                                                                                 cGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAg                                                        <  6:144162/143‑1 (MQ=255)
                                                                                                 cGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAg                                                        >  5:144162/1‑143 (MQ=255)
                                                                                                 cGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGAc                                                   <  8:24770/148‑1 (MQ=255)
                                                                                                        ggcCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGccc                                                                                                                    <  4:62960/76‑1 (MQ=255)
                                                                                                        ggcCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGccc                                                                                                                    >  3:62960/1‑76 (MQ=255)
                                                                                                          ccAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCg                                         <  4:241792/149‑1 (MQ=255)
                                                                                                           cAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGAACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCg                                         >  4:65378/1‑148 (MQ=255)
                                                                                                              gTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAgc                                     >  2:275894/1‑149 (MQ=255)
                                                                                                                ccTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAgcgc                                   <  8:286636/149‑1 (MQ=255)
                                                                                                                  tGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCaa                                 >  2:100914/1‑149 (MQ=255)
                                                                                                                     aCGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAGGCGCAATAg                              >  2:202774/1‑149 (MQ=255)
                                                                                                                      cGGGCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGc                             <  7:173211/149‑1 (MQ=255)
                                                                                                                         tcatcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGccc                           >  2:18654/1‑148 (MQ=255)
                                                                                                                           atcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGag                        >  3:198946/1‑149 (MQ=255)
                                                                                                                           atcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGag                        <  6:55954/149‑1 (MQ=255)
                                                                                                                           atcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGag                        <  8:58488/149‑1 (MQ=255)
                                                                                                                            tcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGagc                       <  3:274881/149‑1 (MQ=255)
                                                                                                                               gcgcGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGagcagc                    <  5:204532/149‑1 (MQ=255)
                                                                                                                               gcgcGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGagcag                     >  3:133607/1‑148 (MQ=255)
                                                                                                                               gcgcGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCA‑‑GCGCCGAAGCGCAATAGCCCGagcagcag                  >  7:226133/1‑149 (MQ=255)
                                                                                                                                cgcgGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGagcagca                   <  4:106072/149‑1 (MQ=255)
                                                                                                                                  cgGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCTCAATAGCCCGagcagcagc                 <  4:133607/149‑1 (MQ=255)
                                                                                                                                    gTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGACGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATg                                                                  <  2:269733/98‑1 (MQ=255)
                                                                                                                                    gTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGACGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATg                                                                  >  1:269733/1‑98 (MQ=255)
                                                                                                                                       aacctgctggGCATGGCCATGCCCAGCCGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGATGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGa             >  5:4819/6‑148 (MQ=255)
                                                                                                                                          tggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAAt         >  7:13977/1‑149 (MQ=255)
                                                                                                                                            gctggGCATGGCCCTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATgc       >  3:51014/1‑149 (MQ=255)
                                                                                                                                            gctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATg        <  8:13977/148‑1 (MQ=255)
                                                                                                                                              tggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATgccg     >  8:218312/1‑149 (MQ=255)
                                                                                                                                                ggCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAGGCGCAATAGCCCGAGCAGCAGCTGGAAAATgccgc    <  1:202774/148‑1 (MQ=255)
                                                                                                                                                 gCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGccgccg  <  7:218312/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCG  >  NC_002947/86092‑86385

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