Predicted mutation
evidence seq id position mutation annotation gene description
RA AP012306 1,438,157 G→A A293T (GCA→ACA)  pykF → pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,438,1570GA93.7% 1299.9 / 53.7 367A293T (GCA→ACA) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (16/7);  new base A (215/129);  total (231/136)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

Only 100 of 379 total aligned reads displayed.

GTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGC  >  AP012306/1438011‑1438303
                                                                                                                                                  |                                                                                                                                                  
 tGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGa                                                                                                                                                <  1:4142668/150‑1 (MQ=255)
    cGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGc                                                                                                                                             <  2:3354118/150‑1 (MQ=255)
        cTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGt                                                                                                                                         <  2:3446876/150‑1 (MQ=255)
         tGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTg                                                                                                                                        <  2:3793813/150‑1 (MQ=255)
         tGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTg                                                                                                                                        >  2:1659565/1‑150 (MQ=255)
         tGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATCCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTg                                                                                                                                        <  2:3894021/150‑1 (MQ=255)
           ggTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGAc                                                                                                                                      >  2:3733254/1‑150 (MQ=255)
            gtgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACg                                                                                                                                     <  2:2839719/150‑1 (MQ=255)
             tgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGt                                                                                                                                    >  1:3307012/1‑150 (MQ=255)
              gtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGtt                                                                                                                                   <  1:951063/150‑1 (MQ=255)
                 gAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCa                                                                                                                                >  1:1016951/1‑150 (MQ=255)
                  aaaTCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCaa                                                                                                                               >  1:3627801/1‑150 (MQ=255)
                   aaTCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCaaa                                                                                                                              >  1:3214314/1‑150 (MQ=255)
                    aTCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAAc                                                                                                                             >  2:2932673/1‑150 (MQ=255)
                     tCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACg                                                                                                                            >  1:3501750/1‑150 (MQ=255)
                         ggTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCAt                                                                                                                        <  1:624400/150‑1 (MQ=255)
                           tAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCAATcc                                                                                                                      >  2:2648481/1‑150 (MQ=255)
                            agaagaAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCt                                                                                                                     <  2:679324/150‑1 (MQ=255)
                             gaagaaGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTc                                                                                                                    >  2:523931/1‑150 (MQ=255)
                               agaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGa                                                                                                                  >  1:2104050/1‑150 (MQ=255)
                                  agTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACgg                                                                                                               >  1:1086848/1‑150 (MQ=255)
                                  agTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACgg                                                                                                               >  2:3995208/1‑150 (MQ=255)
                                   gTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGt                                                                                                              >  2:1418737/1‑150 (MQ=255)
                                      aTCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACt                                                                                                           >  1:2440739/1‑150 (MQ=255)
                                        cTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGa                                                                                                         >  1:564819/1‑150 (MQ=255)
                                           cGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGc                                                                                                      <  1:1729504/150‑1 (MQ=255)
                                               cAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTg                                                                                                  <  2:3642377/150‑1 (MQ=255)
                                                agaagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGa                                                                                                 >  1:2901752/1‑150 (MQ=255)
                                                  aagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATg                                                                                               <  2:2429845/150‑1 (MQ=255)
                                                  aagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATg                                                                                               >  2:878303/1‑150 (MQ=255)
                                                  aagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATg                                                                                               <  1:2674106/150‑1 (MQ=255)
                                                    gatgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCt                                                                                             <  1:358346/150‑1 (MQ=255)
                                                      tgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGt                                                                                           <  2:4261274/150‑1 (MQ=255)
                                                       gatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCTCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTc                                                                                          >  1:4051372/1‑150 (MQ=255)
                                                        atgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCt                                                                                         >  1:3458797/1‑150 (MQ=255)
                                                        atgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCt                                                                                         <  2:82113/150‑1 (MQ=255)
                                                         tgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTg                                                                                        >  2:2195999/1‑150 (MQ=255)
                                                         tgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTg                                                                                        <  1:749526/150‑1 (MQ=255)
                                                         tgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTg                                                                                        <  2:1833269/150‑1 (MQ=255)
                                                          gatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTgg                                                                                       <  2:1080577/150‑1 (MQ=255)
                                                           atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGt                                                                                      <  2:1285898/150‑1 (MQ=255)
                                                             cGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGa                                                                                    >  1:1726281/1‑150 (MQ=255)
                                                              gAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGaa                                                                                   >  1:1695627/1‑150 (MQ=255)
                                                                  aaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCg                                                                               >  2:1723402/1‑150 (MQ=255)
                                                                  aaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCg                                                                               >  2:126318/1‑150 (MQ=255)
                                                                     gTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaa                                                                            >  2:329348/1‑150 (MQ=255)
                                                                     gTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaa                                                                            >  1:2305241/1‑150 (MQ=255)
                                                                     gTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaa                                                                            >  1:968642/1‑150 (MQ=255)
                                                                      tATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaaa                                                                           >  1:385658/1‑150 (MQ=255)
                                                                       aTCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaaaa                                                                          <  2:19487/150‑1 (MQ=255)
                                                                       aTCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaaaa                                                                          >  2:1561765/1‑150 (MQ=255)
                                                                       aTCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaaaa                                                                          <  2:926361/150‑1 (MQ=255)
                                                                        tCCGTGCACGTAAAGTCGTTACCACTGCGACCAAGATGCTGGATTCCATGATCAAAACCCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAg                                                                         <  2:2049955/150‑1 (MQ=255)
                                                                          cGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGt                                                                       >  2:1188215/1‑150 (MQ=255)
                                                                           gTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTa                                                                      >  1:2585506/1‑150 (MQ=255)
                                                                            tGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTaa                                                                     >  1:671096/1‑150 (MQ=255)
                                                                             gCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTaaa                                                                    >  2:354200/1‑150 (MQ=255)
                                                                              cACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAAt                                                                   >  2:1359343/1‑150 (MQ=255)
                                                                              cACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAAt                                                                   <  2:3070616/150‑1 (MQ=255)
                                                                                cGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAc                                                                 >  1:2202362/1‑150 (MQ=255)
                                                                                cGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAc                                                                 >  2:3756650/1‑150 (MQ=255)
                                                                                cGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAc                                                                 >  2:415210/1‑150 (MQ=255)
                                                                                 gTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAcc                                                                >  1:447102/1‑150 (MQ=255)
                                                                                 gTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAcc                                                                >  1:186517/1‑150 (MQ=255)
                                                                                    aaGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGc                                                             <  1:1249903/150‑1 (MQ=255)
                                                                                    aaGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGc                                                             >  1:2747238/1‑150 (MQ=255)
                                                                                        cGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGa                                                         <  1:3671418/150‑1 (MQ=255)
                                                                                          ttATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAg                                                       >  2:4058924/1‑150 (MQ=255)
                                                                                           tATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGc                                                      <  2:1485615/150‑1 (MQ=255)
                                                                                           tATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGc                                                      >  1:2975889/1‑150 (MQ=255)
                                                                                            aTCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCg                                                     >  2:2124041/1‑150 (MQ=255)
                                                                                             tCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCgg                                                    >  1:854997/1‑150 (MQ=255)
                                                                                              cACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGt                                                   <  2:2279900/150‑1 (MQ=255)
                                                                                                    gACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTAt                                             <  1:1001088/150‑1 (MQ=255)
                                                                                                    gACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTAt                                             >  2:1755449/1‑150 (MQ=255)
                                                                                                      cccAGATGCTGGATTCCATGATCAAAACCCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCa                                           <  2:209002/150‑1 (MQ=255)
                                                                                                      cccAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCa                                           >  1:3790177/1‑150 (MQ=255)
                                                                                                       ccAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCAt                                          >  2:610284/1‑150 (MQ=255)
                                                                                                          gATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGc                                       <  1:3155541/150‑1 (MQ=255)
                                                                                                           aTGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCg                                      >  1:2634917/1‑150 (MQ=255)
                                                                                                                 gATTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATc                                <  1:1139490/150‑1 (MQ=255)
                                                                                                                  aTTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCt                               >  2:1799573/1‑150 (MQ=255)
                                                                                                                  aTTCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCCCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCt                               <  1:255228/150‑1 (MQ=255)
                                                                                                                   ttCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTg                              >  2:2627425/1‑150 (MQ=255)
                                                                                                                   ttCCATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTg                              >  2:2347545/1‑150 (MQ=255)
                                                                                                                     ccATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCg                            >  2:2979029/1‑150 (MQ=255)
                                                                                                                     ccATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCg                            >  1:2830471/1‑150 (MQ=255)
                                                                                                                      cATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGa                           <  1:2941350/150‑1 (MQ=255)
                                                                                                                      cATGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGa                           >  1:1503744/1‑150 (MQ=255)
                                                                                                                        tGATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAAc                         <  2:3180841/150‑1 (MQ=255)
                                                                                                                         gATCAAAAACCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACg                        >  1:1347234/1‑150 (MQ=255)
                                                                                                                             aaaaaCCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTacc                    >  2:3395711/1‑150 (MQ=255)
                                                                                                                                aaCCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTaccgac                 >  1:420346/1‑150 (MQ=255)
                                                                                                                                 aCCCACGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTaccgacc                >  1:1591669/1‑150 (MQ=255)
                                                                                                                                     aCGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGt            <  1:1955549/150‑1 (MQ=255)
                                                                                                                                     aCGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGt            >  2:2527989/1‑150 (MQ=255)
                                                                                                                                      cGCCCGACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATGCCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGtg           <  2:3950507/150‑1 (MQ=255)
                                                                                                                                           gACTCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGaa      >  2:1631541/1‑150 (MQ=255)
                                                                                                                                              tCGCACAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAg   >  1:374737/1‑150 (MQ=255)
                                                                                                                                               cgcaCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGc  >  1:3018181/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGC  >  AP012306/1438011‑1438303

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