Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,601,871 C→A 16.4% P90Q (CCG→CAG)  tam → trans‑aconitate methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,601,8710CA16.4% 779.4 / 90.6 262P90Q (CCG→CAG) tamtrans‑aconitate methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (83/136);  new base A (14/29);  total (97/165)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 268 total aligned reads displayed.

GGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGACCCAT  >  CP009273/1601723‑1602019
                                                                                                                                                    |                                                                                                                                                    
ggCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAAcc                                                                                                                                                      <  1:1348593/149‑1 (MQ=255)
   ccAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGt                                                                                                                                                   >  2:3625612/1‑149 (MQ=255)
      ggTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACa                                                                                                                                                >  2:1253320/1‑149 (MQ=255)
      ggTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACa                                                                                                                                                >  2:2673074/1‑149 (MQ=255)
       gTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAg                                                                                                                                               >  2:681708/1‑149 (MQ=255)
           cAGCACCGCCCTTCTACAACAGCGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCAc                                                                                                                                           <  2:1575233/149‑1 (MQ=255)
           cAGCACCGCCCTTCTACAACACCGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCAc                                                                                                                                           <  1:2096828/149‑1 (MQ=255)
            aGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACt                                                                                                                                          >  1:912774/1‑149 (MQ=255)
                 cGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATc                                                                                                                                     >  2:2021248/1‑149 (MQ=255)
                 cGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATc                                                                                                                                     <  1:3496132/149‑1 (MQ=255)
                 cGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATc                                                                                                                                     <  1:1836414/149‑1 (MQ=255)
                 cGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATc                                                                                                                                     <  1:313373/149‑1 (MQ=255)
                   cccTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTg                                                                                                                                   >  1:1869656/1‑149 (MQ=255)
                     cTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAt                                                                                                                                 <  1:2914192/149‑1 (MQ=255)
                       tCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAttt                                                                                                                               <  2:1691958/149‑1 (MQ=255)
                        cTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAtttt                                                                                                                              <  2:2385399/149‑1 (MQ=255)
                        cTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAtttt                                                                                                                              <  1:3963146/149‑1 (MQ=255)
                          acaacaACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTg                                                                                                                            >  2:2202898/1‑149 (MQ=255)
                          acaacaACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTg                                                                                                                            <  2:3060168/149‑1 (MQ=255)
                            aacaacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCt                                                                                                                          <  2:3332812/149‑1 (MQ=255)
                             acaacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTCCAGGCACTCGATCTGATTTTTGCTa                                                                                                                         >  2:3815725/1‑149 (MQ=255)
                               aacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAAt                                                                                                                       >  1:2247314/1‑149 (MQ=255)
                                acGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATg                                                                                                                      >  2:12861/1‑149 (MQ=255)
                                  gTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGcc                                                                                                                    <  2:181575/149‑1 (MQ=255)
                                  gTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGcc                                                                                                                    <  1:42574/149‑1 (MQ=255)
                                  gTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGcc                                                                                                                    <  1:1760598/149‑1 (MQ=255)
                                   ttGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCt                                                                                                                   >  1:2041316/1‑149 (MQ=255)
                                   ttGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCt                                                                                                                   <  1:1642928/149‑1 (MQ=255)
                                    tGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTc                                                                                                                  <  1:2455264/149‑1 (MQ=255)
                                      gCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCAc                                                                                                                <  2:1320058/149‑1 (MQ=255)
                                      gCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCAc                                                                                                                <  1:1402518/149‑1 (MQ=255)
                                       ccTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACt                                                                                                               <  2:1774407/149‑1 (MQ=255)
                                            ggCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCTGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAAt                                                                                                          <  2:3348060/149‑1 (MQ=255)
                                            ggCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAAt                                                                                                          <  1:1300493/149‑1 (MQ=255)
                                             gCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATg                                                                                                         <  2:2826689/149‑1 (MQ=255)
                                             gCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATg                                                                                                         <  2:4000576/149‑1 (MQ=255)
                                             gCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATg                                                                                                         <  2:2919598/149‑1 (MQ=255)
                                                aGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCt                                                                                                      <  1:753268/149‑1 (MQ=255)
                                                aGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCt                                                                                                      >  1:1720009/1‑149 (MQ=255)
                                                 ggATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTg                                                                                                     >  2:57015/1‑149 (MQ=255)
                                                           cATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCACACTGCAATGGCTGCCCGACCACt                                                                                           <  2:957635/149‑1 (MQ=255)
                                                            aTCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTa                                                                                          <  1:1251059/149‑1 (MQ=255)
                                                             tCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTAc                                                                                         >  1:1559768/1‑149 (MQ=255)
                                                             tCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTAc                                                                                         <  1:2118827/149‑1 (MQ=255)
                                                             tCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTAc                                                                                         <  2:2477525/149‑1 (MQ=255)
                                                                aCTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGaa                                                                                      >  2:3459434/1‑149 (MQ=255)
                                                                 cTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAAt                                                                                     <  2:3075302/149‑1 (MQ=255)
                                                                  tCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAAtt                                                                                    <  2:2427347/149‑1 (MQ=255)
                                                                      ctcCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGttt                                                                                <  1:1078313/149‑1 (MQ=255)
                                                                      ctcCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGttt                                                                                <  2:1032510/149‑1 (MQ=255)
                                                                       tcCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTc                                                                               <  2:3942474/149‑1 (MQ=255)
                                                                        ccGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTcc                                                                              <  2:764637/149‑1 (MQ=255)
                                                                        ccGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTcc                                                                              <  1:2580815/149‑1 (MQ=255)
                                                                          ggCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTc                                                                            <  1:682364/149‑1 (MQ=255)
                                                                             gatgatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATc                                                                         <  1:3309567/149‑1 (MQ=255)
                                                                              atgatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCt                                                                        >  1:1481021/1‑149 (MQ=255)
                                                                               tgatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTg                                                                       >  2:3809210/1‑149 (MQ=255)
                                                                               tgatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTg                                                                       >  1:3749181/1‑149 (MQ=255)
                                                                                gatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTgg                                                                      >  1:2014337/1‑149 (MQ=255)
                                                                                  ttGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGtt                                                                    >  1:1321558/1‑149 (MQ=255)
                                                                                  ttGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGtt                                                                    >  1:2317596/1‑149 (MQ=255)
                                                                                  tggCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCCTCTGGtt                                                                    <  2:2150708/147‑1 (MQ=255)
                                                                                     cTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCt                                                                 <  1:769363/149‑1 (MQ=255)
                                                                                      tGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCtt                                                                >  1:2370348/1‑149 (MQ=255)
                                                                                       gAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCttt                                                               <  2:1268986/149‑1 (MQ=255)
                                                                                          gcgcgcAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACt                                                            >  1:2781545/1‑149 (MQ=255)
                                                                                          gcgcgcAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACt                                                            <  1:167841/149‑1 (MQ=255)
                                                                                             cgcAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTaa                                                         <  2:2041520/149‑1 (MQ=255)
                                                                                             cgcAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTaa                                                         <  2:2969165/149‑1 (MQ=255)
                                                                                                   gCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCa                                                   >  1:539084/1‑149 (MQ=255)
                                                                                                   gCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCa                                                   >  1:1021311/1‑149 (MQ=255)
                                                                                                      ttGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAggg                                                <  1:3691678/149‑1 (MQ=255)
                                                                                                      ttGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAggg                                                <  1:3497191/149‑1 (MQ=255)
                                                                                                         ccAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGt                                             <  2:1462794/149‑1 (MQ=255)
                                                                                                          cAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTg                                            >  2:2990691/1‑149 (MQ=255)
                                                                                                           aGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCACTGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGc                                           >  2:3843517/1‑149 (MQ=255)
                                                                                                           aGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGc                                           <  2:348240/149‑1 (MQ=255)
                                                                                                              cTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCCGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTgg                                        >  1:3317476/1‑149 (MQ=255)
                                                                                                               tGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGc                                       <  2:3659148/149‑1 (MQ=255)
                                                                                                               tGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGc                                       >  1:627076/1‑149 (MQ=255)
                                                                                                                gCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCa                                      >  1:2961636/1‑149 (MQ=255)
                                                                                                                 ccAGTTCGTGGACGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAg                                     >  2:804895/1‑149 (MQ=255)
                                                                                                                 ccAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAg                                     <  2:3763759/149‑1 (MQ=255)
                                                                                                                     ttCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACa                                 >  1:861405/1‑149 (MQ=255)
                                                                                                                       cGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGa                               <  2:4098989/149‑1 (MQ=255)
                                                                                                                          ggAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGc                            <  2:154498/149‑1 (MQ=255)
                                                                                                                                 gACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATaa                     >  2:2441280/1‑149 (MQ=255)
                                                                                                                                  aCATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAAc                    <  2:912774/149‑1 (MQ=255)
                                                                                                                                   cATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACGGATGCCAGATAAct                   <  1:4071539/149‑1 (MQ=255)
                                                                                                                                   cATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAAct                   <  1:3459434/149‑1 (MQ=255)
                                                                                                                                    aTCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAActg                  >  2:1852566/1‑149 (MQ=255)
                                                                                                                                          aaCTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGa            <  1:3836591/149‑1 (MQ=255)
                                                                                                                                            cTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGc          <  2:476949/149‑1 (MQ=255)
                                                                                                                                            cTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGc          <  2:1139608/149‑1 (MQ=255)
                                                                                                                                               gCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGa       <  2:1047732/149‑1 (MQ=255)
                                                                                                                                               gCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTCCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGa       >  2:255132/1‑149 (MQ=255)
                                                                                                                                                 aaCCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGAcc     <  1:3266399/149‑1 (MQ=255)
                                                                                                                                                  aCCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGAccc    <  1:2673070/149‑1 (MQ=255)
                                                                                                                                                   ccGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGACCCa   >  1:2119817/1‑149 (MQ=255)
                                                                                                                                                    cGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCTCTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGACCCAt  <  2:1055185/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGACCCAT  >  CP009273/1601723‑1602019

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: