Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 381,994 A→T 13.9% Q111L (CAG→CTG)  tauB → taurine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273381,9940AT13.9% 546.4 / 69.6 194Q111L (CAG→CTG) tauBtaurine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (75/92);  new base T (12/15);  total (87/107)
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.76e-01

Only 100 of 195 total aligned reads displayed.

GCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGAACCGTTTGGTGCGCTGGA  >  CP009273/381851‑382141
                                                                                                                                               |                                                                                                                                                   
gCGTTTTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTTTTGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATg                                                                                                                                                 <  2:1281668/148‑1 (MQ=255)
gCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGATGGGCTACTACCGTGGCGCAATGTACAGCACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATg                                                                                                                                                 <  1:2649478/148‑1 (MQ=255)
gCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCTGATg                                                                                                                                                 <  2:563128/148‑1 (MQ=255)
gCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATg                                                                                                                                                 <  1:3179814/148‑1 (MQ=255)
gCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATg                                                                                                                                                 >  2:2725662/1‑148 (MQ=255)
 cGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGc                                                                                                                                                >  2:2905160/1‑148 (MQ=255)
    aTTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGa                                                                                                                                             >  2:335039/1‑148 (MQ=255)
    aTTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGa                                                                                                                                             >  1:3392505/1‑148 (MQ=255)
         gggACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGaaaaaa                                                                                                                                        <  2:1527637/148‑1 (MQ=255)
           gACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCATTGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGt                                                                                                                                      <  2:3402715/148‑1 (MQ=255)
            aCCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCTGATGCTGAAAAAAGTg                                                                                                                                     >  1:2391946/1‑148 (MQ=255)
            aCCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTg                                                                                                                                     <  2:3575864/148‑1 (MQ=255)
            aCCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTg                                                                                                                                     <  2:3184161/148‑1 (MQ=255)
            aCCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTg                                                                                                                                     <  2:1601375/148‑1 (MQ=255)
             ccGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGTTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTgg                                                                                                                                    <  1:89843/148‑1 (MQ=255)
             ccGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTgg                                                                                                                                    <  2:1335810/148‑1 (MQ=255)
             ccGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTgg                                                                                                                                    <  1:1256617/148‑1 (MQ=255)
             ccGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATGGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTgg                                                                                                                                    <  2:3559899/148‑1 (MQ=255)
               gggAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTgggg                                                                                                                                  <  2:492637/148‑1 (MQ=255)
                ggAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCTGATGCTGAAAAAAGTGGGGc                                                                                                                                 >  1:20018/1‑148 (MQ=255)
                 gAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCt                                                                                                                                >  2:3269137/1‑148 (MQ=255)
                          gTGGCGTATTTTTTCAGAATGAAGGGCTACTACCGTGGCGAAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGTTATATAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGgcgc                                                                                                                       <  1:2589238/148‑1 (MQ=255)
                          gTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGgcgc                                                                                                                       <  1:2220709/148‑1 (MQ=255)
                                tAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGaaaa                                                                                                                 <  1:1514332/148‑1 (MQ=255)
                                 aGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCATTTGGCAGTTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGaaaaa                                                                                                                <  1:746826/148‑1 (MQ=255)
                                   tttttCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACg                                                                                                              <  2:2826640/148‑1 (MQ=255)
                                      ttCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTa                                                                                                           <  2:1660676/148‑1 (MQ=255)
                                        cAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCTGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACa                                                                                                         >  2:3694130/1‑148 (MQ=255)
                                         aGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACAt                                                                                                        >  1:3020895/1‑148 (MQ=255)
                                           aaTGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCt                                                                                                      >  2:3141332/1‑148 (MQ=255)
                                           aaTGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTAGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCt                                                                                                      >  2:1237127/1‑148 (MQ=255)
                                              gAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGc                                                                                                   <  1:1882850/148‑1 (MQ=255)
                                                 gggCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGc                                                                                                <  2:2997967/148‑1 (MQ=255)
                                                  ggCTACTCCGGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCTGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCt                                                                                               <  1:996027/148‑1 (MQ=255)
                                                  ggCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCt                                                                                               <  1:3151466/148‑1 (MQ=255)
                                                  ggCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCt                                                                                               >  1:3688181/1‑148 (MQ=255)
                                                  ggCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCt                                                                                               >  2:421823/1‑148 (MQ=255)
                                                   gCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCtt                                                                                              <  2:2883115/148‑1 (MQ=255)
                                                    ctactaCCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCTGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCttt                                                                                             <  2:20018/148‑1 (MQ=255)
                                                    ctactaCCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCttt                                                                                             >  2:2605424/1‑148 (MQ=255)
                                                     tactacCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGTCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTc                                                                                            <  2:2105485/148‑1 (MQ=255)
                                                        tacCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCgg                                                                                         <  1:505199/148‑1 (MQ=255)
                                                         acCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCTGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCggt                                                                                        >  2:2909331/1‑148 (MQ=255)
                                                          ccGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCggtg                                                                                       >  2:858374/1‑148 (MQ=255)
                                                           cGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCggtgg                                                                                      >  1:673189/1‑148 (MQ=255)
                                                             tGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTc                                                                                    >  2:147031/1‑148 (MQ=255)
                                                               gcgcAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCaa                                                                                  <  2:3475082/148‑1 (MQ=255)
                                                                 gcAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACg                                                                                <  1:2741099/148‑1 (MQ=255)
                                                                  cAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGt                                                                               >  2:2349513/1‑148 (MQ=255)
                                                                  cAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGt                                                                               <  2:2029220/148‑1 (MQ=255)
                                                                    aTGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGTTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCa                                                                             <  1:2306627/148‑1 (MQ=255)
                                                                    aTGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCa                                                                             >  2:685560/1‑148 (MQ=255)
                                                                    aTGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCa                                                                             >  1:1078119/1‑148 (MQ=255)
                                                                     tGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAg                                                                            >  2:885065/1‑148 (MQ=255)
                                                                        aCAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCgg                                                                         >  2:678359/1‑148 (MQ=255)
                                                                            gACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTgg                                                                     >  2:3153796/1‑148 (MQ=255)
                                                                            gACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTgg                                                                     >  2:3153710/1‑148 (MQ=255)
                                                                                aCGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGAt                                                                 >  2:2654166/1‑148 (MQ=255)
                                                                                 cGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCATCGGGTGGGGAtt                                                                >  1:1724152/1‑148 (MQ=255)
                                                                                 cGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGAtt                                                                >  1:1269792/1‑148 (MQ=255)
                                                                                  gTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTg                                                               <  1:1923455/148‑1 (MQ=255)
                                                                                      cGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCTGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCg                                                           <  1:863623/148‑1 (MQ=255)
                                                                                      cGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCg                                                           >  1:2107940/1‑148 (MQ=255)
                                                                                           ggCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTgcgc                                                      >  1:1456611/1‑148 (MQ=255)
                                                                                             ccTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTg                                                    <  1:3053425/148‑1 (MQ=255)
                                                                                               tGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCTGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTggc                                                  >  1:3621798/1‑148 (MQ=255)
                                                                                               tGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTggc                                                  >  1:1635944/1‑148 (MQ=255)
                                                                                                gCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTggcg                                                 >  2:1514702/1‑148 (MQ=255)
                                                                                                    ttGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCTGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGa                                             <  1:1081155/148‑1 (MQ=255)
                                                                                                    ttGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGa                                             >  2:1721172/1‑148 (MQ=255)
                                                                                                      ggCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAAt                                           <  1:3469714/148‑1 (MQ=255)
                                                                                                        cAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCTGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATcc                                         >  2:2090284/1‑148 (MQ=255)
                                                                                                         aGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATccc                                        <  1:2862644/148‑1 (MQ=255)
                                                                                                            tataGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAg                                     >  2:2468044/1‑148 (MQ=255)
                                                                                                             ataGAGAAAATGCAGCGACTGGAAATCGCGCACCTGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGc                                    <  2:3390316/148‑1 (MQ=255)
                                                                                                              taGAGAAAATGCAGCGACTGGAAATCGCGCACCTGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCt                                   <  1:2541118/148‑1 (MQ=255)
                                                                                                                    aaaTGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGttatt                             <  2:830287/148‑1 (MQ=255)
                                                                                                                    aaaTGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGttatt                             <  2:1794858/148‑1 (MQ=255)
                                                                                                                     aaTGCAGCGACTGGAAATCGCGCACCTGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCAAGCTGttatta                            <  1:2090284/148‑1 (MQ=255)
                                                                                                                       tGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACt                          >  1:1572154/1‑148 (MQ=255)
                                                                                                                           gCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTcgac                      >  1:2408195/1‑148 (MQ=255)
                                                                                                                           gCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTcgac                      >  1:487134/1‑148 (MQ=255)
                                                                                                                            cGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTcgacg                     >  2:2264204/1‑148 (MQ=255)
                                                                                                                             gACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTcgacga                    <  1:357030/148‑1 (MQ=255)
                                                                                                                             gACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTcgacga                    <  2:1588253/148‑1 (MQ=255)
                                                                                                                              aCTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGaa                   <  2:1217848/148‑1 (MQ=255)
                                                                                                                                  gAAATCGCGCACCTGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGAACCGt               >  1:1008912/1‑148 (MQ=255)
                                                                                                                                  gAAATCGCGCACCTGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGAACCGt               >  1:120613/1‑148 (MQ=255)
                                                                                                                                  gAAATCGCGCACCTGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGAACCGt               >  1:1009875/1‑148 (MQ=255)
                                                                                                                                  gAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGAACCGt               <  1:2670552/148‑1 (MQ=255)
                                                                                                                                   aaaTCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGAACCGtt              <  1:668784/148‑1 (MQ=255)
                                                                                                                                   aaaTCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGAACCGtt              <  1:3269191/148‑1 (MQ=255)
                                                                                                                                    aaTCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGAACCGttt             <  1:795077/148‑1 (MQ=255)
                                                                                                                                     aTCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGAACCGTTTg            <  2:762125/148‑1 (MQ=255)
                                                                                                                                      tCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGAACCGTTTgg           <  1:3651705/148‑1 (MQ=255)
                                                                                                                                         cgcACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGAACCGTTTGGTgc        >  1:2253400/1‑148 (MQ=255)
                                                                                                                                          gcACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGAACCGTTTGGTgcg       <  2:1199333/148‑1 (MQ=255)
                                                                                                                                             ccAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGAACCGTTTGGTGCGCTg    >  1:575181/1‑148 (MQ=255)
                                                                                                                                              cAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGAACCGTTTGGTGCGCTgg   >  2:2440979/1‑148 (MQ=255)
                                                                                                                                               aGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGAACCGTTTGGTGCGCTGGa  >  2:2336068/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GCGTATTGAGGGACCGGGAGCAGAGCGTGGCGTAGTTTTTCAGAATGAAGGGCTACTACCGTGGCGCAATGTACAGGACAACGTGGCGTTCGGCCTGCAATTGGCAGGTATAGAGAAAATGCAGCGACTGGAAATCGCGCACCAGATGCTGAAAAAAGTGGGGCTGGAAGGCGCAGAAAAACGCTACATCTGGCAGCTTTCCGGTGGTCAACGTCAGCGGGTGGGGATTGCTCGTGCGCTGGCGGCGAATCCCCAGCTGTTATTACTCGACGAACCGTTTGGTGCGCTGGA  >  CP009273/381851‑382141

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