Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 288,011 T→C 5.3% T33A (ACC→GCC)  tauA ← taurine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947288,0110TC 5.3% 223.4 / 2.5 76T33A (ACC→GCC) tauAtaurine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base C (2/2);  total (39/37)
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.91e-01

CCCAAGTAGCCGATCTGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCTAAATGA  >  NC_002947/287866‑288159
                                                                                                                                                    |                                                                                                                                                    
cccAAGTAGCCGATCTGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTg                                                                                                                                                   <  3:359690/149‑1 (MQ=255)
cccAAGTAGCCGATCTGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTg                                                                                                                                                   >  4:235348/1‑149 (MQ=255)
  cgcgcttccGATCTGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGt                                                                                                                                                 <  6:150980/142‑1 (MQ=255)
        gCCGATCTGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCa                                                                                                                                           >  2:18712/1‑149 (MQ=255)
        gCCGATCTGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCa                                                                                                                                           >  2:140521/1‑149 (MQ=255)
          cGATCTGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACt                                                                                                                                         >  3:122174/1‑149 (MQ=255)
               tGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAg                                                                                                                                    >  3:236126/1‑149 (MQ=255)
                gCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAgg                                                                                                                                   >  7:37154/1‑149 (MQ=255)
                 gaCGGCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCTGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGt                                                                                                                                  <  8:306955/148‑1 (MQ=255)
                 cACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTc                                                                                                                                      <  6:184268/145‑1 (MQ=255)
                 cACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGCCTACCAGACCACTg                                                                                                                                        <  2:186404/143‑1 (MQ=255)
                 cACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGCCTACCAGACCACTGTc                                                                                                                                      <  7:277568/145‑1 (MQ=255)
                    gTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTTgg                                                                                                                               <  3:235348/149‑1 (MQ=255)
                    gTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCgg                                                                                                                               <  1:18712/149‑1 (MQ=255)
                      cgccgcTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGc                                                                                                                              >  1:282990/1‑149 (MQ=255)
                       gccgcTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCCTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCt                                                                                                                             <  2:244692/148‑1 (MQ=255)
                         cgcTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGcc                                                                                                                          >  4:361489/1‑149 (MQ=255)
                           cTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTg                                                                                                                        <  4:114323/149‑1 (MQ=255)
                                cACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGcc                                                                                                                    >  5:338320/1‑148 (MQ=255)
                                cACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGcc                                                                                                                    >  2:360496/1‑148 (MQ=255)
                                     cAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAg                                                                                                              >  2:356238/1‑149 (MQ=255)
                                      attgttCACGGTTGCGCCGTTATAAAATTTGCGCCCTTGCATTCGGGCCCTCCCTTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCa                                                                                                               <  1:304964/144‑1 (MQ=255)
                                         gATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAAcc                                                                                                                <  4:158372/143‑1 (MQ=255)
                                         gATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAAcc                                                                                                                >  3:158372/1‑143 (MQ=255)
                                          aTCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGt                                                                                                         <  8:37154/149‑1 (MQ=255)
                                          aTCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGt                                                                                                         <  2:282990/149‑1 (MQ=255)
                                           tCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTg                                                                                                        <  1:360496/149‑1 (MQ=255)
                                            cACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTc                                                                                                                                      >  3:40773/1‑118 (MQ=255)
                                            cACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTc                                                                                                                                      <  4:40773/118‑1 (MQ=255)
                                            cACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGc                                                                                                       >  2:97206/1‑149 (MQ=255)
                                                   gcgcCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAg                                                                                                >  6:191371/1‑149 (MQ=255)
                                                    cgcCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGa                                                                                                 <  5:209432/147‑1 (MQ=255)
                                                    cgcCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAg                                                                                                <  4:236126/148‑1 (MQ=255)
                                                     gcCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAggg                                                                                              >  8:29227/1‑149 (MQ=255)
                                                         ttATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAg                                                                                          >  3:26872/1‑149 (MQ=255)
                                                           aTCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAggg                                                                                        >  3:316271/1‑149 (MQ=255)
                                                           aTCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAgg                                                                                         <  6:338320/148‑1 (MQ=255)
                                                            tCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAgg                                                                                         >  2:79339/1‑147 (MQ=255)
                                                            tCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAgg                                                                                         <  1:79339/147‑1 (MQ=255)
                                                                aTTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGggcggc                                                                                    >  2:156958/1‑148 (MQ=255)
                                                                     cgcCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAg                                                                              <  4:109261/149‑1 (MQ=255)
                                                                      gcCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGCCAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGGCAGCGCGGCGCAGAgg                                                                             >  3:281382/1‑149 (MQ=255)
                                                                           gTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAggg                                                                                        <  1:368949/134‑1 (MQ=255)
                                                                           gTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAggg                                                                                        >  2:368949/1‑134 (MQ=255)
                                                                              gATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCg                                                                      >  1:100239/1‑149 (MQ=255)
                                                                               aTGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCg                                                                      >  3:318096/1‑148 (MQ=255)
                                                                                          gC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGa                                                          >  2:70244/1‑149 (MQ=255)
                                                                                              tGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTgg                                                                                                                       <  2:275789/86‑1 (MQ=255)
                                                                                              tGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTgg                                                                                                                       >  1:275789/1‑86 (MQ=255)
                                                                                                gCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGTCAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCAt                                                      >  1:261209/1‑149 (MQ=255)
                                                                                                gCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGGCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGa                                                                    <  3:25193/135‑1 (MQ=255)
                                                                                                gCTTTTTCGTAGTCGCCATAGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCAACTGTCAGGTCGGCTGCATGGGCCAACCAGCTGGTGCTGGAGAGGGTCAGGTCGGAGAACAGACGCAGCGGa                                                                    >  4:25193/1‑135 (MQ=255)
                                                                                                 cTTTTTCGTAGTCGCCATCGACCTGTGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATg                                                     <  4:126998/149‑1 (MQ=255)
                                                                                                  tttttAGTAGTCGCCATCTACCTGCGCCACTTTGGCCGTGTCGACAGTGGTCTGGTAGTACAATGTAAGGTCGTCTGCCTGGGACAACAAGCTGGTGCTGGCGAGGGGAATTGCGGCGAACAGGCGCAGCGTAGCGTGCGGGATCATgg                                                    >  3:50709/1‑149 (MQ=255)
                                                                                                    tttCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATgg                                                    <  3:361489/147‑1 (MQ=255)
                                                                                                      tCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTg                                                 <  1:140521/148‑1 (MQ=255)
                                                                                                       cGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGaa                                               >  7:167724/1‑149 (MQ=255)
                                                                                                        gTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAAc                                              <  4:26872/149‑1 (MQ=255)
                                                                                                         tAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGCGCGTGCGGGCTCATGGTTgcagg                                             >  2:266575/1‑145 (MQ=255)
                                                                                                            tCGCCATGGACCTGCGCGAGTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCTCAGCGGAGCGTGCGGGATCATGGTTGaa                                               <  7:305250/144‑1 (MQ=255)
                                                                                                            tCGCCATCGACCTGCGCCACTTTGGCCCGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGCCGCAGCGGAGCGTGCGGGATCAAGCTAGAAcaca                                           >  8:305250/1‑145 (MQ=255)
                                                                                                                cATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGt                                       >  7:31065/1‑148 (MQ=255)
                                                                                                                cATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTc                                      <  5:191371/149‑1 (MQ=255)
                                                                                                                cATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTc                                      <  1:70244/149‑1 (MQ=255)
                                                                                                                        tGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGc                                  <  8:167013/145‑1 (MQ=255)
                                                                                                                        tGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGc                                  >  7:167013/1‑145 (MQ=255)
                                                                                                                            ccACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGAGCTGGCGAGGGTCATGGCGGCCAACAGGCGCAGCGGAGCGTGCGGGATCATGGTAGAACCTCACGTCATGCAAGTCgtg                          >  5:240103/1‑149 (MQ=255)
                                                                                                                             cACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCgtgt                         >  3:227338/1‑149 (MQ=255)
                                                                                                                              aCTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGtt                        <  2:100239/149‑1 (MQ=255)
                                                                                                                                    gCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACGGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATaa                  <  4:122174/149‑1 (MQ=255)
                                                                                                                                    gCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATaa                  <  4:227338/149‑1 (MQ=255)
                                                                                                                                          ctgACAGTGGCCTACCAGACCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACgag              >  1:346963/3‑147 (MQ=255)
                                                                                                                                          ctgACAGTGGCCTACCAGACCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACgag              >  5:291341/3‑147 (MQ=255)
                                                                                                                                                 tGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCTaaa     >  6:89977/1‑149 (MQ=255)
                                                                                                                                                 tGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCTaaa     <  8:31065/149‑1 (MQ=255)
                                                                                                                                                   gTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCTAAATg   <  4:318096/149‑1 (MQ=255)
                                                                                                                                                    tCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCTAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCTAAATGa  <  1:356238/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCAAGTAGCCGATCTGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGC‑‑TGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCTAAATGA  >  NC_002947/287866‑288159

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: