Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 456,588 G→T 93.4% 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,5880GT93.4% 891.9 / 47.8 272D749Y (GAC→TAC) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base G (11/7);  new base T (121/133);  total (132/140)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 272 total aligned reads displayed.

CGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCT  >  CP009273/456440‑456734
                                                                                                                                                    |                                                                                                                                                  
     tgGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACg                                                                                                                                              >  2:897232/1‑150 (MQ=255)
     tgGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACg                                                                                                                                              <  2:2288729/150‑1 (MQ=255)
      ggCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGt                                                                                                                                             >  1:779421/1‑150 (MQ=255)
        cAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTaa                                                                                                                                           <  2:255317/150‑1 (MQ=255)
         aaTGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAAt                                                                                                                                          <  2:3081032/150‑1 (MQ=255)
          aTGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAAtt                                                                                                                                         <  1:2022975/150‑1 (MQ=255)
           tGACCGGTGAGAACACTCTGCGTGGACAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGAGCATCGCGGCGGGATTAAAACAGTGCTAATTCAGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTg                                                                                                                                        >  2:2471035/1‑150 (MQ=255)
             aCCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGcc                                                                                                                                      >  1:1210968/1‑150 (MQ=255)
              ccGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCg                                                                                                                                     <  1:1528707/150‑1 (MQ=255)
               cGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGa                                                                                                                                    <  1:182680/150‑1 (MQ=255)
               cGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGa                                                                                                                                    <  1:2557199/150‑1 (MQ=255)
                        cACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACAt                                                                                                                           >  2:3683549/1‑150 (MQ=255)
                         aCTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGACAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACAtt                                                                                                                          <  2:503858/150‑1 (MQ=255)
                          ctctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTc                                                                                                                         <  1:4175928/150‑1 (MQ=255)
                          ctctGCGTGGTCAGGGACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTc                                                                                                                         <  2:3907940/150‑1 (MQ=255)
                            ctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCAt                                                                                                                       <  1:1859467/150‑1 (MQ=255)
                            ctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCAt                                                                                                                       <  2:2737219/150‑1 (MQ=255)
                            ctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCAt                                                                                                                       <  1:1248247/150‑1 (MQ=255)
                            ctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCAt                                                                                                                       >  1:484571/1‑150 (MQ=255)
                                gTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCtg                                                                                                                   <  1:1920108/150‑1 (MQ=255)
                                 tGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCtgt                                                                                                                  <  1:4203804/150‑1 (MQ=255)
                                   gTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGa                                                                                                                >  1:354476/1‑150 (MQ=255)
                                    tCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGaa                                                                                                               >  2:2912022/1‑150 (MQ=255)
                                     cAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAg                                                                                                              >  2:3473244/1‑150 (MQ=255)
                                     cAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAg                                                                                                              >  1:3307290/1‑150 (MQ=255)
                                      aGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAgc                                                                                                             <  1:1036286/150‑1 (MQ=255)
                                      aGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAgc                                                                                                             >  2:3271204/1‑150 (MQ=255)
                                       ggTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAgcg                                                                                                            >  2:137963/1‑150 (MQ=255)
                                        gTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAgcgc                                                                                                           <  2:1451438/150‑1 (MQ=255)
                                        gTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGAGTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAgcgc                                                                                                           <  2:2488770/150‑1 (MQ=255)
                                         tACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCa                                                                                                          >  1:3454059/1‑150 (MQ=255)
                                          aCTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCAt                                                                                                         <  2:4713882/150‑1 (MQ=255)
                                              ccGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAg                                                                                                     <  2:830387/150‑1 (MQ=255)
                                                 aTCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGaa                                                                                                  >  1:1392082/1‑150 (MQ=255)
                                                   cGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGt                                                                                                >  1:4475671/1‑150 (MQ=255)
                                                    ggtggtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGtt                                                                                               >  2:4070572/1‑150 (MQ=255)
                                                      tggtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCt                                                                                             >  2:2233958/1‑150 (MQ=255)
                                                      tggtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCt                                                                                             >  1:285825/1‑150 (MQ=255)
                                                            gaaagaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGActct                                                                                       <  1:1189280/150‑1 (MQ=255)
                                                             aaagaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTg                                                                                      <  2:1443383/150‑1 (MQ=255)
                                                             aaagaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTg                                                                                      >  2:2081920/1‑150 (MQ=255)
                                                               agaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGgc                                                                                    >  2:1861829/1‑150 (MQ=255)
                                                               agaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGgc                                                                                    <  2:2716103/150‑1 (MQ=255)
                                                                     aCTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCa                                                                              >  1:3915733/1‑150 (MQ=255)
                                                                     aCTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCa                                                                              >  1:2665884/1‑150 (MQ=255)
                                                                      cTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCaa                                                                             >  2:1095899/1‑150 (MQ=255)
                                                                      cTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCaa                                                                             >  2:4443116/1‑150 (MQ=255)
                                                                       tCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAACTAAACGCGATCTGGACGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCaaa                                                                            >  2:4480711/1‑150 (MQ=255)
                                                                          tGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATg                                                                         >  2:3768773/1‑150 (MQ=255)
                                                                          tGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATg                                                                         >  1:3988898/1‑150 (MQ=255)
                                                                              aGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAAcc                                                                     <  1:2524289/150‑1 (MQ=255)
                                                                              aGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAAcc                                                                     <  1:3474705/150‑1 (MQ=255)
                                                                              aGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGAGTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAAcc                                                                     <  2:4173889/150‑1 (MQ=255)
                                                                               gcgcATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCg                                                                    <  2:240306/150‑1 (MQ=255)
                                                                                 gcATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTc                                                                  >  1:963765/1‑150 (MQ=255)
                                                                                 gcATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTc                                                                  >  2:1138060/1‑150 (MQ=255)
                                                                                   aTCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTg                                                                >  1:3117282/1‑150 (MQ=255)
                                                                                   aTCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTg                                                                <  1:2869103/150‑1 (MQ=255)
                                                                                   aTCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTg                                                                <  1:1265852/150‑1 (MQ=255)
                                                                                       cggcggGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTAt                                                            <  1:2907269/150‑1 (MQ=255)
                                                                                        ggcggGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATg                                                           >  2:754829/1‑150 (MQ=255)
                                                                                         gcggTATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGc                                                          <  1:4632413/150‑1 (MQ=255)
                                                                                            ggATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAgg                                                       <  2:1212199/150‑1 (MQ=255)
                                                                                            ggATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAgg                                                       <  2:1211770/150‑1 (MQ=255)
                                                                                              tttAAAACAGTGCTTATTCCGTTCGAAAATAATCGCGATCTGGACGAGAGTCCTTACAACGTAAGTGCCGATCTGGACATTCATCCTTTGAAGCTCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGtt                                                     <  2:828667/149‑1 (MQ=255)
                                                                                              aTTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGtt                                                     <  2:2758673/150‑1 (MQ=255)
                                                                                               ttAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTtg                                                    >  2:540666/1‑150 (MQ=255)
                                                                                               ttAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTtg                                                    >  2:4638322/1‑150 (MQ=255)
                                                                                                tAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTtgt                                                   <  2:700472/150‑1 (MQ=255)
                                                                                                tAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTtgt                                                   >  2:1814087/1‑150 (MQ=255)
                                                                                                 aaaaCAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTtgtg                                                  >  1:596471/1‑150 (MQ=255)
                                                                                                 aaaaCAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTtgtg                                                  >  1:172208/1‑150 (MQ=255)
                                                                                                  aaaCAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGa                                                 >  1:658679/1‑150 (MQ=255)
                                                                                                  aaaCAGTGCTAATTCCGTTCGAAAATAAACGAGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGa                                                 <  1:4008289/150‑1 (MQ=255)
                                                                                                            aaTTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATa                                       <  1:1036047/150‑1 (MQ=255)
                                                                                                                ccGTTAGAAAATAAACGAGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTTGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGa                                   <  2:2921414/150‑1 (MQ=255)
                                                                                                                  gTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGAcc                                 >  1:4383684/1‑150 (MQ=255)
                                                                                                                     cGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTcg                              <  1:4379024/150‑1 (MQ=255)
                                                                                                                     cGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTcg                              >  1:506074/1‑150 (MQ=255)
                                                                                                                        aaaTAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCgcgc                           <  2:1057114/150‑1 (MQ=255)
                                                                                                                         aaTAAACTCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCa                          <  2:1198614/150‑1 (MQ=255)
                                                                                                                          aTAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCaa                         >  2:2606131/1‑150 (MQ=255)
                                                                                                                            aaaCGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCaaaa                       >  2:867707/1‑150 (MQ=255)
                                                                                                                             aaCGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAAt                      <  1:1646913/150‑1 (MQ=255)
                                                                                                                              aCGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATg                     <  1:3784633/150‑1 (MQ=255)
                                                                                                                                gcgATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCa                   >  1:3395472/1‑150 (MQ=255)
                                                                                                                                  gATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACt                 <  2:3700343/150‑1 (MQ=255)
                                                                                                                                  gATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACt                 >  2:4393166/1‑150 (MQ=255)
                                                                                                                                     cTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACtaat              >  2:3320828/1‑150 (MQ=255)
                                                                                                                                     cTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACtaat              >  1:3221945/1‑150 (MQ=255)
                                                                                                                                       ggAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACtaataa            >  1:2938793/1‑150 (MQ=255)
                                                                                                                                          agagATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATaaaaa         <  2:3819450/150‑1 (MQ=255)
                                                                                                                                           gagaTTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAAc        >  2:2658851/1‑150 (MQ=255)
                                                                                                                                           gagaTTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAAc        <  1:1469002/150‑1 (MQ=255)
                                                                                                                                            agaTTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACa       >  2:454086/1‑150 (MQ=255)
                                                                                                                                              aTTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAgg     >  2:2064161/1‑150 (MQ=255)
                                                                                                                                              aTTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAgg     <  1:1436854/150‑1 (MQ=255)
                                                                                                                                               ttCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAggg    <  2:4611081/150‑1 (MQ=255)
                                                                                                                                                tCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGc   >  1:2266921/1‑150 (MQ=255)
                                                                                                                                                tCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGc   >  2:1807/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCT  >  CP009273/456440‑456734

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: