Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,944,2880AC38.2% 348.5 / 59.3 206G361G (GGT→GGGmenEo‑succinylbenzoate‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (90/36);  new base C (0/78);  total (90/116)
Fisher's exact test for biased strand distribution p-value = 2.44e-28
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 280 total aligned reads displayed.

ACCCATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCTGCCAT  >  AP012306/1944139‑1944437
                                                                                                                                                     |                                                                                                                                                     
aCCCATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATa                                                                                                                                                       >  1:2423768/1‑150 (MQ=255)
  ccATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAcc                                                                                                                                                     <  2:3371846/150‑1 (MQ=255)
   cATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAccc                                                                                                                                                    <  2:4107173/150‑1 (MQ=255)
    aTTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCct                                                                                                                                                   <  2:3823271/150‑1 (MQ=255)
        aCTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTcc                                                                                                                                               <  2:1655695/150‑1 (MQ=255)
            aGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCa                                                                                                                                           >  1:3657390/1‑150 (MQ=255)
              accAACGCTCTCGTGGTCATACCCCATCCCCGCCCCCGGTCAATGACCAACCTCTTTGCCGGCAACGGGGACGATAAACCCCTCCAAAACCGCAGGATGTCCAGCAATTCCCCGCCCGACTTCTCCCGGCTAAATCCCCTCTCCGCCACt                                                                                                                                         <  2:1219688/148‑1 (MQ=255)
               ccatcGCTCTCGTGATCATCCTCCACCACCCCCCCCGGTCGATGCCCAACCTCCTTGACGCCAAGGGGGACGATAAACTCCTCCAAATCCGCAGGATGTGCAGCAATTTCGCGCTGGACTTCTTCCGGTTAAATACCCTCTCCGCCACTg                                                                                                                                        <  2:912534/146‑1 (MQ=255)
                cAACCCTCTCGTGGTCATACCCCATCCCCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTga                                                                                                                                       <  2:1750583/150‑1 (MQ=255)
                 aaCGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTgaa                                                                                                                                      >  2:548351/1‑150 (MQ=255)
                 aaCGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTgaa                                                                                                                                      >  1:22992/1‑150 (MQ=255)
                  aCGCTCTCGTGGTCATACTCCATCCCCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTgaag                                                                                                                                     <  1:2608079/150‑1 (MQ=255)
                  aCGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTgaag                                                                                                                                     <  1:1546602/150‑1 (MQ=255)
                   cGCTCTCGTGGTCATACTCCATCCCCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTgaaga                                                                                                                                    <  2:600159/150‑1 (MQ=255)
                   cGCCCTCGTGTTCATCCCCCATCCCCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTgaaga                                                                                                                                    <  2:342093/150‑1 (MQ=255)
                    gCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTgaagaa                                                                                                                                   >  1:3727825/1‑150 (MQ=255)
                      tctcGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATa                                                                                                                                 >  1:3231599/1‑150 (MQ=255)
                      tctcGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATa                                                                                                                                 >  1:268931/1‑150 (MQ=255)
                      tctcGTGGTCATACTCCATCACCGCCACCGGTCGATGACAAACCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATa                                                                                                                                 <  1:3774171/150‑1 (MQ=255)
                      tctcGTGGTCATACTCCACCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATa                                                                                                                                 <  2:704662/150‑1 (MQ=255)
                       ctcGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAg                                                                                                                                <  1:1126441/150‑1 (MQ=255)
                         cGTGGTCATCCTCCATCCCCGCCACCGGTCGATGACAAAACTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt                                                                                                                              <  2:491626/150‑1 (MQ=255)
                           tGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTg                                                                                                                            >  2:1075482/1‑150 (MQ=255)
                            ggTCATACTCCATCACCGCCACCGGTCGATGACCAACCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGt                                                                                                                           <  1:2054819/150‑1 (MQ=255)
                            ggTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGt                                                                                                                           >  2:3169052/1‑150 (MQ=255)
                               cATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTa                                                                                                                        >  1:498135/1‑150 (MQ=255)
                                atcccccATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAACCACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGATTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTaa                                                                                                                       <  2:2776758/145‑1 (MQ=255)
                                aTACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTaa                                                                                                                       >  2:648073/1‑150 (MQ=255)
                                 tACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTaaa                                                                                                                      <  1:105157/150‑1 (MQ=255)
                                  aCTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAAc                                                                                                                     >  2:1169145/1‑150 (MQ=255)
                                  aCTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAAc                                                                                                                     >  1:1650577/1‑150 (MQ=255)
                                   cTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATGCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACg                                                                                                                    <  1:2870740/150‑1 (MQ=255)
                                        tCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGa                                                                                                               <  1:3521993/150‑1 (MQ=255)
                                         ccccGCCCCCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAACCCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTAAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGAc                                                                                                              <  2:3737822/148‑1 (MQ=255)
                                          aCCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTAAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACa                                                                                                             <  2:3193801/150‑1 (MQ=255)
                                            cGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAAt                                                                                                           <  2:3999166/150‑1 (MQ=255)
                                            cGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAAt                                                                                                           >  2:1411533/1‑150 (MQ=255)
                                             gCCCCCGGTCGATGACCAACCTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAACCCGCAGGATGGGCAGCAATTACGCGCTCGACTTCTCCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATg                                                                                                          <  2:3655260/150‑1 (MQ=255)
                                              ccACCGGTCGATGACCAACCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGATTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATgg                                                                                                         <  1:3407617/150‑1 (MQ=255)
                                                 ccGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGAAAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCa                                                                                                      <  2:2895176/150‑1 (MQ=255)
                                                     ccGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCtt                                                                                                  <  2:3436288/149‑1 (MQ=255)
                                                      cgttgcccAACCCCCTTGTTGGCACCGGGGAACATAAACTCCTGCAAAACCGCAGGATGGGCAATAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCCCCCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTg                                                                                                 <  2:3292193/144‑1 (MQ=255)
                                                      cGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTg                                                                                                 >  1:1388166/1‑150 (MQ=255)
                                                            ccAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTa                                                                                           >  2:2243504/1‑150 (MQ=255)
                                                             cAAACTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTCCCGGCTAAATCCCCTCCCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTAt                                                                                          <  2:234304/150‑1 (MQ=255)
                                                             cAAACCCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTAt                                                                                          <  2:779418/150‑1 (MQ=255)
                                                                     ttGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCa                                                                                  <  1:2125276/150‑1 (MQ=255)
                                                                       gTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCAcc                                                                                >  1:1215985/1‑150 (MQ=255)
                                                                        tCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACcg                                                                               >  2:1793671/1‑150 (MQ=255)
                                                                        tCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGCCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACcg                                                                               <  2:3137675/150‑1 (MQ=255)
                                                                         cGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGCCTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACcgc                                                                              <  2:190844/150‑1 (MQ=255)
                                                                           gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGa                                                                            >  2:3034661/1‑150 (MQ=255)
                                                                               cGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATcgc                                                                        >  1:695483/1‑150 (MQ=255)
                                                                                gggaatgATAAACCCCTGCAAAACGGAAGGTTGTGAAGAAATTACGCGTTCCATTTCTTCCGGTGGAATTCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATcgcg                                                                       <  2:62693/144‑1 (MQ=255)
                                                                                 gggACGATAAACCCCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGTCTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATcgcgc                                                                      <  2:2506680/150‑1 (MQ=255)
                                                                                   gACGAAAACCCCCTGAAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTTTTCGGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgt                                                                    <  2:3527753/150‑1 (MQ=255)
                                                                                    aCGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgta                                                                   >  1:4399681/1‑150 (MQ=255)
                                                                                     cGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTCCGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag                                                                  <  1:171910/150‑1 (MQ=255)
                                                                                                tGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCt                                                       <  2:686532/150‑1 (MQ=255)
                                                                                                  cAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTc                                                     >  2:1849011/1‑150 (MQ=255)
                                                                                                      aCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCa                                                 <  1:184987/150‑1 (MQ=255)
                                                                                                       ccGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCAt                                                >  1:2801420/1‑150 (MQ=255)
                                                                                                       ccGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCAt                                                >  1:1782915/1‑150 (MQ=255)
                                                                                                       ccGCAGGATGTGCAGCAATTACGCCCTCGACTTCTTCCGGCTGAATCCCCCCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCAt                                                <  2:1540700/150‑1 (MQ=255)
                                                                                                        cGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCAtt                                               >  2:4328089/1‑150 (MQ=255)
                                                                                                        cGCAGGAGGTGCGCCATTTACGCGCTCGCCTTCTTCCGACTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCAtt                                               <  2:2128516/150‑1 (MQ=255)
                                                                                                         gCAGGATGTGCAGCAATTCCGCTCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTa                                              <  1:266212/150‑1 (MQ=255)
                                                                                                         gCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTa                                              >  1:2543774/1‑150 (MQ=255)
                                                                                                            ggATGTGCAGCAATTACGCGTTCGACTTCTCCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAAcc                                           <  1:5310/150‑1 (MQ=255)
                                                                                                               tgtgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGt                                        >  2:4051669/1‑150 (MQ=255)
                                                                                                               tgtgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGt                                        <  1:3289039/150‑1 (MQ=255)
                                                                                                                gtgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTg                                       >  1:2333217/1‑150 (MQ=255)
                                                                                                                 tgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGa                                      <  1:427814/150‑1 (MQ=255)
                                                                                                                  gcagcaATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGaa                                     <  2:1553300/150‑1 (MQ=255)
                                                                                                                  gcagcaATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGaa                                     <  2:1286331/150‑1 (MQ=255)
                                                                                                                    agcaATTACGCGCTCGACTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAAc                                   <  2:3114779/150‑1 (MQ=255)
                                                                                                                    agcaATTACGCGCTCGACTTCTCCCGGCTAAATCCCCCCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAAc                                   <  2:86006/150‑1 (MQ=255)
                                                                                                                      taaTTACGCGCTCGCTTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCa                                 <  2:388494/149‑1 (MQ=255)
                                                                                                                      caATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCa                                 <  1:861578/150‑1 (MQ=255)
                                                                                                                       aaTTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAg                                <  1:1264518/150‑1 (MQ=255)
                                                                                                                        aTTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGt                               >  2:4246064/1‑150 (MQ=255)
                                                                                                                         ttACGCGCTCGACTTCTTCCGGCGGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGtt                              <  1:1915935/150‑1 (MQ=255)
                                                                                                                            cgcgCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGcc                           >  1:2891037/1‑150 (MQ=255)
                                                                                                                            cgcgCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGCCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTccc                           >  2:3749093/1‑150 (MQ=255)
                                                                                                                             gcgcTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGccc                          <  1:3006617/150‑1 (MQ=255)
                                                                                                                             gcgcTCGACTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGccc                          <  1:3199527/150‑1 (MQ=255)
                                                                                                                              cgcTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCg                         >  2:4228992/1‑150 (MQ=255)
                                                                                                                                cTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGtt                       >  1:4196653/1‑150 (MQ=255)
                                                                                                                                  cGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGCTGCCCGTTAc                     >  1:3117984/1‑150 (MQ=255)
                                                                                                                                      ttcttcCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCa                 >  2:2625205/1‑150 (MQ=255)
                                                                                                                                         ttcCGGCTGAATGCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTTCCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTa              <  1:4047593/150‑1 (MQ=255)
                                                                                                                                         ttcCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTa              >  1:3868004/1‑150 (MQ=255)
                                                                                                                                          tcCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTaa             <  1:2242993/150‑1 (MQ=255)
                                                                                                                                           ccGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAAc            >  1:118521/1‑150 (MQ=255)
                                                                                                                                            cGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAAcc           <  1:596704/150‑1 (MQ=255)
                                                                                                                                            cGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAAcc           <  1:4101682/150‑1 (MQ=255)
                                                                                                                                             ggCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCt          <  1:2202611/150‑1 (MQ=255)
                                                                                                                                               cTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTc        >  1:2962081/1‑150 (MQ=255)
                                                                                                                                                   aTACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCTGcc    >  1:1969756/1‑150 (MQ=255)
                                                                                                                                                    tACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCTGCCa   <  1:3612996/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCCATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCTGCCAT  >  AP012306/1944139‑1944437

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