Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 381,815 G→T 9.6% V51V (GTG→GTT)  tauB → taurine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273381,8150GT 9.6% 905.2 / 48.0 255V51V (GTG→GTT) tauBtaurine transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (122/108);  new base T (12/13);  total (134/121)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 268 total aligned reads displayed.

CTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAG  >  CP009273/381666‑381963
                                                                                                                                                     |                                                                                                                                                    
ctgcAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTtgtg                                                                                                                                                      <  1:3649826/150‑1 (MQ=255)
ctgcAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTtgtg                                                                                                                                                      <  1:1169881/150‑1 (MQ=255)
ctgcAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTtgtg                                                                                                                                                      >  2:2583043/1‑150 (MQ=255)
 tgcAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGc                                                                                                                                                     <  1:4493888/150‑1 (MQ=255)
  gcAAAGCTCTCATCTTTACGCCGAATATGTCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGcc                                                                                                                                                    <  2:2632687/150‑1 (MQ=255)
      aTCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTAt                                                                                                                                                <  2:3383400/150‑1 (MQ=255)
       tctctcATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATc                                                                                                                                               <  1:1423556/150‑1 (MQ=255)
       tctctcATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATc                                                                                                                                               >  1:707146/1‑150 (MQ=255)
           tcATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCa                                                                                                                                           >  2:471200/1‑150 (MQ=255)
           tcATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCa                                                                                                                                           >  2:471769/1‑150 (MQ=255)
           tcATCTTTACGCCGAGTATGGCGGCAAACCGGCACTGGAAGATATCATCCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCa                                                                                                                                           >  2:4400587/1‑150 (MQ=255)
            cATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCAt                                                                                                                                          >  1:2197263/1‑150 (MQ=255)
            cATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCAt                                                                                                                                          <  2:2740282/150‑1 (MQ=255)
            cATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCAt                                                                                                                                          <  2:2740476/150‑1 (MQ=255)
            cATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCAt                                                                                                                                          <  2:2740644/150‑1 (MQ=255)
               cTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGgc                                                                                                                                       >  1:3608951/1‑150 (MQ=255)
                 ttACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGgcag                                                                                                                                     <  1:2254039/150‑1 (MQ=255)
                  tACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGgcagc                                                                                                                                    >  1:4177582/1‑150 (MQ=255)
                     gCCGATTATGGCGGCAAACCGGCACTGGAATATAACAAGCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCAtt                                                                                                                                 <  2:2910720/150‑1 (MQ=255)
                     gCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCAtt                                                                                                                                 <  2:3987967/150‑1 (MQ=255)
                      ccGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTc                                                                                                                                >  2:289910/1‑150 (MQ=255)
                       cGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCa                                                                                                                               >  2:3171008/1‑150 (MQ=255)
                       cGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCa                                                                                                                               <  1:4533828/150‑1 (MQ=255)
                       cGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCa                                                                                                                               >  2:4347891/1‑150 (MQ=255)
                        gATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTTCCTTATCAGCATGGCAGCATTCaa                                                                                                                              >  2:4054554/1‑150 (MQ=255)
                          ttATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTTCCTTATCAGCATGGCAGCATTCAACt                                                                                                                            <  1:933135/150‑1 (MQ=255)
                          ttATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTTCCTTATCAGCATGGCAGCATTCAACt                                                                                                                            <  1:932986/150‑1 (MQ=255)
                           tATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTg                                                                                                                           >  2:2323556/1‑150 (MQ=255)
                              ggcggcAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCg                                                                                                                        <  1:3048941/150‑1 (MQ=255)
                                   cAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTaa                                                                                                                   >  1:3998988/1‑150 (MQ=255)
                                    aaaCCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAg                                                                                                                  >  2:917913/1‑150 (MQ=255)
                                    aaaCCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAg                                                                                                                  >  2:4252272/1‑150 (MQ=255)
                                      aCCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCg                                                                                                                <  1:139030/150‑1 (MQ=255)
                                         ggCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTTCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTAt                                                                                                             <  2:3739156/150‑1 (MQ=255)
                                         ggCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTAt                                                                                                             <  2:4009225/150‑1 (MQ=255)
                                         ggCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCATATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTAt                                                                                                             >  1:1064355/1‑150 (MQ=255)
                                          gCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTAtt                                                                                                            >  2:1127840/1‑150 (MQ=255)
                                            aCTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGa                                                                                                          >  2:1058603/1‑150 (MQ=255)
                                              tGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAgg                                                                                                        >  1:2773082/1‑150 (MQ=255)
                                                gAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGa                                                                                                      >  1:487286/1‑150 (MQ=255)
                                                       tCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAg                                                                                               <  1:2161197/150‑1 (MQ=255)
                                                        cAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGc                                                                                              >  1:4458258/1‑150 (MQ=255)
                                                        cAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGc                                                                                              <  2:1997167/150‑1 (MQ=255)
                                                         aaCCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCa                                                                                             >  1:306684/1‑150 (MQ=255)
                                                          aCCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCag                                                                                            <  1:1000979/150‑1 (MQ=255)
                                                           ccTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCaga                                                                                           <  2:1074694/150‑1 (MQ=255)
                                                           ccTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCaga                                                                                           <  1:4211003/150‑1 (MQ=255)
                                                            cTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCagag                                                                                          >  1:1317772/1‑150 (MQ=255)
                                                            cTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCagag                                                                                          >  1:2087251/1‑150 (MQ=255)
                                                              gACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCg                                                                                        <  2:2333865/150‑1 (MQ=255)
                                                                  cTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTTCCTTATCAGCATGGCAGCATTCTACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGc                                                                                    <  1:2933727/150‑1 (MQ=255)
                                                                  cTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGc                                                                                    <  2:3867064/150‑1 (MQ=255)
                                                                    ggAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTTCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGt                                                                                  >  1:3871869/1‑150 (MQ=255)
                                                                     gAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTa                                                                                 >  1:2063717/1‑150 (MQ=255)
                                                                       aaGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTTCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGt                                                                               <  2:2389316/150‑1 (MQ=255)
                                                                        aGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGtt                                                                              >  1:3659010/1‑150 (MQ=255)
                                                                          cggcgAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGtttt                                                                            <  2:2439033/150‑1 (MQ=255)
                                                                           ggcgAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGttttt                                                                           >  2:565832/1‑150 (MQ=255)
                                                                              gAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAg                                                                        <  2:1025351/150‑1 (MQ=255)
                                                                              gAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAg                                                                        <  2:1602909/150‑1 (MQ=255)
                                                                               aGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGa                                                                       <  1:587142/150‑1 (MQ=255)
                                                                                gCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGaa                                                                      <  1:2279114/150‑1 (MQ=255)
                                                                                  tACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATg                                                                    <  2:4453305/150‑1 (MQ=255)
                                                                                   aCTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGa                                                                   >  1:3895055/1‑150 (MQ=255)
                                                                                   aCTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGa                                                                   >  2:186470/1‑150 (MQ=255)
                                                                                   aCTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGa                                                                   >  1:3894651/1‑150 (MQ=255)
                                                                                         ggtgCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGct                                                             <  2:3946815/150‑1 (MQ=255)
                                                                                             cTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGctacta                                                         <  1:1706221/150‑1 (MQ=255)
                                                                                             cTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGctacta                                                         >  1:329651/1‑150 (MQ=255)
                                                                                               ggggCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTAcc                                                       <  1:1751577/150‑1 (MQ=255)
                                                                                                  gCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTg                                                    <  1:3722805/150‑1 (MQ=255)
                                                                                                   ccGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTgg                                                   >  2:4153605/1‑150 (MQ=255)
                                                                                                    cGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGgc                                                  <  2:1914704/150‑1 (MQ=255)
                                                                                                     gTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGgcg                                                 >  2:1712074/1‑150 (MQ=255)
                                                                                                      tCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGgcgc                                                >  1:4118889/1‑150 (MQ=255)
                                                                                                       ccGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCa                                               >  2:912976/1‑150 (MQ=255)
                                                                                                         ggCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAAt                                             >  1:320071/1‑150 (MQ=255)
                                                                                                          gCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTTCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATg                                            >  2:616666/1‑150 (MQ=255)
                                                                                                              cGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACa                                        >  2:3408454/1‑150 (MQ=255)
                                                                                                                    aaCCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACaa                                  >  2:2264281/1‑150 (MQ=255)
                                                                                                                         cccTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTgg                             >  2:418077/1‑150 (MQ=255)
                                                                                                                         cccTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTgg                             <  2:1317772/150‑1 (MQ=255)
                                                                                                                           ctgctgAATCTGATTGCCGGTTTTGTTCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCg                           <  1:3029725/150‑1 (MQ=255)
                                                                                                                           ctgctgAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCg                           <  2:4150408/150‑1 (MQ=255)
                                                                                                                             gctgAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGtt                         <  2:3822474/150‑1 (MQ=255)
                                                                                                                              ctgAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTc                        <  2:4020738/150‑1 (MQ=255)
                                                                                                                              ctgAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTc                        >  2:2108954/1‑150 (MQ=255)
                                                                                                                              ctgAATCTGATTGCCGGATTTGTGCCTTATCAGCATGGCTGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTc                        <  2:756489/150‑1 (MQ=255)
                                                                                                                               tgAATCTGATTGCCGGTTTTGTTCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCg                       >  1:1261957/1‑150 (MQ=255)
                                                                                                                                gAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACATTACAACGTGGCGTTctg                      >  2:168681/1‑148 (MQ=255)
                                                                                                                                gAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCgg                      >  1:989429/1‑150 (MQ=255)
                                                                                                                                      gATTGCCGGTTTTGTTCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCa                >  2:2979805/1‑150 (MQ=255)
                                                                                                                                       aTTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCaa               >  2:2590580/1‑150 (MQ=255)
                                                                                                                                          gCCGGTTTTGTTCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTg            <  1:2170650/150‑1 (MQ=255)
                                                                                                                                          gCCGGTTTTGTTCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTg            <  1:2570852/150‑1 (MQ=255)
                                                                                                                                          gCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTg            <  1:3171007/150‑1 (MQ=255)
                                                                                                                                          gCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTg            <  1:4347891/150‑1 (MQ=255)
                                                                                                                                             ggTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCa         >  2:1298328/1‑150 (MQ=255)
                                                                                                                                             ggTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCa         >  2:3567232/1‑150 (MQ=255)
                                                                                                                                              gTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGTCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAg        >  2:1455010/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGCAAATCTCTCATCTTTACGCCGATTATGGCGGCAAACCGGCACTGGAAGATATCAACCTGACGCTGGAAAGCGGCGAGCTACTGGTGGTGCTGGGGCCGTCCGGCTGCGGTAAAACCACCCTGCTGAATCTGATTGCCGGTTTTGTGCCTTATCAGCATGGCAGCATTCAACTGGCGGGTAAGCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAG  >  CP009273/381666‑381963

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑