Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,374,303 T→C 5.7% R721R (CGT→CGC ycjT → putative family 65 glycosyl hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,374,3030TC 5.7% 832.1 / 28.0 250R721R (CGT→CGCycjTputative family 65 glycosyl hydrolase
Reads supporting (aligned to +/- strand):  ref base T (111/124);  new base C (7/7);  total (118/132)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 253 total aligned reads displayed.

GTGCGTGACGGTGAATTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGGATGGTGTAATCACCGA  >  CP009273/1374154‑1374451
                                                                                                                                                     |                                                                                                                                                    
 tgcgtgACGGTGAATTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGtt                                                                                                                                                     <  1:2329672/150‑1 (MQ=255)
   cgtgACGGTGAATTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTATGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTg                                                                                                                                                   <  1:2309580/150‑1 (MQ=255)
    gtgACGGTGAATTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGa                                                                                                                                                  <  1:153492/150‑1 (MQ=255)
     tgACGGTGAATTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGaa                                                                                                                                                 <  2:1585064/150‑1 (MQ=255)
     tgACGGTGAATTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGaa                                                                                                                                                 <  2:1585634/150‑1 (MQ=255)
           tGAATTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCa                                                                                                                                           >  2:2478593/1‑150 (MQ=255)
             aaTTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGc                                                                                                                                         >  1:3215958/1‑150 (MQ=255)
             aaTTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGc                                                                                                                                         >  2:2667289/1‑150 (MQ=255)
              aTTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCt                                                                                                                                        <  1:662930/150‑1 (MQ=255)
               ttGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCtt                                                                                                                                       >  1:2699673/1‑150 (MQ=255)
               ttGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCtt                                                                                                                                       >  1:1499634/1‑150 (MQ=255)
               ttGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCtt                                                                                                                                       <  1:4572782/150‑1 (MQ=255)
                tGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGCTTGAACGGTCAGCTta                                                                                                                                      >  1:1288891/1‑150 (MQ=255)
                  cATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTtata                                                                                                                                    >  2:574095/1‑150 (MQ=255)
                  cATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGCTTGAACGGTCAGCTtata                                                                                                                                    <  1:4426567/150‑1 (MQ=255)
                     ctcAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAAcc                                                                                                                                 >  2:907044/1‑150 (MQ=255)
                      tcAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTtata                                                                                                                                    <  1:3418103/146‑1 (MQ=255)
                       cAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGt                                                                                                                               <  2:4353893/150‑1 (MQ=255)
                        aaTCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTg                                                                                                                              >  2:3496754/1‑150 (MQ=255)
                           ccGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCt                                                                                                                           <  2:6520/150‑1 (MQ=255)
                           ccGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCt                                                                                                                           <  1:1284738/150‑1 (MQ=255)
                             ggCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTga                                                                                                                         <  1:4138655/150‑1 (MQ=255)
                             ggCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTga                                                                                                                         >  1:4005602/1‑150 (MQ=255)
                              gCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTgaa                                                                                                                        <  2:1688501/150‑1 (MQ=255)
                                  tACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAAt                                                                                                                    <  1:3888429/150‑1 (MQ=255)
                                     cTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAAtctg                                                                                                                 <  2:4362191/150‑1 (MQ=255)
                                     cTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAAtctg                                                                                                                 <  1:4432896/150‑1 (MQ=255)
                                      tGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAAtctgt                                                                                                                <  1:926613/150‑1 (MQ=255)
                                       gAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAAtctgtt                                                                                                               >  2:3086853/1‑150 (MQ=255)
                                        aGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAAtctgttt                                                                                                              <  2:3310967/150‑1 (MQ=255)
                                        aGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAAtctgttt                                                                                                              >  1:2672903/1‑150 (MQ=255)
                                           aGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAAtctgttttct                                                                                                           <  2:3206988/150‑1 (MQ=255)
                                           aGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAAtctgttttct                                                                                                           >  1:995185/1‑150 (MQ=255)
                                                cAACAGTTGTCTTTCCCTCTGTTCTGGCAGGTCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTg                                                                                                      <  1:4178/150‑1 (MQ=255)
                                                cAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTg                                                                                                      >  2:2638465/1‑150 (MQ=255)
                                                cAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTg                                                                                                      >  2:2151830/1‑150 (MQ=255)
                                                cAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTg                                                                                                      >  1:2515745/1‑150 (MQ=255)
                                                  aCAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTTTTTTCTGTTTtgg                                                                                                    >  1:2458124/1‑150 (MQ=255)
                                                     gTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGa                                                                                                 >  2:4582732/1‑150 (MQ=255)
                                                     gTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGa                                                                                                 <  1:858589/150‑1 (MQ=255)
                                                       tGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGatt                                                                                               >  1:4581369/1‑150 (MQ=255)
                                                       tGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGagt                                                                                               >  2:464411/1‑148 (MQ=255)
                                                        gTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGattt                                                                                              <  2:3578/150‑1 (MQ=255)
                                                          cTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGatttta                                                                                            >  1:1128653/1‑150 (MQ=255)
                                                           tttCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGattttat                                                                                           <  1:3144763/150‑1 (MQ=255)
                                                            gtCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGattttatt                                                                                          <  1:2491945/149‑1 (MQ=255)
                                                                ctctGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGc                                                                                      <  1:539171/150‑1 (MQ=255)
                                                                     ttCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTc                                                                                 <  1:4614889/150‑1 (MQ=255)
                                                                        tGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAAt                                                                              >  1:3081510/1‑150 (MQ=255)
                                                                         ggCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATg                                                                             >  1:388057/1‑150 (MQ=255)
                                                                          gCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATgg                                                                            >  2:1970247/1‑150 (MQ=255)
                                                                           cAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATggg                                                                           >  1:1380659/1‑150 (MQ=255)
                                                                           cAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATggg                                                                           <  2:3685771/150‑1 (MQ=255)
                                                                            aGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGa                                                                          <  2:4208932/150‑1 (MQ=255)
                                                                            aGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATGTTGCCCTTCAATGGGa                                                                          >  1:1027366/1‑150 (MQ=255)
                                                                                 tGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCt                                                                     >  1:4475568/1‑150 (MQ=255)
                                                                                  gCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTa                                                                    <  1:951246/150‑1 (MQ=255)
                                                                                       ttACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCaaa                                                               >  2:4424339/1‑150 (MQ=255)
                                                                                          cAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACAt                                                            >  2:3180881/1‑150 (MQ=255)
                                                                                            ggTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCa                                                          <  1:669852/150‑1 (MQ=255)
                                                                                              tCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCagg                                                        <  2:3991657/150‑1 (MQ=255)
                                                                                                 ctctTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCaggagg                                                     <  2:3956444/150‑1 (MQ=255)
                                                                                                      gACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGatgaa                                                >  2:1311791/1‑150 (MQ=255)
                                                                                                      gACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTGGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGatgaa                                                <  2:3874477/150‑1 (MQ=255)
                                                                                                       aCGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGatgaat                                               <  2:457903/150‑1 (MQ=255)
                                                                                                       aCGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGatgaat                                               <  1:4333530/150‑1 (MQ=255)
                                                                                                          cgcAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGAtgaatgaa                                            <  1:4381841/150‑1 (MQ=255)
                                                                                                             aGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACt                                         >  2:397645/1‑150 (MQ=255)
                                                                                                               cGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGc                                       >  2:4707480/1‑150 (MQ=255)
                                                                                                                    tGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAggg                                  <  1:1993095/150‑1 (MQ=255)
                                                                                                                       gATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTa                               >  2:2055444/1‑150 (MQ=255)
                                                                                                                       gATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTa                               >  2:4482021/1‑150 (MQ=255)
                                                                                                                       gATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTa                               >  2:1943051/1‑150 (MQ=255)
                                                                                                                        aTTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTaa                              <  1:4249110/150‑1 (MQ=255)
                                                                                                                        aTTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGACTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTaa                              >  1:1766907/1‑150 (MQ=255)
                                                                                                                          tCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAAtt                            <  1:3448937/150‑1 (MQ=255)
                                                                                                                           cGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAAttt                           <  1:191835/150‑1 (MQ=255)
                                                                                                                            gAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCACGGGGTAAtatt                          >  1:8434/1‑150 (MQ=255)
                                                                                                                              aCTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCg                        <  2:2326354/150‑1 (MQ=255)
                                                                                                                              aCTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCg                        <  1:1665401/150‑1 (MQ=255)
                                                                                                                                ttCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGAt                      <  1:3313715/150‑1 (MQ=255)
                                                                                                                                ttCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGAt                      >  2:619907/1‑150 (MQ=255)
                                                                                                                                ttCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGAt                      <  2:1640640/150‑1 (MQ=255)
                                                                                                                                   tGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTg                   >  1:3928787/1‑150 (MQ=255)
                                                                                                                                    gcgcCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTgg                  <  1:1348655/150‑1 (MQ=255)
                                                                                                                                     cgcCCGTTTCACTGCGCTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGGa                 <  2:660049/150‑1 (MQ=255)
                                                                                                                                      gcCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGGAt                >  1:21073/1‑150 (MQ=255)
                                                                                                                                       cccGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGGATg               >  2:3734402/1‑150 (MQ=255)
                                                                                                                                           tttCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGGATGgtgt           <  2:4153027/150‑1 (MQ=255)
                                                                                                                                            ttCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGGATGGTGTa          <  2:2319051/150‑1 (MQ=255)
                                                                                                                                              cACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGGATGGTGTAAt        <  2:1137035/150‑1 (MQ=255)
                                                                                                                                              cACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGGATGGTGTAAt        <  2:3539334/150‑1 (MQ=255)
                                                                                                                                              cACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGGATGGTGTAAt        <  2:3460698/150‑1 (MQ=255)
                                                                                                                                              cACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGAATGGTGTAAt        <  2:1137129/150‑1 (MQ=255)
                                                                                                                                               aCTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGGATGGTGTAATc       <  2:1176411/150‑1 (MQ=255)
                                                                                                                                               aCTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGGATGGTGTAATc       >  1:2059858/1‑150 (MQ=255)
                                                                                                                                                   cGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGGATGGTGTAATCACCg   >  1:3501876/1‑150 (MQ=255)
                                                                                                                                                   cGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGGATGGTGTAATCACCg   >  2:3654712/1‑150 (MQ=255)
                                                                                                                                                    gTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGGATGGTGTAATCACCGa  >  1:3730794/1‑150 (MQ=255)
                                                                                                                                                    gTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGGATGGTGTAATCACCGa  <  1:568575/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTGCGTGACGGTGAATTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCAGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTTTCACTGCGTTTGAACGGTCAGCTTATAACCGTGGCTGAAGAATCTGTTTTCTGTTTGGGTGATTTTATTTTGCCCTTCAATGGGACCGCTACCAAACATCAGGAGGATGAATGAAACTGCAAGGGGTAATTTTCGATCTGGATGGTGTAATCACCGA  >  CP009273/1374154‑1374451

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: