Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP012306800,6750AC26.2% 617.3 / 29.6 263T192P (ACC→CCC) lolAchaperone for lipoproteins
Reads supporting (aligned to +/- strand):  ref base A (145/49);  new base C (0/69);  total (145/118)
Fisher's exact test for biased strand distribution p-value = 1.49e-31
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 322 total aligned reads displayed.

GCAATGGCAATCTGAAGCAGTTCACCATTAACGTGGGACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCGGCCAGCAACA  >  AP012306/800526‑800824
                                                                                                                                                     |                                                                                                                                                     
      gCAATCTGAAGCAGTTCACCATTAACGTGGGACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCccgc                                                                                                                                                 >  2:1060948/1‑150 (MQ=255)
       cAATCTGAAGCAGTTCACCATTAACGTGGGACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCccgcc                                                                                                                                                <  2:544092/150‑1 (MQ=255)
        aaTCTGAAGCAGTTCACCATTAACGTGGGACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCccgccg                                                                                                                                               >  1:3112731/1‑150 (MQ=255)
             gAAGCAGTTCACCATTAACGTGGGACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAgg                                                                                                                                          >  1:2212074/1‑150 (MQ=255)
              aaGCAGTTCACCATTAACGTGGGACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGc                                                                                                                                         >  2:3765656/1‑150 (MQ=255)
                     tccccATTAACGGGGGACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCATTTATCAACGGAAACCCCAGCAAAATGGGGCTGGGGATGCAGCAAAATTTCCCTTCCCCCCGCCGCAAGGCGTCACgg                                                                                                                                  <  2:2490684/147‑1 (MQ=255)
                     tCACCATTAACGTGGGACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACgg                                                                                                                                  >  1:1648319/1‑150 (MQ=255)
                     tCACCATTAACGTGGGACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACgg                                                                                                                                  >  2:234709/1‑150 (MQ=255)
                          aTTAACGTGGGACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAAGGGGGCTGTGGATGCAGCGAAATTTACCTCCACCCCGCCGCAAGGCGTCACGGTAgat                                                                                                                             <  1:2917206/150‑1 (MQ=255)
                           ttAACGTGGGACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAgatg                                                                                                                            >  1:957936/1‑150 (MQ=255)
                            tAACGTGGGACGTGATGGCACAATCCATCAGTTTAGCGGGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTCCCTTCCCCCCGCCGCAAGGCGTCACGGTAgatga                                                                                                                           <  1:3380150/150‑1 (MQ=255)
                            tAACGTGGGACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAgatga                                                                                                                           >  2:3479293/1‑150 (MQ=255)
                             aaggggggAGGGGATGGCAAAATCCATCATTTTGGCGGGGTGGAGCAGGACGATCGGCGCAGCATTTATCAACTAAAACCCCAGAAAAAGGGGGTTGGGGATGCAGCAAATTTTCCCTTCCCCCCGCCGCAAGGCGTCACGGTAgatgat                                                                                                                          <  2:2132770/145‑1 (MQ=255)
                              aCGTGGGACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATc                                                                                                                         >  1:1285251/1‑150 (MQ=255)
                                gTGGGACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCaa                                                                                                                       >  1:3449583/1‑150 (MQ=255)
                                  gggACGTGATGGCACAATCCATCAGTTTAGCGCGGGGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAACCCCAGAAAAAGGGGGCTGTGGATGCAGCGAAATTTCCCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACg                                                                                                                     <  2:3404531/150‑1 (MQ=255)
                                   ggACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTCCCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGt                                                                                                                    <  1:1681331/150‑1 (MQ=255)
                                   ggACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGt                                                                                                                    >  1:322876/1‑150 (MQ=255)
                                   ggACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGt                                                                                                                    >  1:4017534/1‑150 (MQ=255)
                                   ggACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGt                                                                                                                    >  1:3773024/1‑150 (MQ=255)
                                    gACGTGATGGCACAATCCATCAGTTTAGCGGGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAACCCCAGCAAAATGGGGCTGTGGATGCAGCAAATTTTACCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTa                                                                                                                   <  2:3729537/150‑1 (MQ=255)
                                      cGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAg                                                                                                                 >  2:1861252/1‑150 (MQ=255)
                                          aTGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTaga                                                                                                             >  2:1362617/1‑150 (MQ=255)
                                             gCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGc                                                                                                          >  1:2360796/1‑150 (MQ=255)
                                              cacaATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCa                                                                                                         >  2:3220908/1‑150 (MQ=255)
                                              cacaATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAACCCCAGAAAAATGGGGCTGTGGATGCAGCGAAATTTACCTCCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCa                                                                                                         <  2:2967347/150‑1 (MQ=255)
                                                 aaTCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCt                                                                                                      >  2:480550/1‑150 (MQ=255)
                                                 aaTCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCt                                                                                                      >  1:3209538/1‑150 (MQ=255)
                                                     cATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGAAAAAGGGGGCTGTGGAGGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCtgagt                                                                                                  <  1:1755492/150‑1 (MQ=255)
                                                         aGTTTAGCGGGGGGGAGAAGGACGATCAGCGCAGCATTTATCAACTAAAACCCCAGCAAAATGGGGCTGTGGATGCAGCGAAACTTCCCTCCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGc                                                                                              <  2:3331326/150‑1 (MQ=255)
                                                          gTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCa                                                                                             >  2:192157/1‑150 (MQ=255)
                                                            ttAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAAt                                                                                           <  1:1213098/150‑1 (MQ=255)
                                                             tAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATc                                                                                          >  2:947066/1‑150 (MQ=255)
                                                             tAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAAGGGGGCTGTGGATGCAGCGAAATTTACCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATc                                                                                          <  2:4205232/150‑1 (MQ=255)
                                                              aGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCt                                                                                         <  1:2821594/150‑1 (MQ=255)
                                                                 gcgGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTACACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTc                                                                                      <  2:3576783/150‑1 (MQ=255)
                                                                  cgGTGGAGCAGGAAGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCg                                                                                     >  2:2594651/1‑150 (MQ=255)
                                                                   ggTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGc                                                                                    >  1:278614/1‑150 (MQ=255)
                                                                   ggTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTAAAACCCCAGAAAAAGGGGGCTGTGGATGCAGCGAAATTTCCCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGc                                                                                    <  2:2848827/150‑1 (MQ=255)
                                                                    gTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCt                                                                                   <  1:300262/150‑1 (MQ=255)
                                                                     tGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTc                                                                                  >  2:360952/1‑150 (MQ=255)
                                                                      ggAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTCCCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCg                                                                                 <  1:3984278/150‑1 (MQ=255)
                                                                           aGGACGATCAGCGCAGCATTTATCAACTGAAACCCCAGCAAAAGGGGGCTGTGGATGCAGCGAAATTTACCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGAtttt                                                                            <  2:3346918/150‑1 (MQ=255)
                                                                           aGGACGATCAGCGCAGCAGTTAACAACCAAAATCCAAGAAAAAGGGGGCTGGGGAGGCGGCAAACTTTCCCTCCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGAtttt                                                                            <  2:2259682/150‑1 (MQ=255)
                                                                            ggACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTCCCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGAttttt                                                                           <  1:4071197/150‑1 (MQ=255)
                                                                            ggACGATCAGCGCAGCAGTTATCAACTAAAATCCCAGCAAAAGGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGAttttt                                                                           <  1:3319892/150‑1 (MQ=255)
                                                                             gACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTACACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTc                                                                          <  1:2042637/150‑1 (MQ=255)
                                                                               cGATCAGCGCAGCAGTTATCAACTGAAACCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCgg                                                                        <  1:1017795/150‑1 (MQ=255)
                                                                                  tCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGata                                                                     <  1:197885/150‑1 (MQ=255)
                                                                                    aGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGataat                                                                   >  1:212354/1‑150 (MQ=255)
                                                                                      cgcAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATAc                                                                 >  2:3033533/1‑150 (MQ=255)
                                                                                       gcagcaGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACt                                                                >  1:3709928/1‑150 (MQ=255)
                                                                                       gcagcaGTTATCAACTGAAACCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACt                                                                <  1:3539778/150‑1 (MQ=255)
                                                                                        cagcagTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACtt                                                               >  2:728065/1‑150 (MQ=255)
                                                                                         agcatTTATCAACTGAAACCCCAGAAAAATGGGGCTGTGGATGCAGCGAAATTTCCCTCCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACttt                                                              <  2:4380720/150‑1 (MQ=255)
                                                                                         agcagTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACttt                                                              >  1:3485653/1‑150 (MQ=255)
                                                                                            agTTATCAACTGAAATCCCGGCAAAAGGGGGCTGTGGATGCAGCGAAATTTCCCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCa                                                           <  1:1922298/150‑1 (MQ=255)
                                                                                            agTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCa                                                           <  1:2270433/150‑1 (MQ=255)
                                                                                            agTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCa                                                           <  1:2661191/150‑1 (MQ=255)
                                                                                             tttATCAACTGAAACCCCAGCAAAAGGGGGTTGTGGATGCGGCAAAATTTCCCTCCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCaa                                                          <  2:181318/149‑1 (MQ=255)
                                                                                              ttATAAACTGAAACCCCAGCAAAATGGGGGTGGGGATGCAGCGAAATTTCCCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAAc                                                         <  2:3858884/150‑1 (MQ=255)
                                                                                               tATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAAcc                                                        >  2:1629172/1‑150 (MQ=255)
                                                                                               tATCAACTAAAATCCCAGCAAAATGGGTGTGTGGATGCAGCGAATTTTACCTCCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAAcc                                                        <  2:1441875/150‑1 (MQ=255)
                                                                                                aTCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACct                                                       >  2:2700287/1‑150 (MQ=255)
                                                                                                    aCTAAAACCCCAGCAAAATGGGGCTGGGGATGCAGCGAAATTTCCCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTgg                                                   <  1:3856668/150‑1 (MQ=255)
                                                                                                       gAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCcg                                                >  2:1145571/1‑150 (MQ=255)
                                                                                                       gAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATATGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCcg                                                <  1:3911475/150‑1 (MQ=255)
                                                                                                        aaaTCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCcgc                                               >  1:2540212/1‑150 (MQ=255)
                                                                                                           tCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCcgcgcg                                            >  2:4008799/1‑150 (MQ=255)
                                                                                                               aGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATg                                        >  2:1683571/1‑150 (MQ=255)
                                                                                                                  aaaaTGGGGCTGTGGATGCAGCGAAATTTCCCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCgg                                     <  1:552539/150‑1 (MQ=255)
                                                                                                                    aaTGGGGCTGTGGATCCAGCGAAATTTACCTCCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGcc                                   <  2:3655521/150‑1 (MQ=255)
                                                                                                                     aTGGGGCTGTGGATGCAGCGAAATTTCCCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCa                                  <  2:3935909/150‑1 (MQ=255)
                                                                                                                     aTGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCa                                  >  2:2780844/1‑150 (MQ=255)
                                                                                                                        ggggTGTGGATGCAGCAAATTTTACTTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGaa                               <  1:2349698/150‑1 (MQ=255)
                                                                                                                        gggCTGTGGATGCAGCAAAATTTCCCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGaa                               <  2:4282114/150‑1 (MQ=255)
                                                                                                                         ggCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGaaa                              >  2:1221867/1‑150 (MQ=255)
                                                                                                                          gCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGCGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACATCTGGCCGCGCGTATGCGGCCAGaaaa                             >  2:4128512/1‑150 (MQ=255)
                                                                                                                          gCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGaaaa                             >  2:1321408/1‑150 (MQ=255)
                                                                                                                            tgtgGATGCAGCGAATTTTCCCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAAtt                           <  2:2527410/150‑1 (MQ=255)
                                                                                                                            tgtgGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAAtt                           >  1:4162422/1‑150 (MQ=255)
                                                                                                                               ggATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAg                        >  1:4178395/1‑150 (MQ=255)
                                                                                                                               ggATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAg                        >  1:4396327/1‑150 (MQ=255)
                                                                                                                               ggATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTAGGCGGCCAGAAAATTTAg                        >  2:1394130/1‑150 (MQ=255)
                                                                                                                                gATGCAGCGAATTTTCCCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGc                       <  1:1394024/150‑1 (MQ=255)
                                                                                                                                gATGCAGCGAAATTTCCCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGc                       <  2:2743790/150‑1 (MQ=255)
                                                                                                                                   gCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGcaca                    <  1:3072675/150‑1 (MQ=255)
                                                                                                                                    cAGGGAAATTTACCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAg                   <  1:1850023/150‑1 (MQ=255)
                                                                                                                                    cAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAg                   <  1:775131/150‑1 (MQ=255)
                                                                                                                                     aacaaaaTTTCCCTCCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGt                  <  2:3467567/146‑1 (MQ=255)
                                                                                                                                      gCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGta                 >  2:1012085/1‑150 (MQ=255)
                                                                                                                                      gCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGta                 >  2:2007923/1‑150 (MQ=255)
                                                                                                                                         aaaTTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTatat              >  2:2454066/1‑150 (MQ=255)
                                                                                                                                          attttACCTCCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATc             <  2:686867/148‑1 (MQ=255)
                                                                                                                                          aaTTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATc             <  1:4327946/150‑1 (MQ=255)
                                                                                                                                           aTTTACCTTCCCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCg            <  1:611160/150‑1 (MQ=255)
                                                                                                                                             ttcccTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCGGc          <  2:83424/147‑1 (MQ=255)
                                                                                                                                              tACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCGGcc         <  1:1114074/150‑1 (MQ=255)
                                                                                                                                                    cACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCGGCCAGCAAc   >  2:3166412/1‑150 (MQ=255)
                                                                                                                                                     aCCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACGGTATATCGGCCAGCAACa  >  2:2420699/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAATGGCAATCTGAAGCAGTTCACCATTAACGTGGGACGTGATGGCACAATCCATCAGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCGGCCAGCAACA  >  AP012306/800526‑800824

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