Predicted mutation
evidence seq id position mutation annotation gene description
RA AP012306 1,438,241 C→A P321T (CCG→ACG)  pykF → pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,438,2410CA99.4% 1925.3 / ‑0.9 477P321T (CCG→ACG) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (274/200);  total (275/202)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 482 total aligned reads displayed.

AAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTG  >  AP012306/1438094‑1438387
                                                                                                                                                   |                                                                                                                                                  
 aaGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGc                                                                                                                                                 >  1:3639494/1‑150 (MQ=255)
  aGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCt                                                                                                                                                >  1:4715897/1‑150 (MQ=255)
   gTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTg                                                                                                                                               >  1:2624807/1‑150 (MQ=255)
     cGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGa                                                                                                                                             >  2:2920920/1‑150 (MQ=255)
     cGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGa                                                                                                                                             >  1:1241425/1‑150 (MQ=255)
      gTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGaa                                                                                                                                            >  2:5036595/1‑150 (MQ=255)
      gTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGaa                                                                                                                                            >  1:5253744/1‑150 (MQ=255)
        tATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGc                                                                                                                                          >  2:2724052/1‑150 (MQ=255)
          tCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCgg                                                                                                                                        >  1:5084285/1‑150 (MQ=255)
            aCTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGtt                                                                                                                                      <  1:746954/150‑1 (MQ=255)
             cTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGttt                                                                                                                                     <  1:4870959/150‑1 (MQ=255)
                 gACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTAt                                                                                                                                 <  2:3329487/150‑1 (MQ=255)
                 gACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTAt                                                                                                                                 >  2:5249245/1‑150 (MQ=255)
                   cccAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCa                                                                                                                               >  1:4621627/1‑150 (MQ=255)
                     cAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATg                                                                                                                             >  2:124570/1‑150 (MQ=255)
                     cAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATg                                                                                                                             >  1:2351359/1‑150 (MQ=255)
                     cAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATg                                                                                                                             >  2:3499353/1‑150 (MQ=255)
                      aGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATgg                                                                                                                            <  1:629377/150‑1 (MQ=255)
                      aGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATgg                                                                                                                            <  2:67419/150‑1 (MQ=255)
                      aGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATgg                                                                                                                            >  1:2987932/1‑150 (MQ=255)
                          gCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGAc                                                                                                                        >  1:3585389/1‑150 (MQ=255)
                              gATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATc                                                                                                                    >  1:1000725/1‑150 (MQ=255)
                               aTTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCt                                                                                                                   >  1:3829510/1‑150 (MQ=255)
                                ttCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGAGCATCTg                                                                                                                  >  1:4760305/1‑150 (MQ=255)
                                    aTGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGaa                                                                                                              >  1:2533238/1‑150 (MQ=255)
                                     tGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAAc                                                                                                             >  1:879800/1‑150 (MQ=255)
                                     tGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAAc                                                                                                             <  1:1486215/150‑1 (MQ=255)
                                       aTCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGt                                                                                                           >  1:240901/1‑150 (MQ=255)
                                        tCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTa                                                                                                          <  1:1612148/150‑1 (MQ=255)
                                           aaaaCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTaccg                                                                                                       <  2:3817600/150‑1 (MQ=255)
                                           aaaaCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTaccg                                                                                                       <  2:4223510/150‑1 (MQ=255)
                                            aaaCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTaccga                                                                                                      >  1:1071107/1‑150 (MQ=255)
                                             aaCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTaccgac                                                                                                     >  2:3997452/1‑150 (MQ=255)
                                                ccACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACcgc                                                                                                  >  1:4563512/1‑150 (MQ=255)
                                                    gCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGtga                                                                                              >  1:3418993/1‑150 (MQ=255)
                                                      ccGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACcgcagatgg                                                                                            >  1:137317/1‑144 (MQ=255)
                                                       cGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGtgatga                                                                                           >  2:4099167/1‑150 (MQ=255)
                                                        gACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGaa                                                                                          >  2:1672148/1‑150 (MQ=255)
                                                          cTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACa                                                                                        >  1:35490/1‑150 (MQ=255)
                                                          cTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACa                                                                                        >  2:3183607/1‑150 (MQ=255)
                                                           tCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGc                                                                                      >  1:2822475/1‑150 (MQ=255)
                                                             gcgcAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGcc                                                                                     <  1:3693513/150‑1 (MQ=255)
                                                               gcAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGt                                                                                   >  2:4707384/1‑150 (MQ=255)
                                                                cAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGtc                                                                                  >  2:2571942/1‑150 (MQ=255)
                                                                 aGAAGCCGGTGACTTTGCAAACGCCATCCTCGACGGTCCTGACGCAGTGATGCTGTCTGGTGATTCCGCAAAAGGTAAATACACGCTGGAAGCAGTTTCTATCATGCCGACCAGCTGCAAACGTACTGACCGCGTGATGAACAGCCGtct                                                                                 <  2:2501119/150‑1 (MQ=255)
                                                                  gAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGtctc                                                                                <  2:2177719/150‑1 (MQ=255)
                                                                   aaGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCg                                                                               <  2:2086391/150‑1 (MQ=255)
                                                                     gCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAg                                                                             <  2:4968138/150‑1 (MQ=255)
                                                                      ccGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGt                                                                            <  2:2389858/150‑1 (MQ=255)
                                                                      ccGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGt                                                                            <  2:4648162/150‑1 (MQ=255)
                                                                       cGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGtt                                                                           <  2:4442716/150‑1 (MQ=255)
                                                                          tGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTcaa                                                                        <  1:257773/150‑1 (MQ=255)
                                                                          tGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTcaa                                                                        <  2:1922184/150‑1 (MQ=255)
                                                                          tGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTcaa                                                                        <  1:132567/150‑1 (MQ=255)
                                                                             cGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTcaacaa                                                                     <  2:4828509/150‑1 (MQ=255)
                                                                              gTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAAt                                                                    <  1:4977207/150‑1 (MQ=255)
                                                                                tGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGa                                                                  >  2:2292912/1‑150 (MQ=255)
                                                                                 gCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGAc                                                                 <  1:1871252/150‑1 (MQ=255)
                                                                                  cAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACa                                                                <  1:2053596/150‑1 (MQ=255)
                                                                                    aaCGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAAc                                                              <  1:5165743/150‑1 (MQ=255)
                                                                                        ccATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTa                                                          >  1:2869772/1‑150 (MQ=255)
                                                                                            ccTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACt                                                      >  2:1536805/1‑150 (MQ=255)
                                                                                             cTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTg                                                     <  2:3918477/150‑1 (MQ=255)
                                                                                              tCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTgc                                                    <  1:212400/150‑1 (MQ=255)
                                                                                                  cGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCAt                                                <  2:4493321/150‑1 (MQ=255)
                                                                                                   ggTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCAtt                                               <  1:1627502/150‑1 (MQ=255)
                                                                                                    gTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTa                                              <  2:1378012/150‑1 (MQ=255)
                                                                                                     tACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTAc                                             <  2:3447994/150‑1 (MQ=255)
                                                                                                     tACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTAc                                             >  2:3086947/1‑150 (MQ=255)
                                                                                                        tGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGGTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGa                                          >  2:5117261/1‑150 (MQ=255)
                                                                                                          aCGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAg                                        <  1:271272/150‑1 (MQ=255)
                                                                                                            gCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCg                                      <  2:1442169/150‑1 (MQ=255)
                                                                                                            gCAGGGATGCGGTCTGTTCAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACGATCTGCGAGCGTACGGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCg                                      <  2:3718395/150‑1 (MQ=255)
                                                                                                              aGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGt                                    >  1:102357/1‑150 (MQ=255)
                                                                                                                 gATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACACTGACAACCGTAAACTGCGCATTACCGAAGCGGTATg                                 >  1:4248629/1‑150 (MQ=255)
                                                                                                                 gATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATg                                 <  2:876413/150‑1 (MQ=255)
                                                                                                                 gATGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATg                                 >  2:77193/1‑150 (MQ=255)
                                                                                                                   tGCTGTCTGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGcc                               >  1:5147972/1‑150 (MQ=255)
                                                                                                                          tGGTGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGc                        <  2:4548134/150‑1 (MQ=255)
                                                                                                                             tGAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGt                     <  1:3964395/150‑1 (MQ=255)
                                                                                                                              gAATCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGtt                    <  2:1333697/150‑1 (MQ=255)
                                                                                                                               aaTCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTg                   >  1:4104780/1‑150 (MQ=255)
                                                                                                                                aTCCGCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGa                  <  2:4899695/150‑1 (MQ=255)
                                                                                                                                    gCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAAct              <  2:119054/150‑1 (MQ=255)
                                                                                                                                    gCAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAAct              <  1:551981/150‑1 (MQ=255)
                                                                                                                                     cAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAActg             >  1:3256860/1‑150 (MQ=255)
                                                                                                                                     cAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAActg             <  1:3183238/150‑1 (MQ=255)
                                                                                                                                     cAAAAGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAActg             >  1:308369/1‑150 (MQ=255)
                                                                                                                                      aaaaGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAActgc            >  1:4389184/1‑150 (MQ=255)
                                                                                                                                      aaaaGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAActgc            >  1:4473165/1‑150 (MQ=255)
                                                                                                                                        aaGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAActgctg          >  1:177551/1‑150 (MQ=255)
                                                                                                                                         aGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGa         >  2:5197257/1‑150 (MQ=255)
                                                                                                                                         aGGTAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGa         <  1:1407190/150‑1 (MQ=255)
                                                                                                                                            tAAATACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGaaaa      <  2:178134/150‑1 (MQ=255)
                                                                                                                                             aaaTACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGaaaaa     >  2:3255105/1‑150 (MQ=255)
                                                                                                                                             aaaTACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTCGAAACTGCTGaaaaa     >  2:4328777/1‑150 (MQ=255)
                                                                                                                                              aaTACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAAc    <  2:3569989/150‑1 (MQ=255)
                                                                                                                                               aTACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACt   <  2:1018579/150‑1 (MQ=255)
                                                                                                                                               aTACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACt   >  1:1934012/1‑150 (MQ=255)
                                                                                                                                                tACACGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTg  <  2:27797/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTG  >  AP012306/1438094‑1438387

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: