Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,835,413 C→A 7.2% I355I (ATC→ATA ascB → cryptic 6‑phospho‑beta‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,835,4130CA 7.2% 1417.7 / 56.4 374I355I (ATC→ATAascBcryptic 6‑phospho‑beta‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (172/175);  new base A (17/10);  total (189/185)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 375 total aligned reads displayed.

TCTCTTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGC  >  CP009273/2835265‑2835560
                                                                                                                                                    |                                                                                                                                                   
 ctctTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTTAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCa                                                                                                                                                    <  1:1835825/149‑1 (MQ=255)
  tctTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGACTTGATCCACTCGGTCTGCGTATCAc                                                                                                                                                   >  2:3660081/1‑149 (MQ=255)
   ctTTCAGCTTTTACGCCTGGCGCTGCGCCTGGGCGGAGATTAACGCCAACAACTGCTGTGGGGCGAAGGTGGTGAAATCGCTGCGTAATCCGTATCTACTGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCAcc                                                                                                                                                  <  1:2269448/149‑1 (MQ=255)
   ctTTCAGCTATTCCGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCAcc                                                                                                                                                  <  2:2253616/149‑1 (MQ=255)
   ctTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCAcc                                                                                                                                                  >  1:1061971/1‑149 (MQ=255)
    tttCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCa                                                                                                                                                 <  2:43113/149‑1 (MQ=255)
     ttCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCAt                                                                                                                                                <  2:3977871/149‑1 (MQ=255)
           tATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAatat                                                                                                                                          <  1:3045534/149‑1 (MQ=255)
            aTTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATatg                                                                                                                                         >  2:4666530/1‑149 (MQ=255)
            aTTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATatg                                                                                                                                         >  2:3932486/1‑149 (MQ=255)
             ttACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGGATatga                                                                                                                                        >  2:4379863/1‑149 (MQ=255)
             ttACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATatga                                                                                                                                        >  2:111165/1‑149 (MQ=255)
             ttACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATatga                                                                                                                                        <  2:697459/149‑1 (MQ=255)
             ttACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATatga                                                                                                                                        <  2:1947126/149‑1 (MQ=255)
             ttACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATatga                                                                                                                                        >  1:3508499/1‑149 (MQ=255)
                cGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGt                                                                                                                                     <  1:2886008/149‑1 (MQ=255)
                   cTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACg                                                                                                                                  >  2:5239395/1‑149 (MQ=255)
                    tCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGa                                                                                                                                 >  1:1670876/1‑149 (MQ=255)
                      gcgcTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATc                                                                                                                               <  1:18822/149‑1 (MQ=255)
                      gcgcTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATc                                                                                                                               <  1:4493146/149‑1 (MQ=255)
                         cTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGtt                                                                                                                            <  1:1945288/149‑1 (MQ=255)
                         cTGCGCCTCGGCGGAGATGAACGCAAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGtt                                                                                                                            <  1:4561654/149‑1 (MQ=255)
                            cgcCTCGGCGGAGATGAACGCCAACAACAGCAGTGTGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATc                                                                                                                         <  2:3560968/149‑1 (MQ=255)
                             gcCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCa                                                                                                                        >  1:834167/1‑149 (MQ=255)
                                  ggcggAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGc                                                                                                                   >  1:1248386/1‑149 (MQ=255)
                                  ggcggAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGc                                                                                                                   <  1:2731047/149‑1 (MQ=255)
                                  ggcggAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGc                                                                                                                   <  1:4678810/149‑1 (MQ=255)
                                  ggcggAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGc                                                                                                                   >  1:3866435/1‑149 (MQ=255)
                                  ggcggAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGc                                                                                                                   <  1:3766123/149‑1 (MQ=255)
                                          tGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATAGAACCACTCGGTCAGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGttt                                                                                                           >  1:217485/1‑149 (MQ=255)
                                          tGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGCGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCAGGAATATGATGTACGATCGTTATCAGAAGCCGCTGatt                                                                                                           >  2:5021365/1‑149 (MQ=255)
                                           gAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTc                                                                                                          <  2:5163166/149‑1 (MQ=255)
                                           gAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTc                                                                                                          <  1:114318/149‑1 (MQ=255)
                                           gAACGCAAACAACAGCTGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTc                                                                                                          <  1:2739483/149‑1 (MQ=255)
                                            aaCGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCt                                                                                                         <  2:2929510/149‑1 (MQ=255)
                                            aaCGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCt                                                                                                         >  2:2735234/1‑149 (MQ=255)
                                             aCGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCtg                                                                                                        <  1:2182793/149‑1 (MQ=255)
                                             aCGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCtg                                                                                                        >  1:2969353/1‑149 (MQ=255)
                                                ccAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCtggtg                                                                                                     >  2:2970467/1‑149 (MQ=255)
                                                ccAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCtggtg                                                                                                     >  2:2583961/1‑149 (MQ=255)
                                                ccAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCtggtg                                                                                                     <  2:144283/149‑1 (MQ=255)
                                                  aacaacAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGa                                                                                                   >  2:3205360/1‑149 (MQ=255)
                                                  aacaacAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGa                                                                                                   >  2:1180176/1‑149 (MQ=255)
                                                   acaacaGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGaa                                                                                                  <  1:4356131/149‑1 (MQ=255)
                                                   acaacaGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGaa                                                                                                  >  1:505639/1‑149 (MQ=255)
                                                    caacaGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGaaa                                                                                                 <  2:111959/149‑1 (MQ=255)
                                                      acaGCAGTGCGGCTAACGTGGTGATATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCTGTCTGCGTATAACAATTAATATGATGTACGATCTTTATCATAAGCCGCTGTTTCTGTTGGAAAAc                                                                                               <  2:1992234/149‑1 (MQ=255)
                                                      acaGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAAc                                                                                               <  1:3829337/149‑1 (MQ=255)
                                                      acaGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAAc                                                                                               >  2:4546471/1‑149 (MQ=255)
                                                      acaGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAAc                                                                                               <  1:3240490/149‑1 (MQ=255)
                                                        agcagTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACgg                                                                                             <  2:5119520/149‑1 (MQ=255)
                                                          ctgTGCGGCTAACGGGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGcc                                                                                           <  1:4445567/147‑1 (MQ=255)
                                                          cagTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGcc                                                                                           <  1:278542/149‑1 (MQ=255)
                                                           agTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCt                                                                                          <  2:3380521/149‑1 (MQ=255)
                                                               cggcgAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGgc                                                                                      <  1:2583961/149‑1 (MQ=255)
                                                                    aaCGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCaaa                                                                                 >  2:3264052/1‑149 (MQ=255)
                                                                          gtgAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGa                                                                           >  2:4490519/1‑149 (MQ=255)
                                                                           tgAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGaa                                                                          <  2:1452599/149‑1 (MQ=255)
                                                                           tgAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGaa                                                                          >  1:3154613/1‑149 (MQ=255)
                                                                             aaaTCGCTGCGTAATCCTTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAAtt                                                                        <  1:561928/149‑1 (MQ=255)
                                                                              aaTCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAAttt                                                                       >  2:2176320/1‑149 (MQ=255)
                                                                                    tGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                                                                 >  2:2359125/1‑149 (MQ=255)
                                                                                     gCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCa                                                                >  1:631304/1‑149 (MQ=255)
                                                                                       gTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAAt                                                              <  1:346113/149‑1 (MQ=255)
                                                                                          aTCCGTATCTACTGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGc                                                           <  1:4063542/149‑1 (MQ=255)
                                                                                            ccGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCAGGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCga                                                         >  2:2049233/1‑149 (MQ=255)
                                                                                              gTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCgaga                                                       >  1:909961/1‑149 (MQ=255)
                                                                                               tATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGAt                                                      >  1:3426590/1‑149 (MQ=255)
                                                                                                aTCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGAtt                                                     >  2:4088959/1‑149 (MQ=255)
                                                                                                  cTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTaa                                                   <  1:4978217/149‑1 (MQ=255)
                                                                                                   tACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAac                                                  >  1:2683071/1‑149 (MQ=255)
                                                                                                    aCAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacg                                                 >  1:3970689/1‑149 (MQ=255)
                                                                                                       ggTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacgacg                                              >  1:3293764/1‑149 (MQ=255)
                                                                                                       ggTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacgacg                                              <  1:985826/149‑1 (MQ=255)
                                                                                                             cGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATc                                        <  2:2490198/149‑1 (MQ=255)
                                                                                                             cGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATc                                        <  2:4231922/149‑1 (MQ=255)
                                                                                                                ctggggctggggAATTGCGCCCCTCGGTCTGCGTATCACCATGAAGATGATGTACGCTCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCa                                     <  2:4576894/149‑1 (MQ=255)
                                                                                                                 tggggctggggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCAt                                    <  2:573666/149‑1 (MQ=255)
                                                                                                                  ggggctggggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATc                                   >  2:4795392/1‑149 (MQ=255)
                                                                                                                  ggggctggggAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATc                                   >  2:1806867/1‑149 (MQ=255)
                                                                                                                     gctggggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGc                                >  1:161861/1‑149 (MQ=255)
                                                                                                                       tggggAATTGATCCACTCGGTCTGCGTATCACCATGACTATGATGTACGATCGGTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTa                              <  2:508571/149‑1 (MQ=255)
                                                                                                                       tggggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTCTCGCATCAGCTa                              >  2:1725766/1‑149 (MQ=255)
                                                                                                                        ggggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAACGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTAc                             <  1:1777620/149‑1 (MQ=255)
                                                                                                                         gggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACt                            <  2:1546448/149‑1 (MQ=255)
                                                                                                                         gggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACt                            >  2:536453/1‑149 (MQ=255)
                                                                                                                          gtcaTTGATCCTCTCGGTCTGCGTATAACAATGATTTTGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGGGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACtt                           <  1:4154438/146‑1 (MQ=255)
                                                                                                                           gAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTa                          <  1:242316/149‑1 (MQ=255)
                                                                                                                            aaTTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTAc                         >  1:2054251/1‑149 (MQ=255)
                                                                                                                                 aTCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGaa                    >  1:3200821/1‑149 (MQ=255)
                                                                                                                                 aTCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGaa                    >  2:3851958/1‑149 (MQ=255)
                                                                                                                                 aTCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGaa                    >  2:3110757/1‑149 (MQ=255)
                                                                                                                                     aCTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACata                <  1:3609990/149‑1 (MQ=255)
                                                                                                                                      cTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACatat               >  2:2853903/1‑149 (MQ=255)
                                                                                                                                       tCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATc              <  2:1314048/149‑1 (MQ=255)
                                                                                                                                        cGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATcc             <  1:2824407/149‑1 (MQ=255)
                                                                                                                                           tCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCcgcg          >  1:759520/1‑149 (MQ=255)
                                                                                                                                              gCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCaa       >  1:1424046/1‑149 (MQ=255)
                                                                                                                                               cGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCTCATCAGCTACTTACGCGAACATATCCGCGCAAt      >  1:4830721/1‑149 (MQ=255)
                                                                                                                                                gTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATg     >  2:1356521/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCTCTTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGC  >  CP009273/2835265‑2835560

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: