Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 88,991 C→T 8.2% L109L (TTG→TTA)  betC ← choline‑sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294788,9910CT 8.2% 174.3 / 4.7 61L109L (TTG→TTA) betCcholine‑sulfatase
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base T (3/2);  total (32/29)
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.37e-01

CGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGCCGCGTTGTCGTAGGCGCCAATGCGGCTGGGCAACT  >  NC_002947/88854‑89130
                                                                                                                                         |                                                                                                                                           
cgcacgcaCGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCg                                                                                                                                   <  2:152468/148‑1 (MQ=255)
   acgcacgGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGcc                                                                                                                               >  5:247407/1‑149 (MQ=255)
     gcacgGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCg                                                                                                                              >  4:75710/1‑148 (MQ=255)
     gcacgGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCg                                                                                                                              >  3:273961/1‑148 (MQ=255)
       acgGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCg                                                                                                                              >  1:100276/1‑146 (MQ=255)
               gCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCAt                                                                                                                   >  6:27741/1‑149 (MQ=255)
                       caCCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGccc                                                                                                            >  4:122818/1‑148 (MQ=255)
                                 aTGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAATCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCgcggt                                                                                                 >  4:12407/1‑149 (MQ=255)
                                  tgttgtGGTAACAGCTTGGACGCACATCGGGTTCGTCCAAGTTCACCGCCCAACCGAAGTCGGACCGGTAGATGTCGCTGCTCAGGCGTTCTTCTTACCAGTGTAACTGGTCCGGGCCGCTGAAGTGCATCTTGACCGACAGCgaggtg                                                                                                >  5:2568/1‑149 (MQ=255)
                                   gttgtGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCgcggtgc                                                                                               >  5:239838/1‑149 (MQ=255)
                                    ttgtGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCgcggtgcg                                                                                              >  2:236230/1‑149 (MQ=255)
                                      gtgGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCgcggtg                                                                                                <  8:78301/145‑1 (MQ=255)
                                      gtgGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCgcggtg                                                                                                >  7:78301/1‑145 (MQ=255)
                                         gTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCaa                                                                                                                                           <  3:161525/99‑1 (MQ=255)
                                         gTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCaa                                                                                                                                           >  4:157382/1‑99 (MQ=255)
                                         gTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCaa                                                                                                                                           >  4:161525/1‑99 (MQ=255)
                                         gTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCaa                                                                                                                                           <  3:157382/99‑1 (MQ=255)
                                          tACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGTGTAGATCTCGCTGGTCACCCTTTCTTAAAAAAAATACAACTGGTCTGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGcc                                                                                        >  4:146127/1‑149 (MQ=255)
                                           aCCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGcc                                                                                        >  3:180456/1‑148 (MQ=255)
                                            ccAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGccc                                                                                       <  2:100276/148‑1 (MQ=255)
                                              aGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGTAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAgg                                                                                    >  2:150927/1‑149 (MQ=255)
                                              aGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGTAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAgg                                                                                    >  2:146049/1‑149 (MQ=255)
                                              aGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAgg                                                                                    <  4:273961/149‑1 (MQ=255)
                                               gCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAgg                                                                                    <  1:120229/148‑1 (MQ=255)
                                               gCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCATGCGCGTCTGCGGCGTGCCCTAGCCCAgg                                                                                    >  2:120229/1‑148 (MQ=255)
                                                      cGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAg                                                                            <  4:135530/149‑1 (MQ=255)
                                                       gCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAgg                                                                           >  5:273964/1‑149 (MQ=255)
                                                        cacaTCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGt                                                                          <  6:273964/149‑1 (MQ=255)
                                                         acaTCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGt                                                                          >  4:91710/1‑148 (MQ=255)
                                                               ggTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGc                                                                   <  3:281425/149‑1 (MQ=255)
                                                                    gTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCAGGTCAGGCGTTCTTCATAGCCGTGCAACTGCTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGCTGCGGTAGCCCAGGCGACGCAGGTAATGGCCGTAgg                                                              >  3:222391/1‑149 (MQ=255)
                                                                       ccAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCg                                                           >  6:310364/1‑149 (MQ=255)
                                                                        cAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCgg                                                          <  6:239838/149‑1 (MQ=255)
                                                                         aGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCgg                                                          >  8:98906/1‑148 (MQ=255)
                                                                         aGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGc                                                         <  5:27741/149‑1 (MQ=255)
                                                                          gTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGca                                                        <  3:122818/149‑1 (MQ=255)
                                                                             cACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAgg                                                                                    >  8:238432/1‑118 (MQ=255)
                                                                             cACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAgg                                                                                    <  7:238432/118‑1 (MQ=255)
                                                                                cggccATCCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCATATc                                                    <  4:222391/144‑1 (MQ=255)
                                                                                cGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCgg                                                  >  2:229594/1‑149 (MQ=255)
                                                                                  cccAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGGGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCa                                                <  3:204405/149‑1 (MQ=255)
                                                                                       ccGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAggg                                             >  3:74884/1‑147 (MQ=255)
                                                                                       ccGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGaa                                           <  5:310364/149‑1 (MQ=255)
                                                                                        cGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAg                                          >  5:112086/1‑149 (MQ=255)
                                                                                            gTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTc                                        <  3:75710/147‑1 (MQ=255)
                                                                                              cGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGTAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGcc                                    <  1:150927/149‑1 (MQ=255)
                                                                                               ggCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGTAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGcc                                    <  1:146049/148‑1 (MQ=255)
                                                                                                 ccGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGTCGACGCAGGTAATGGGCGTAGGTCCGCACATCGGCAGGGAAGTCGGCcgag                                 >  3:18740/1‑147 (MQ=255)
                                                                                                 ccGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGCcgcg                                 <  4:74884/149‑1 (MQ=255)
                                                                                                   gggTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGCCGCGtt                               <  6:247407/149‑1 (MQ=255)
                                                                                                        gATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGCCGCGTTGTCGt                          <  1:236230/149‑1 (MQ=255)
                                                                                                          tGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGCCGCGTTGTCGTAg                        >  4:193108/1‑149 (MQ=255)
                                                                                                              gCTGGTCAGGCGTTCTTGATAGCCGTGCAAGTGGTCCTGGCCGCAGAAGTGCATCTTGCCCGACAGCGGGGTGCGGTAGCGCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGCCGCGTTGTCGTAGgcgc                    <  4:18740/149‑1 (MQ=255)
                                                                                                                      ggCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGCCGCGTTGTCGTAGGCGCCAATGCgg            <  3:91710/149‑1 (MQ=255)
                                                                                                                        cGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGCCGCGTTGTCGTAGGCGCCAATGCGGCt          <  3:193108/149‑1 (MQ=255)
                                                                                                                            cttcATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGCCGCGTTGTCGTAGGCGCCAATGCGGCTggg       <  3:12407/148‑1 (MQ=255)
                                                                                                                             ttcATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCTTCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGCCGCGTTGTCGTAGGCGCCAATGCGGCTGGGc      >  3:199179/1‑148 (MQ=255)
                                                                                                                             ttcATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGCCGCGTTGTCGTAGGCGCCAATGCGGCTGGGc      >  4:113273/1‑148 (MQ=255)
                                                                                                                               cATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCTTCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGCCGCGTTGTCGTAGGCGCCAATGCGGCTGGGc      <  4:199179/146‑1 (MQ=255)
                                                                                                                                aTAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGCCGCGTTGTCGTAGGCGCCAATGCGGCTGGGCAAc   <  2:116028/148‑1 (MQ=255)
                                                                                                                                aTAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGCCGCGTTGTCGTAGGCGCCAATGCGGCTGGGCAACt  >  1:116028/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                           
CGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGTTCACCGCCCAACCGTAGTCGGCCGGGTAGATGTCGCTGGTCAGGCGTTCTTCATAGCCGTGCAACTGGTCCGGGCCGCAGAAGTGCATCTTGCCCGACAGCGCGGTGCGGTAGCCCAGGCGACGCAGGTAATGGGCGTAGGTCGGCACATCGGCAGGGAAGTCGGCCGCGTTGTCGTAGGCGCCAATGCGGCTGGGCAACT  >  NC_002947/88854‑89130

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