Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 325,936 A→T 7.3% K255* (AAG→TAG)  betT → choline transporter of high affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273325,9360AT 7.3% 1264.4 / 65.4 350K255* (AAG→TAG) betTcholine transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base A (171/153);  new base T (17/9);  total (188/162)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 358 total aligned reads displayed.

TATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCTTTATGGGCATGACGCT  >  CP009273/325788‑326084
                                                                                                                                                    |                                                                                                                                                    
 aTTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTNTGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGNTCTCTGTCACCTCCGGTGTCGATaa                                                                                                                                                     >  1:3124659/1‑149 (MQ=255)
 aTTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATaa                                                                                                                                                     <  1:550622/149‑1 (MQ=255)
 aTTGCCACCACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATaa                                                                                                                                                     >  2:4137777/1‑149 (MQ=255)
    gCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAggg                                                                                                                                                  >  1:3975171/1‑149 (MQ=255)
    gCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAggg                                                                                                                                                  >  1:2618830/1‑149 (MQ=255)
       aCTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCAt                                                                                                                                               <  1:3340856/149‑1 (MQ=255)
        cTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCAtt                                                                                                                                              <  2:2302252/149‑1 (MQ=255)
              tCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATTAGGGCATTCGCgtg                                                                                                                                        >  1:4925981/1‑149 (MQ=255)
               cGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCgtgt                                                                                                                                       <  2:598294/149‑1 (MQ=255)
                   aTCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATc                                                                                                                                   >  2:4494773/1‑149 (MQ=255)
                    tCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCg                                                                                                                                  <  1:1304654/149‑1 (MQ=255)
                        tgtgGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGc                                                                                                                              <  1:2065938/149‑1 (MQ=255)
                        tgtgGTGCAGCCTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGc                                                                                                                              >  1:3232151/1‑149 (MQ=255)
                                  cTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTcgc                                                                                                                    >  2:3802189/1‑149 (MQ=255)
                                   ttAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTcgcg                                                                                                                   >  1:3418239/1‑149 (MQ=255)
                                   ttAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTcgcg                                                                                                                   >  2:292843/1‑149 (MQ=255)
                                   ttAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTcgcg                                                                                                                   >  1:1901349/1‑149 (MQ=255)
                                    tAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCgcgc                                                                                                                  <  1:2117383/149‑1 (MQ=255)
                                     aaCTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCgcgct                                                                                                                 >  2:1071157/1‑149 (MQ=255)
                                     aaCTATGGATTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCgcgct                                                                                                                 <  2:2792388/149‑1 (MQ=255)
                                        tATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCgcgctggc                                                                                                              >  1:4862370/1‑149 (MQ=255)
                                           ggCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCgcgctggcgct                                                                                                           >  2:5078278/1‑149 (MQ=255)
                                              ttGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTggg                                                                                                        >  1:703845/1‑149 (MQ=255)
                                                gAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGAt                                                                                                      <  1:3513867/149‑1 (MQ=255)
                                                   cGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGa                                                                                                   <  2:4500148/149‑1 (MQ=255)
                                                   cGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGa                                                                                                   <  2:957374/149‑1 (MQ=255)
                                                      aCTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATcc                                                                                                >  2:3082075/1‑149 (MQ=255)
                                                      aCTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATcc                                                                                                >  2:5188255/1‑149 (MQ=255)
                                                       cTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCt                                                                                               >  1:3679746/1‑149 (MQ=255)
                                                        tGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATTAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTg                                                                                              >  1:1515368/1‑149 (MQ=255)
                                                        tGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGACAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTg                                                                                              >  2:4188830/1‑149 (MQ=255)
                                                         gTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGt                                                                                             >  2:3390024/1‑149 (MQ=255)
                                                         gTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGt                                                                                             <  2:4971470/149‑1 (MQ=255)
                                                           ttGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTc                                                                                           <  2:765941/149‑1 (MQ=255)
                                                           ttGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTc                                                                                           <  2:2175417/149‑1 (MQ=255)
                                                             gATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGt                                                                                         >  2:5043519/1‑149 (MQ=255)
                                                                 ttCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGTGATTGATCCTGTGCGTATTg                                                                                     >  2:4369721/1‑149 (MQ=255)
                                                                  tCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGt                                                                                    <  2:429300/149‑1 (MQ=255)
                                                                  tCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGt                                                                                    >  1:1009805/1‑149 (MQ=255)
                                                                    ccGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGttt                                                                                  <  1:4360169/149‑1 (MQ=255)
                                                                     cGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTa                                                                                 <  1:1404480/149‑1 (MQ=255)
                                                                     cGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTa                                                                                 >  1:414451/1‑149 (MQ=255)
                                                                      gATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATTAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTAt                                                                                >  1:2495027/1‑149 (MQ=255)
                                                                        ttCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATgg                                                                              >  2:140741/1‑149 (MQ=255)
                                                                         tCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATggg                                                                             <  1:3487858/149‑1 (MQ=255)
                                                                          cGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGc                                                                            <  1:4429135/149‑1 (MQ=255)
                                                                             tGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGac                                                                         >  2:4730644/1‑149 (MQ=255)
                                                                                 ggcgAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGAGTAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCGGTGCGGAGTGTTTATGGGCGACACtt                                                                     >  2:4785832/1‑149 (MQ=255)
                                                                                  gcgAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTc                                                                    >  1:4705987/1‑149 (MQ=255)
                                                                                   cgAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCg                                                                   <  2:1247092/149‑1 (MQ=255)
                                                                                    gAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGt                                                                  <  2:3882895/149‑1 (MQ=255)
                                                                                     aaaGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGtt                                                                 >  2:1196004/1‑149 (MQ=255)
                                                                                      aaGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTc                                                                >  1:3283102/1‑149 (MQ=255)
                                                                                       aGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTcc                                                               <  2:1646460/149‑1 (MQ=255)
                                                                                       aGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTcc                                                               >  1:4568613/1‑149 (MQ=255)
                                                                                          ggCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGc                                                            <  2:3538579/149‑1 (MQ=255)
                                                                                          ggCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGc                                                            >  1:920954/1‑149 (MQ=255)
                                                                                           gCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCt                                                           >  1:365597/1‑149 (MQ=255)
                                                                                              cTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTaa                                                        >  2:1881416/1‑149 (MQ=255)
                                                                                              cTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTaa                                                        <  2:2439640/149‑1 (MQ=255)
                                                                                               tGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAAt                                                       <  2:3749003/149‑1 (MQ=255)
                                                                                                gATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATg                                                      >  2:2857799/1‑149 (MQ=255)
                                                                                                     ccTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGAACTg                                                 >  1:233970/1‑149 (MQ=255)
                                                                                                      cTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTgg                                                >  1:1200106/1‑149 (MQ=255)
                                                                                                      cTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTgg                                                <  2:2618791/149‑1 (MQ=255)
                                                                                                       gtGTCGGGGATAATCGCCACGATCTCGGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGt                                               <  1:4973379/148‑1 (MQ=255)
                                                                                                         gTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGc                                             >  2:3183344/1‑149 (MQ=255)
                                                                                                            ggTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGa                                          >  2:692920/1‑149 (MQ=255)
                                                                                                             gTGATAATCGCCACGATCTTTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGaa                                         >  2:3820041/1‑149 (MQ=255)
                                                                                                             gTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGaa                                         >  1:1509275/1‑149 (MQ=255)
                                                                                                              tGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAAt                                        <  1:2913298/149‑1 (MQ=255)
                                                                                                              tGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAAt                                        >  1:319703/1‑149 (MQ=255)
                                                                                                              tGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAAt                                        <  2:1998656/149‑1 (MQ=255)
                                                                                                                aTAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGt                                      <  2:1901349/149‑1 (MQ=255)
                                                                                                                  aaTCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTg                                    >  1:4297556/1‑149 (MQ=255)
                                                                                                                   aTCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTgg                                   <  1:3836479/149‑1 (MQ=255)
                                                                                                                      gCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGa                                <  1:4063858/149‑1 (MQ=255)
                                                                                                                          cGATCTCTGTCACCTCCGGTGTCGATTAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTAt                            >  2:820640/1‑149 (MQ=255)
                                                                                                                          cGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTAt                            <  2:4329905/149‑1 (MQ=255)
                                                                                                                           gATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTAtg                           >  2:3883066/1‑149 (MQ=255)
                                                                                                                             tctcTGTCACCTCCGGTGTCGATTAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTAtgtg                         <  1:3805751/149‑1 (MQ=255)
                                                                                                                              ctctGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGa                        >  2:182331/1‑149 (MQ=255)
                                                                                                                                ctGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAAt                      <  2:2191979/149‑1 (MQ=255)
                                                                                                                                   tCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGc                   <  1:1641881/149‑1 (MQ=255)
                                                                                                                                   tCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGc                   <  1:760584/149‑1 (MQ=255)
                                                                                                                                     aCCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCtt                 >  2:911006/1‑149 (MQ=255)
                                                                                                                                      ccTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCttt                >  1:3786785/1‑149 (MQ=255)
                                                                                                                                      ccTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCttt                <  2:4583991/149‑1 (MQ=255)
                                                                                                                                      ccTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCttt                >  1:2448518/1‑149 (MQ=255)
                                                                                                                                       cTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCTTTa               >  2:3291999/1‑149 (MQ=255)
                                                                                                                                            gtgtCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCTTTATGGGc          <  2:3656554/149‑1 (MQ=255)
                                                                                                                                             tgtCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCTTTATGGGCa         <  1:1906092/149‑1 (MQ=255)
                                                                                                                                             tgtCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCTTTATGGGCa         <  2:2344320/149‑1 (MQ=255)
                                                                                                                                              gtCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCTTTATGGGCAt        >  2:1804391/1‑149 (MQ=255)
                                                                                                                                               tCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCTTTATGGGCATg       <  1:2824968/149‑1 (MQ=255)
                                                                                                                                                cGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCTTTATGGGCATGa      >  2:3641426/1‑149 (MQ=255)
                                                                                                                                                 gATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCTTTATGGGCATGAc     <  1:2510073/149‑1 (MQ=255)
                                                                                                                                                 gATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCTTTATGGGCATGAc     >  2:2243024/1‑149 (MQ=255)
                                                                                                                                                   tAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCTTTATGGGCATGACGc   <  1:5160846/149‑1 (MQ=255)
                                                                                                                                                    aaGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCTTTATGGGCATGACGCt  <  2:3595230/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGCCTTGTCGGTGATAATCGCCACGATCTCTGTCACCTCCGGTGTCGATAAGGGCATTCGCGTGTTATCGGAGCTTAATGTCGCGCTGGCGCTGGGATTGATCCTGTTCGTATTGTTTATGGGCGACACTTCGTTCCTGCTTAATGCACTGGTGCTGAATGTTGGCGACTATGTGAATCGCTTTATGGGCATGACGCT  >  CP009273/325788‑326084

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: