Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 456,588 G→T 88.6% D749Y (GAC→TAC)  lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273456,5880GT88.6% 701.1 / 84.7 262D749Y (GAC→TAC) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base G (20/10);  new base T (119/113);  total (139/123)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 262 total aligned reads displayed.

GATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGG  >  CP009273/456441‑456736
                                                                                                                                                   |                                                                                                                                                    
gATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTAc                                                                                                                                                    >  2:3273170/1‑150 (MQ=255)
 aTGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACa                                                                                                                                                   <  1:3393403/150‑1 (MQ=255)
    tgGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACg                                                                                                                                                >  1:1239425/1‑150 (MQ=255)
    tgGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACg                                                                                                                                                >  1:1241433/1‑150 (MQ=255)
       cAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTaa                                                                                                                                             <  2:1351137/150‑1 (MQ=255)
       cAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTaa                                                                                                                                             <  1:97006/150‑1 (MQ=255)
       aaaTGACCGGTGAGATCTCTCTGCGTTGTAAGGTACTTCCGATAGGTGGTTTGATAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAATACAGTGCTAATTCCTTTCGAAAATAAACGCGATCTGGAAGAGATTCATTTCAACGTaa                                                                                                                                             <  2:3618019/149‑1 (MQ=255)
         aTGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAAtt                                                                                                                                           <  2:2137261/150‑1 (MQ=255)
         aTGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAAtt                                                                                                                                           <  1:3622199/150‑1 (MQ=255)
          tGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTg                                                                                                                                          <  1:668858/150‑1 (MQ=255)
          tGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTg                                                                                                                                          >  1:4486888/1‑150 (MQ=255)
             ccGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCg                                                                                                                                       >  2:413950/1‑150 (MQ=255)
             ccGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCg                                                                                                                                       >  2:411216/1‑150 (MQ=255)
                  gagaTCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTg                                                                                                                                  >  1:3172799/1‑150 (MQ=255)
                   agaTCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTgg                                                                                                                                 >  1:1362515/1‑150 (MQ=255)
                      tCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACa                                                                                                                              >  1:946606/1‑150 (MQ=255)
                        aCTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACAtt                                                                                                                            >  1:4228522/1‑150 (MQ=255)
                         ctctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTc                                                                                                                           >  1:1783513/1‑150 (MQ=255)
                         ctctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTc                                                                                                                           >  1:4036930/1‑150 (MQ=255)
                         ctctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTc                                                                                                                           <  1:3029892/150‑1 (MQ=255)
                           ctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCAt                                                                                                                         <  2:3485696/150‑1 (MQ=255)
                           ctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCAt                                                                                                                         <  1:1167977/150‑1 (MQ=255)
                             gCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAACACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATcc                                                                                                                       <  2:1134622/150‑1 (MQ=255)
                               gTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCtg                                                                                                                     >  1:2611728/1‑150 (MQ=255)
                                tGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCtgt                                                                                                                    <  2:829981/150‑1 (MQ=255)
                                 ggTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCtgtg                                                                                                                   >  2:2895634/1‑150 (MQ=255)
                                   tCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGaa                                                                                                                 >  2:3095825/1‑150 (MQ=255)
                                   tCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGaa                                                                                                                 <  2:3140075/150‑1 (MQ=255)
                                    cAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAg                                                                                                                <  2:4244033/150‑1 (MQ=255)
                                     aGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAgc                                                                                                               <  1:380430/150‑1 (MQ=255)
                                      ggTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGACAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAgcg                                                                                                              >  1:3704792/1‑150 (MQ=255)
                                        tACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCa                                                                                                            >  2:2649971/1‑150 (MQ=255)
                                          cTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCAtt                                                                                                          <  2:2599020/150‑1 (MQ=255)
                                           tGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTg                                                                                                         <  2:1503242/150‑1 (MQ=255)
                                           tGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTg                                                                                                         <  2:1838210/150‑1 (MQ=255)
                                                 tCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAg                                                                                                   >  2:2812142/1‑150 (MQ=255)
                                                   ggtggtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGtt                                                                                                 <  1:1656101/150‑1 (MQ=255)
                                                      ggtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTg                                                                                              <  2:3927978/150‑1 (MQ=255)
                                                      ggtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTg                                                                                              <  2:3317954/150‑1 (MQ=255)
                                                       gtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGa                                                                                             >  1:80852/1‑150 (MQ=255)
                                                        tttGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGACGTTCTGAc                                                                                            >  2:3900654/1‑150 (MQ=255)
                                                           gaaagaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGActct                                                                                         <  1:1373890/150‑1 (MQ=255)
                                                                aaaaaCTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGACGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGgcgc                                                                                    >  2:3172009/1‑150 (MQ=255)
                                                                aaaaaCTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGgcgc                                                                                    >  2:887130/1‑150 (MQ=255)
                                                                 aaaaCTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCt                                                                                   >  2:3784255/1‑150 (MQ=255)
                                                                   aaCTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGc                                                                                 >  1:1444795/1‑150 (MQ=255)
                                                                     cTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCaa                                                                               >  1:2623134/1‑150 (MQ=255)
                                                                     cTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCaa                                                                               >  1:324686/1‑150 (MQ=255)
                                                                        cTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAAt                                                                            <  2:4119290/150‑1 (MQ=255)
                                                                           gCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGaa                                                                         <  1:2710582/150‑1 (MQ=255)
                                                                            cAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAAc                                                                        <  2:1019729/150‑1 (MQ=255)
                                                                            cAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAAc                                                                        <  1:229012/150‑1 (MQ=255)
                                                                             aGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAAcc                                                                       >  2:1657119/1‑150 (MQ=255)
                                                                                 cATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATGGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCt                                                                   >  1:1163717/1‑150 (MQ=255)
                                                                                   tCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTgg                                                                 >  2:1900148/1‑150 (MQ=255)
                                                                                   tCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTgg                                                                 <  2:1944783/150‑1 (MQ=255)
                                                                                       ggcggGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATg                                                             <  2:2611728/150‑1 (MQ=255)
                                                                                        gcggGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGc                                                            >  1:2339179/1‑150 (MQ=255)
                                                                                         cggGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCa                                                           >  2:1086166/1‑150 (MQ=255)
                                                                                         cggGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCa                                                           <  1:2200177/150‑1 (MQ=255)
                                                                                             aTTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGtt                                                       >  1:1520251/1‑150 (MQ=255)
                                                                                                aaaaCAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTtgtg                                                    >  2:4323997/1‑150 (MQ=255)
                                                                                                 aaaCAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGa                                                   <  2:2761453/150‑1 (MQ=255)
                                                                                                   aCAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACt                                                 >  1:1420735/1‑150 (MQ=255)
                                                                                                   aCAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACt                                                 >  1:3001104/1‑150 (MQ=255)
                                                                                                    cAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTg                                                <  2:4334120/150‑1 (MQ=255)
                                                                                                    cAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTg                                                >  2:3299948/1‑150 (MQ=255)
                                                                                                     aGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGc                                               >  1:1218259/1‑150 (MQ=255)
                                                                                                        gCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCaaa                                            >  1:4269471/1‑150 (MQ=255)
                                                                                                          tAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAAt                                          <  1:3506783/150‑1 (MQ=255)
                                                                                                          tAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAAt                                          <  1:3066080/150‑1 (MQ=255)
                                                                                                           aaTTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATa                                         <  2:2234201/150‑1 (MQ=255)
                                                                                                                    cGAAAATAAACGCGATCTGGAAGAGAGTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTcg                                <  2:3466222/150‑1 (MQ=255)
                                                                                                                       aaaTAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGACGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCgcgc                             >  2:3217972/1‑150 (MQ=255)
                                                                                                                       aaaTAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCgcgc                             >  2:3218011/1‑150 (MQ=255)
                                                                                                                        aaTAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCa                            <  1:2514730/150‑1 (MQ=255)
                                                                                                                         aTAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCaa                           <  1:2029867/150‑1 (MQ=255)
                                                                                                                          tAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCCTTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCaaa                          >  2:4251166/1‑150 (MQ=255)
                                                                                                                          tAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCaaa                          >  2:202286/1‑150 (MQ=255)
                                                                                                                          tAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCaaa                          >  2:4055855/1‑150 (MQ=255)
                                                                                                                           aaaCGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCaaaa                         <  1:4335781/150‑1 (MQ=255)
                                                                                                                            aaCGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAAt                        >  1:3262708/1‑150 (MQ=255)
                                                                                                                                cgATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCAc                    >  2:33998/1‑150 (MQ=255)
                                                                                                                                cgATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCAc                    <  1:1936848/150‑1 (MQ=255)
                                                                                                                                cgATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCAc                    >  1:1600293/1‑150 (MQ=255)
                                                                                                                                 gATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACt                   >  2:3823877/1‑150 (MQ=255)
                                                                                                                                 gATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACt                   <  2:4461512/150‑1 (MQ=255)
                                                                                                                                  aTCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACta                  >  1:1134859/1‑150 (MQ=255)
                                                                                                                                   tCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACtaa                 >  2:279532/1‑150 (MQ=255)
                                                                                                                                    cTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACtaat                <  1:413927/150‑1 (MQ=255)
                                                                                                                                    cTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACtaat                <  1:413990/150‑1 (MQ=255)
                                                                                                                                    cTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACtaat                >  1:605395/1‑150 (MQ=255)
                                                                                                                                        aaGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATaaaa            >  2:2292439/1‑150 (MQ=255)
                                                                                                                                        aaGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATaaaa            <  2:79487/150‑1 (MQ=255)
                                                                                                                                          gagaTTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAAc          >  1:1300928/1‑150 (MQ=255)
                                                                                                                                           agaTTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACa         >  1:295172/1‑150 (MQ=255)
                                                                                                                                              ttCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATTAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAggg      >  1:2720068/1‑150 (MQ=255)
                                                                                                                                                  tGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTgg  >  2:903531/1‑150 (MQ=255)
                                                                                                                                                  tGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTgg  >  2:904588/1‑150 (MQ=255)
                                                                                                                                                  tGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGAGCAAAATGCACTAATAAAAACAGGGCTgg  >  2:905105/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGG  >  CP009273/456441‑456736

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: