Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,527,739 C→A 6.8% F521L (TTC→TTA)  pck → phosphoenolpyruvate carboxykinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,527,7390CA 6.8% 769.8 / 30.1 222F521L (TTC→TTA) pckphosphoenolpyruvate carboxykinase
Reads supporting (aligned to +/- strand):  ref base C (109/98);  new base A (8/7);  total (117/105)
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.46e-01

Only 100 of 224 total aligned reads displayed.

CCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTACCCTTCTTTTG  >  CP009273/3527591‑3527888
                                                                                                                                                    |                                                                                                                                                     
ccTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACACCGTCg                                                                                                                                                      >  2:4840201/1‑150 (MQ=255)
ccTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCg                                                                                                                                                      >  2:952699/1‑150 (MQ=255)
 cTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGa                                                                                                                                                     >  1:1680578/1‑150 (MQ=255)
  ttCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGAt                                                                                                                                                    <  2:1102679/150‑1 (MQ=255)
   tCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATa                                                                                                                                                   <  1:4186103/150‑1 (MQ=255)
      ctctGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAAt                                                                                                                                                >  1:2989669/1‑150 (MQ=255)
        ctGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAACCTGTTTATCGACAACTTCGATAAATac                                                                                                                                              >  2:4880530/1‑150 (MQ=255)
         tGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATaca                                                                                                                                             <  2:1700494/150‑1 (MQ=255)
          gCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATacac                                                                                                                                            <  1:690391/150‑1 (MQ=255)
              aTGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGac                                                                                                                                        <  1:4665507/150‑1 (MQ=255)
                gTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGacac                                                                                                                                      >  2:153057/1‑150 (MQ=255)
                gTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGacac                                                                                                                                      >  1:4390952/1‑150 (MQ=255)
                 tttAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACAcc                                                                                                                                     >  2:3947770/1‑150 (MQ=255)
                 tttAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACAcc                                                                                                                                     >  1:4457169/1‑150 (MQ=255)
                   tAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACAcccc                                                                                                                                   <  2:2829233/150‑1 (MQ=255)
                    aaCCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCt                                                                                                                                  >  1:2335477/1‑150 (MQ=255)
                    aaCCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCt                                                                                                                                  >  2:2243466/1‑150 (MQ=255)
                     aCCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTAGATAAATACACCGACACCCCTg                                                                                                                                 >  2:92289/1‑150 (MQ=255)
                      ccTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGc                                                                                                                                <  1:3093257/150‑1 (MQ=255)
                       cTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCg                                                                                                                               >  1:4282600/1‑150 (MQ=255)
                          gCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCGGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGt                                                                                                                            <  2:1715451/150‑1 (MQ=255)
                            gATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGc                                                                                                                          <  1:630449/150‑1 (MQ=255)
                             aTCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGcc                                                                                                                         <  1:4048528/150‑1 (MQ=255)
                               cccAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCcgc                                                                                                                       >  1:4391073/1‑150 (MQ=255)
                                  aaCCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCt                                                                                                                    <  2:4697638/150‑1 (MQ=255)
                                    ccGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTgg                                                                                                                  >  2:4964793/1‑150 (MQ=255)
                                       aaCTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGAAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAg                                                                                                               >  1:3053971/1‑150 (MQ=255)
                                        aCTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGc                                                                                                              <  2:3564690/150‑1 (MQ=255)
                                         cTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCg                                                                                                             <  1:3784525/150‑1 (MQ=255)
                                         cTGCCGCGCGAAGACACGAATATTCACGATCCGCGTAACACCTAAGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTTGCGAAACTGTTTAACGTCCCCTTCGATAAAAACACCGACACCCCTGCGGGTGACGCGCTGGTAGCg                                                                                                             >  2:125832/1‑150 (MQ=255)
                                            ccGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTAGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCt                                                                                                          <  1:1020944/150‑1 (MQ=255)
                                              gggCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTgg                                                                                                        >  1:332547/1‑150 (MQ=255)
                                               ggCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGt                                                                                                       <  1:2437726/150‑1 (MQ=255)
                                                 cGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTcc                                                                                                     >  2:4392502/1‑150 (MQ=255)
                                                  gTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCg                                                                                                    >  1:32157/1‑150 (MQ=255)
                                                    aGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGaa                                                                                                  >  1:4124508/1‑150 (MQ=255)
                                                    aGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGaa                                                                                                  >  1:1667537/1‑150 (MQ=255)
                                                     gACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGaaa                                                                                                 <  1:3930181/150‑1 (MQ=255)
                                                     gACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGaaa                                                                                                 <  1:1411173/150‑1 (MQ=255)
                                                      acacGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAAc                                                                                                <  2:3092581/150‑1 (MQ=255)
                                                       cacGAAGAGTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACt                                                                                               <  2:17387/150‑1 (MQ=255)
                                                        acGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTg                                                                                              <  1:1216859/150‑1 (MQ=255)
                                                           aaGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTaa                                                                                           >  2:2719387/1‑150 (MQ=255)
                                                              aTTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTCGTAGCGGCTGGTCCGAAACTGTAATGa                                                                                        >  1:929194/1‑150 (MQ=255)
                                                                   cGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGa                                                                                   >  2:2609921/1‑150 (MQ=255)
                                                                   cGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTAGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGa                                                                                   <  2:4680328/150‑1 (MQ=255)
                                                                       ccGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCt                                                                               <  1:3854664/150‑1 (MQ=255)
                                                                        cgcgTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTg                                                                              >  1:386637/1‑150 (MQ=255)
                                                                         gcgTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTgg                                                                             <  2:2072981/150‑1 (MQ=255)
                                                                          cgTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCt                                                                               >  1:4912231/1‑147 (MQ=255)
                                                                          cgTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGga                                                                            <  1:2675482/150‑1 (MQ=255)
                                                                          cgTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTAGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGga                                                                            >  1:1279914/1‑150 (MQ=255)
                                                                            tAACACCGACGCTTCTCCGGTACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGgaga                                                                          <  2:660955/150‑1 (MQ=255)
                                                                             aaCACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGaa                                                                         >  1:1205905/1‑150 (MQ=255)
                                                                              acacCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAat                                                                        <  1:952699/150‑1 (MQ=255)
                                                                               cacCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAata                                                                       <  2:4482298/150‑1 (MQ=255)
                                                                                   tACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTCATGATTTGAAGCTGGAGAAtatcta                                                                   >  1:742083/1‑150 (MQ=255)
                                                                                      gCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATcc                                                                >  2:1112986/1‑150 (MQ=255)
                                                                                      gCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATcc                                                                >  2:2222354/1‑150 (MQ=255)
                                                                                          ctcCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTa                                                            >  1:4779127/1‑150 (MQ=255)
                                                                                          ctcCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCAGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTa                                                            >  1:2154713/1‑150 (MQ=255)
                                                                                           tcCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTAGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTAt                                                           >  1:5057673/1‑150 (MQ=255)
                                                                                              ggAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCtt                                                        <  2:4620595/150‑1 (MQ=255)
                                                                                              ggAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCAAGTATCtt                                                        <  1:3756210/150‑1 (MQ=255)
                                                                                               gAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTCTCTta                                                       >  2:4935508/1‑150 (MQ=255)
                                                                                                 aCAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTtata                                                     <  2:333620/150‑1 (MQ=255)
                                                                                                     tGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGaaa                                                 <  1:4880530/150‑1 (MQ=255)
                                                                                                     tGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGaaa                                                 >  2:717898/1‑150 (MQ=255)
                                                                                                        cAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCa                                              <  1:2929609/150‑1 (MQ=255)
                                                                                                           gAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCaaaa                                           <  1:58458/150‑1 (MQ=255)
                                                                                                            aaaaaGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAAc                                          <  1:3704492/150‑1 (MQ=255)
                                                                                                             aaaaGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGACAGCAAAACg                                         >  2:1600870/1‑150 (MQ=255)
                                                                                                             aaaaGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTAGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACg                                         <  2:3090008/150‑1 (MQ=255)
                                                                                                             aaaaGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTAGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACg                                         <  2:617042/150‑1 (MQ=255)
                                                                                                              aaaGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACgg                                        >  1:2597502/1‑150 (MQ=255)
                                                                                                               aaGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACggg                                       >  1:3688355/1‑150 (MQ=255)
                                                                                                                 gCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAg                                     <  2:815483/150‑1 (MQ=255)
                                                                                                                        cccTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCtt                              <  1:2490274/150‑1 (MQ=255)
                                                                                                                         ccTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGGGCTGGTTGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCTCCTTc                             <  1:4552013/150‑1 (MQ=255)
                                                                                                                         ccTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTc                             <  2:2883708/150‑1 (MQ=255)
                                                                                                                         ccTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTc                             >  1:2297291/1‑150 (MQ=255)
                                                                                                                          cTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCg                            >  2:4277290/1‑150 (MQ=255)
                                                                                                                                 aaCTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCg                     <  1:512288/150‑1 (MQ=255)
                                                                                                                                 aaCTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCg                     >  1:2548301/1‑150 (MQ=255)
                                                                                                                                  aCTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGt                    <  2:864867/150‑1 (MQ=255)
                                                                                                                                   cTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGtt                   <  2:3483926/150‑1 (MQ=255)
                                                                                                                                      tttATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGtttat                >  2:1938178/1‑150 (MQ=255)
                                                                                                                                       ttATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGtttatt               >  1:279708/1‑150 (MQ=255)
                                                                                                                                         aTCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGtttattta             <  2:3648337/150‑1 (MQ=255)
                                                                                                                                         aTCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTAGCTTATAGAAAGCAAAACGGGAGGCACCTGCGCCTCCCGttttgtta             >  1:772029/1‑145 (MQ=255)
                                                                                                                                          tCGACAACTTCGATAAATACACCGAGACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTAc            <  2:3239211/150‑1 (MQ=255)
                                                                                                                                          tCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTAc            <  1:987245/150‑1 (MQ=255)
                                                                                                                                            gACAACTTAGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGGTTATTTAccc          >  2:3655796/1‑150 (MQ=255)
                                                                                                                                              cAACTTCGATAAATACACCGACACCCCTGCGGGTGTCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTACCctt        >  1:3482392/1‑150 (MQ=255)
                                                                                                                                              cAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTACCctt        >  2:3046339/1‑150 (MQ=255)
                                                                                                                                              cAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTACCctt        <  2:237401/150‑1 (MQ=255)
                                                                                                                                              cAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTACCctt        >  1:840454/1‑150 (MQ=255)
                                                                                                                                                aCTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTACCcttct      >  2:691562/1‑150 (MQ=255)
                                                                                                                                                   tCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTACCCTTCtttt   <  1:1954488/150‑1 (MQ=255)
                                                                                                                                                    cGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTGTATTTACCCTTCTTTTg  >  2:4579666/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTACCCTTCTTTTG  >  CP009273/3527591‑3527888

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: ‑