Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,614,727 G→T 9.6% V200V (GTG→GTT) ‡ pet → protein Pet18

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,614,7270GT 9.6% 179.9 / 6.2 60V200V (GTG→GTT) ‡petprotein Pet18
Reads supporting (aligned to +/- strand):  ref base G (30/24);  new base T (2/4);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGCACGCT  >  NC_002947/3614584‑3614872
                                                                                                                                               |                                                                                                                                                 
ccGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCt                                                                                                                                              <  4:29803/149‑1 (MQ=255)
  ggACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTgg                                                                                                                                            <  7:91057/149‑1 (MQ=255)
  ggACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTg                                                                                                                                             >  4:188375/1‑148 (MQ=255)
   gACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCCACGGTGCGCTggt                                                                                                                                           <  2:165525/149‑2 (MQ=255)
   gACAACTTCGTCCATGCCGACTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTATCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCCACGGTGCGCTggt                                                                                                                                           <  3:126609/149‑2 (MQ=255)
    aCAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCCACGGTGCGCTggta                                                                                                                                          <  7:68826/149‑3 (MQ=255)
        cTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTgg                                                                                                                                      >  1:208959/1‑149 (MQ=255)
            gTCCATGCCGAGTGGATTACCCTGGATGACAACGCGGCATTGCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCttt                                                                                                                                  <  7:203991/149‑1 (MQ=255)
               cATGCCGAGTGGATTACCCTGCATGACACCCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTttctt                                                                                                                               >  3:48965/1‑149 (MQ=255)
                   ccGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGAt                                                                                                                           <  7:248597/149‑1 (MQ=255)
                     gAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGc                                                                                                                         >  6:112916/1‑149 (MQ=255)
                                ccTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCg                                                                                                               >  7:246867/1‑148 (MQ=255)
                                ccTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTACAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCg                                                                                                               >  1:260144/1‑148 (MQ=255)
                                 cTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGa                                                                                                             <  2:22125/149‑1 (MQ=255)
                                  tGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAg                                                                                                            >  4:211185/1‑149 (MQ=255)
                                   gCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAgg                                                                                                           >  5:196502/1‑149 (MQ=255)
                                    cATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAggg                                                                                                          <  2:140191/149‑1 (MQ=255)
                                    cATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAgg                                                                                                           <  6:138548/148‑1 (MQ=255)
                                    cATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAg                                                                                                            <  6:278267/147‑1 (MQ=255)
                                       gACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCa                                                                                                        >  5:817/1‑148 (MQ=255)
                                        aCAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCaaa                                                                                                      >  1:7275/1‑149 (MQ=255)
                                        aCAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCaaa                                                                                                      <  5:112916/149‑1 (MQ=255)
                                        aCAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCaa                                                                                                       <  2:208959/148‑1 (MQ=255)
                                          aaCCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCCACGGTGCGCTg                                                                                                                                             >  2:338707/1‑108 (MQ=255)
                                          aaCCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCCACGGTGCGCTg                                                                                                                                             <  1:338707/108‑1 (MQ=255)
                                                 cATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATgg                                                                                             <  8:246867/149‑1 (MQ=255)
                                                   ttCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGat                                                                                           <  3:211185/149‑1 (MQ=255)
                                                      cAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATc                                                                                        >  1:271/1‑149 (MQ=255)
                                                           aTTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTc                                                                                     >  3:60868/1‑147 (MQ=255)
                                                           aTTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTc                                                                                     >  4:141918/1‑147 (MQ=255)
                                                                      tCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGcc                                                                              <  3:167367/143‑1 (MQ=255)
                                                                      tCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGcc                                                                              >  4:167367/1‑143 (MQ=255)
                                                                        cTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGc                                                                       <  6:196502/148‑1 (MQ=255)
                                                                                cGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTgccac                                                              <  6:817/149‑1 (MQ=255)
                                                                                cGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTgccac                                                              >  6:48089/1‑149 (MQ=255)
                                                                                cGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCAGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTgccac                                                              >  2:360042/1‑149 (MQ=255)
                                                                                        aCCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAAt                                                      <  3:188375/149‑1 (MQ=255)
                                                                                            gttCGGCCTGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGaa                                                  >  8:34076/3‑149 (MQ=255)
                                                                                            cgTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGaa                                                  <  1:360042/149‑1 (MQ=255)
                                                                                                 gCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTgc                                                                 <  6:334111/129‑1 (MQ=255)
                                                                                                 gCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTgc                                                                 >  5:334111/1‑129 (MQ=255)
                                                                                                   ccGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTAc                                           >  6:84263/1‑149 (MQ=255)
                                                                                                   ccGCAGGAGGCCGCTTTGGCAAGTGACTTCTACAAGCGCACCGTGGCGCTGGAGCAGGCTTTCTTCGATGCAGTGAACCCGCAGGGCAAATGAGCATGTATATCATCGACCGCCTCAAGGCCGCTGCCACGCCACAATCGAACCGCTAc                                           >  4:269417/1‑149 (MQ=255)
                                                                                                    cGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACg                                          <  4:48965/149‑1 (MQ=255)
                                                                                                         gaggCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGAc                                     >  6:102993/1‑149 (MQ=255)
                                                                                                                   tGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCg                           >  4:201904/1‑149 (MQ=255)
                                                                                                                      cAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGc                        <  2:271/149‑1 (MQ=255)
                                                                                                                         gTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCgg                      >  3:173826/1‑148 (MQ=255)
                                                                                                                            aCTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCgg                      >  6:289046/1‑145 (MQ=255)
                                                                                                                             cTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAAAGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGAt                 <  2:7275/149‑1 (MQ=255)
                                                                                                                                gctccAGCGCCACGGTGCGCTGGTGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATggg              >  4:302484/4‑149 (MQ=255)
                                                                                                                                  aCCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGgcg            >  1:53502/1‑149 (MQ=255)
                                                                                                                                    cAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCATCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCgg                      <  7:42350/137‑1 (MQ=255)
                                                                                                                                    cAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCATCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCgg                      >  8:42350/1‑137 (MQ=255)
                                                                                                                                     aGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAg         >  4:124953/1‑149 (MQ=255)
                                                                                                                                      gcgcACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGa                                      >  5:168755/1‑119 (MQ=255)
                                                                                                                                      gcgcACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGa                                      <  6:168755/119‑1 (MQ=255)
                                                                                                                                          aCCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGCACg    >  5:278941/1‑149 (MQ=255)
                                                                                                                                          aCCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGCACg    <  3:141918/149‑1 (MQ=255)
                                                                                                                                            cGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGCACGCt  >  1:113049/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGCACGCT  >  NC_002947/3614584‑3614872

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: