Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 288,015 G→C 5.1% Y31* (TAC→TAG) ‡ 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,0150GC 5.1% 254.1 / 3.5 79Y31* (TAC→TAG) ‡tauAtaurine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (35/39);  new base C (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

CAAGTAGCCGATCTGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCTAAAT  >  NC_002947/287868‑288157
                                                                                                                                                   |                                                                                                                                              
cAAGTAGCCGATCTGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTgg                                                                                                                                                >  1:189434/1‑148 (MQ=255)
 aaGTAGCCGATCTGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTa                                                                                                                                              <  5:3106/149‑1 (MQ=255)
   gTAGCCGATCTGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAgg                                                                                                                                            <  2:189434/149‑1 (MQ=255)
         gATCTGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAgg                                                                                                                                            >  1:167086/1‑143 (MQ=255)
         gATCTGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAgg                                                                                                                                            <  2:167086/143‑1 (MQ=255)
            cTGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTc                                                                                                                                    <  2:283435/148‑1 (MQ=255)
              gCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAgg                                                                                                                                 >  1:389639/1‑149 (MQ=255)
                aCGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTc                                                                                                                               <  7:284841/149‑1 (MQ=255)
                 cGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGc                                                                                                                                           <  6:176625/136‑1 (MQ=255)
                 cGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGc                                                                                                                                           >  5:176625/1‑136 (MQ=255)
                   tCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGc                                                                                                                            >  4:253332/1‑149 (MQ=255)
                   tCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGc                                                                                                                            >  2:1891/1‑149 (MQ=255)
                            gCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTgg                                                                                                                     <  2:292643/147‑1 (MQ=255)
                            gCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGc                                                                                                                   <  3:395606/149‑1 (MQ=255)
                            gCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGc                                                                                                                   >  4:281071/1‑149 (MQ=255)
                             ccACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAg                                                                                                                                  <  5:86741/133‑1 (MQ=255)
                             ccACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAg                                                                                                                                  >  6:86741/1‑133 (MQ=255)
                             ccACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTg                                                                                                                      >  6:382386/1‑145 (MQ=255)
                             ccACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTg                                                                                                                      <  5:382386/145‑1 (MQ=255)
                              cACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGcc                                                                                                                  >  6:122376/1‑148 (MQ=255)
                               aCGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCaa                                                                                                                >  3:358251/1‑149 (MQ=255)
                                cGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTTGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCt                                                                                                                           <  6:53992/137‑1 (MQ=255)
                                cGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCt                                                                                                                           >  5:53992/1‑137 (MQ=255)
                                cGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCt                                                                                                                           >  5:75338/1‑137 (MQ=255)
                                cGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCt                                                                                                                           <  6:75338/137‑1 (MQ=255)
                                cGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAAc                                                                                                               <  2:65243/149‑1 (MQ=255)
                                        aTCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGt                                                                                                       <  2:68725/149‑1 (MQ=255)
                                         tCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGt                                                                                                                                               >  4:386698/1‑108 (MQ=255)
                                         tCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGt                                                                                                                                               <  3:386698/108‑1 (MQ=255)
                                            cGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTg                                                                                                   <  4:358251/149‑1 (MQ=255)
                                                          tCATATTTTCGGCAGTGGATGCTGGAGTTGCTGGCTTTTTAGTTGTCGGCATCGTCCTGCGACTCTTTGGCCGGGTCGACAGTGGTCTGTTAGGGGACTGTCAGGGCGGCTGCCTGTTCCAAGCATCTTGTGCTTGCGAGGGTAAgtt                                                                                      <  8:337435/148‑3 (MQ=255)
                                                          tCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTTGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCACCCAGCTGGTGCTGGCGAGGGTCAGggc                                                                                     >  8:294175/1‑149 (MQ=255)
                                                          tCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAggg                                                                                      <  7:95236/148‑1 (MQ=255)
                                                            aaaTTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGggcgg                                                                                   <  5:111221/149‑1 (MQ=255)
                                                            aaaTTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCAGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGggcgg                                                                                   >  5:51689/1‑149 (MQ=255)
                                                              aTTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGgcggcg                                                                                 <  1:281705/149‑1 (MQ=255)
                                                                   cgcCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAg                                                                            <  6:3285/149‑1 (MQ=255)
                                                                    gcCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAgg                                                                           >  6:360933/1‑149 (MQ=255)
                                                                      cAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGgcg                                                                         <  3:253332/149‑1 (MQ=255)
                                                                      cAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGgcg                                                                         <  1:1891/149‑1 (MQ=255)
                                                                       aGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGgcg                                                                         >  7:290302/1‑148 (MQ=255)
                                                                             tGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGa                                                                  <  2:389639/149‑1 (MQ=255)
                                                                             tGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGa                                                                  <  5:122376/149‑1 (MQ=255)
                                                                                 ggCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGTGTCAGGGCGGCGAACAGGCGCAGCGGa                                                                  <  2:262218/145‑1 (MQ=255)
                                                                                 ggCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGa                                                                  >  1:262218/1‑145 (MQ=255)
                                                                                     ttGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCgg                                                          >  2:327796/1‑149 (MQ=255)
                                                                                          gacTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGgcgc                                                                        >  7:113442/3‑130 (MQ=255)
                                                                                          gacTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGgcgc                                                                        <  8:113442/128‑1 (MQ=255)
                                                                                           gCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCa                                                     <  1:215761/148‑1 (MQ=255)
                                                                                            cTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATg                                                   <  8:290302/149‑1 (MQ=255)
                                                                                               tttCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGt                                                 <  5:360933/148‑1 (MQ=255)
                                                                                                       tCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCg                                                           <  3:102994/130‑1 (MQ=255)
                                                                                                       tCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCg                                                           >  4:102994/1‑130 (MQ=255)
                                                                                                         gCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCtctc                                       <  3:205572/148‑1 (MQ=255)
                                                                                                           cATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGt                                     >  5:39284/1‑148 (MQ=255)
                                                                                                           cATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCGGGTGCTGGCGAGGGTCAGGGCGGCGCCCAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTc                                    >  7:162529/1‑149 (MQ=255)
                                                                                                            aTCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTTCGGGATCATGGTTGAACCTCTCGt                                     >  7:107048/1‑147 (MQ=255)
                                                                                                            aTCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTGTCGt                                     >  2:344936/1‑147 (MQ=255)
                                                                                                            aTCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGt                                     >  5:49562/1‑147 (MQ=255)
                                                                                                               gACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAgg                                 <  3:281071/148‑1 (MQ=255)
                                                                                                                aCCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGc                                >  2:275581/1‑148 (MQ=255)
                                                                                                                  cTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAg                             >  2:34789/1‑149 (MQ=255)
                                                                                                                  cTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAg                             <  1:2445/149‑1 (MQ=255)
                                                                                                                   tGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGt                            >  1:26104/1‑149 (MQ=255)
                                                                                                                   tGCGCCACTTTGGCCGGGTCGACAGTGGTATGGTAGGCCACTGTCAGGTCGGCTGCATGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCACGCAAGt                            >  2:133629/1‑149 (MQ=255)
                                                                                                                         aCTTTGGCCGGGTCGACAGTGGTCTGGTAGGCGACTGTCAGGTCGTCTGCCTGGTCCAACGAGCTGGTGCTGGGGAGTGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCTTGATCATGGTTGAACCTGTCGTCAGGCAAGTCGTGtt                      <  1:133629/149‑1 (MQ=255)
                                                                                                                           tttGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTg                     >  5:166314/1‑148 (MQ=255)
                                                                                                                                ccGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGTTCATGGTTGAACctc                                         >  5:354772/1‑123 (MQ=255)
                                                                                                                                ccGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACctc                                         <  6:354772/123‑1 (MQ=255)
                                                                                                                                ccGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAAc               <  1:327796/149‑1 (MQ=255)
                                                                                                                                    cctgACAGTGGCCTACCAGACCACTGTCGGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCa                                   <  6:213938/122‑1 (MQ=255)
                                                                                                                                    cctgACAGTGGCCTACCAGACCACTGTCGGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCa                                   >  5:213938/4‑125 (MQ=255)
                                                                                                                                       gACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAgg        <  6:166314/149‑1 (MQ=255)
                                                                                                                                       gACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAgg        >  3:306800/1‑149 (MQ=255)
                                                                                                                                       gACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCCTCAGGCAAGTCGTGTTGGATAACGAGAAgg        <  1:34789/149‑1 (MQ=255)
                                                                                                                                        aCAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGc       <  6:39284/149‑1 (MQ=255)
                                                                                                                                        aCAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGc       >  5:361629/1‑149 (MQ=255)
                                                                                                                                        aCAGTGGCCTACCAGACCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACga             >  6:210070/1‑143 (MQ=255)
                                                                                                                                        aCAGTGGCCTACCAGACCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACga             <  5:210070/143‑1 (MQ=255)
                                                                                                                                              gTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCTAAAt  >  5:19856/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CAAGTAGCCGATCTGCACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCTAAAT  >  NC_002947/287868‑288157

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: