Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,944,2880AC35.1% 446.8 / 53.0 225G361G (GGT→GGG) menEo‑succinylbenzoate‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (108/38);  new base C (0/79);  total (108/117)
Fisher's exact test for biased strand distribution p-value = 1.34e-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 286 total aligned reads displayed.

ACCCATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCTGCCAT  >  AP012306/1944139‑1944437
                                                                                                                                                     |                                                                                                                                                     
aCCCATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATa                                                                                                                                                       >  1:3125885/1‑150 (MQ=255)
 cccATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAc                                                                                                                                                      <  1:335522/150‑1 (MQ=255)
    aTTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCct                                                                                                                                                   <  2:2843902/150‑1 (MQ=255)
     ttCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCctc                                                                                                                                                  <  2:2438920/150‑1 (MQ=255)
      tCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCctct                                                                                                                                                 <  1:3680690/150‑1 (MQ=255)
      tCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCctct                                                                                                                                                 <  1:3844725/150‑1 (MQ=255)
        aCTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTcc                                                                                                                                               <  1:4423964/150‑1 (MQ=255)
         cTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCg                                                                                                                                              <  1:1585877/150‑1 (MQ=255)
           aaGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGcc                                                                                                                                            <  2:1715936/150‑1 (MQ=255)
            aGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCa                                                                                                                                           >  2:2363705/1‑150 (MQ=255)
             gATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCAc                                                                                                                                          >  1:1465250/1‑150 (MQ=255)
             gATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCAc                                                                                                                                          >  1:4070008/1‑150 (MQ=255)
                cAACGCTCTCGTGGTCATCCTCCATCCCCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTga                                                                                                                                       <  2:2304404/150‑1 (MQ=255)
                cAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTga                                                                                                                                       >  1:2961494/1‑150 (MQ=255)
                cAACGCTCTCGGGGTCATCCTCCCCCCCCACCACCGGTCGAGGACCAACCTCCTTGTCGCCAACGGGGACTAAAAACCCCTGCAAAACCGCAGGATGTGCACCATTTACGCTCTAGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTga                                                                                                                                       <  2:498158/150‑1 (MQ=255)
                 aaCGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTgaa                                                                                                                                      >  1:755510/1‑150 (MQ=255)
                   cGCTCTCGTGGTCATACCCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTgaaga                                                                                                                                    <  2:1600758/150‑1 (MQ=255)
                    gCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTgaagaa                                                                                                                                   >  2:2918692/1‑150 (MQ=255)
                     ctctCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAAt                                                                                                                                  >  2:2541503/1‑150 (MQ=255)
                      tctcGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATa                                                                                                                                 >  1:2213312/1‑150 (MQ=255)
                         ggTGGTCCTCCCCCATCCCCGCCTCCGGTCGATGACCAATCTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGAGGTGCAGCAATTACGCGTTCAACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt                                                                                                                              <  2:513803/149‑1 (MQ=255)
                         cGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt                                                                                                                              <  1:1175880/150‑1 (MQ=255)
                           tGGTCATCCTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTg                                                                                                                            <  1:3062284/150‑1 (MQ=255)
                              tCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCt                                                                                                                         >  1:943587/1‑150 (MQ=255)
                               cATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTa                                                                                                                        <  1:1424961/150‑1 (MQ=255)
                                 tACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTATCGGCAACGGGGACGATAAACACCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTaaa                                                                                                                      <  1:1878339/150‑1 (MQ=255)
                                    tCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGt                                                                                                                   >  2:2380834/1‑150 (MQ=255)
                                      cATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc                                                                                                                 >  1:3469299/1‑150 (MQ=255)
                                         cACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGAc                                                                                                              >  1:1924125/1‑150 (MQ=255)
                                            cGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAAt                                                                                                           <  1:2824213/150‑1 (MQ=255)
                                             gCCCCCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATg                                                                                                          <  1:1544314/150‑1 (MQ=255)
                                               cACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGt                                                                                                        >  2:447048/1‑150 (MQ=255)
                                                 ccGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCa                                                                                                      >  2:3901775/1‑150 (MQ=255)
                                                  cGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTAAGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAg                                                                                                     <  1:693303/150‑1 (MQ=255)
                                                            ccAAACTCCTTGTCGGAAACGGGGACGATAAACCCCTGCAAAACCGCGGGATGTGCAGCATTTCCGCCCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTa                                                                                           <  2:769957/150‑1 (MQ=255)
                                                             aaaaCTCCTTCTCGACAACGGGGACGATAACCACCTGCAAACCCGGGGGATGTGCAGCAATCACGCGCTCGACTTCTCCGGGCAAATTCCCCCCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTAt                                                                                          <  2:1857962/149‑1 (MQ=255)
                                                              aaccTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATg                                                                                         <  1:2863333/147‑1 (MQ=255)
                                                              aaccTCCTTGTCGGCAACGGGGACGAAAAACCCCTGAAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCCCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATg                                                                                         <  2:2986586/147‑1 (MQ=255)
                                                              aaaCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAACCCGCAGGATGTGCAGCAATTCCGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATg                                                                                         <  2:2049221/150‑1 (MQ=255)
                                                                 cTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAt                                                                                      >  1:1640160/1‑150 (MQ=255)
                                                                   ccTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtct                                                                                    >  2:2882196/1‑150 (MQ=255)
                                                                     ttGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCa                                                                                  <  1:123372/150‑1 (MQ=255)
                                                                       gTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCAcc                                                                                <  2:2653217/150‑1 (MQ=255)
                                                                       gTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCAcc                                                                                >  1:2019160/1‑150 (MQ=255)
                                                                        tCGGAAACGGGGACGATAAACCCCTGAAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTCCCGGTTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACcg                                                                               <  2:2740971/150‑1 (MQ=255)
                                                                         cGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACcgc                                                                              <  2:2248589/150‑1 (MQ=255)
                                                                         cGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACcgc                                                                              >  1:3653747/1‑150 (MQ=255)
                                                                            cAACGGGGACGATAAACACCTGAAAAACCGCAGGATGTGCAGCAATTACGCGCTCGCCTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGAt                                                                           <  2:1891870/150‑1 (MQ=255)
                                                                             aaCGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATc                                                                          >  1:2488606/1‑150 (MQ=255)
                                                                                    aCGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgta                                                                   <  2:1523620/150‑1 (MQ=255)
                                                                                    aCGAAAACCACCTGAAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgta                                                                   <  2:1547878/150‑1 (MQ=255)
                                                                                     cGATAACCCCCTGCAAACCCGCAGGATGTGCAGCAATTCCGCGCTCGCCTCCTCCCGGCTGAATCCCCCCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag                                                                  <  2:195936/150‑1 (MQ=255)
                                                                                      gATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACCCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtagc                                                                 <  2:3736714/150‑1 (MQ=255)
                                                                                      gATAAACACCTGAAAAACCGCAGCATGTGCACAAATTACGGGCTCGACTACTCCCGGCTCAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtagc                                                                 <  2:3466875/150‑1 (MQ=255)
                                                                                        tAAACACCTGAAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtagcgt                                                               <  2:1202944/150‑1 (MQ=255)
                                                                                         aaaCCCCTCCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtagcgta                                                              <  2:1389814/150‑1 (MQ=255)
                                                                                          accacCTGAAAACCCGCAGGATGTGCAGCAATTACGCGCTCGCCTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTAc                                                             <  2:2169700/148‑1 (MQ=255)
                                                                                             acCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAg                                                          <  1:2267122/150‑1 (MQ=255)
                                                                                               cTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGcc                                                        >  1:3563637/1‑150 (MQ=255)
                                                                                                tGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCt                                                       <  2:1265788/150‑1 (MQ=255)
                                                                                                tGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCt                                                       >  1:1355572/1‑150 (MQ=255)
                                                                                                 gCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCtt                                                      <  2:3122844/150‑1 (MQ=255)
                                                                                                  cAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTc                                                     >  1:1006166/1‑150 (MQ=255)
                                                                                                   aaaaCCGCAGGATGTGCAGCAATTACGCGCTCGATTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCg                                                    <  2:3034746/150‑1 (MQ=255)
                                                                                                     acccGAAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATCCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTc                                                  <  1:4098972/148‑1 (MQ=255)
                                                                                                     aaCCGCAGGCTGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAAAGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTc                                                  >  1:3489582/1‑150 (MQ=255)
                                                                                                       ccGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCAt                                                >  1:435555/1‑150 (MQ=255)
                                                                                                       ccGCAGGAGGTGCAGAAATTACGCTCTCGTCTTCTTCCGGCGGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCAt                                                <  2:1647980/150‑1 (MQ=255)
                                                                                                         gCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTa                                              >  1:1905758/1‑150 (MQ=255)
                                                                                                          cAGGATGTGCAGCAATTACGCGCTCTACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTaa                                             <  2:3216395/150‑1 (MQ=255)
                                                                                                            ggATGTGCAGCAATTCCGCGCTCGCCTTCTTCCGGCTAAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAAcc                                           <  2:2116380/150‑1 (MQ=255)
                                                                                                              gtgggcaccaATCGCGCGCTCGCTTTTTTCCGGCGGAACTCCTTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAg                                         <  2:7020/146‑1 (MQ=255)
                                                                                                              aggtgAAGCAATTACGACGTGCACTTCTCCGGGCTGAATACCAGCCCCGCCACTGAAGAATAGATTGTCAAAACGACCGACAATGGTCAGCTTGCCATTATGCATCTCTCCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAg                                         <  2:3886170/148‑1 (MQ=255)
                                                                                                                  gcagcaATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGaa                                     <  1:3197627/150‑1 (MQ=255)
                                                                                                                    agcaATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAAc                                   <  1:1774512/150‑1 (MQ=255)
                                                                                                                     gcaATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAAcc                                  >  1:4158612/1‑150 (MQ=255)
                                                                                                                      caATTACGCGCTCGACTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCa                                 <  1:4263574/150‑1 (MQ=255)
                                                                                                                       aaTTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAg                                <  1:2585827/150‑1 (MQ=255)
                                                                                                                       aaTTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAg                                >  1:2076868/1‑150 (MQ=255)
                                                                                                                          tACGCGCTCGACTCCTCCCGGCGGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTg                             <  1:3412708/150‑1 (MQ=255)
                                                                                                                              cgcTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCg                         <  1:1704952/150‑1 (MQ=255)
                                                                                                                              cgcTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCg                         >  2:1794537/1‑150 (MQ=255)
                                                                                                                              cgcTCGACTTCTTCCGGCTAAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCg                         <  2:4043848/150‑1 (MQ=255)
                                                                                                                               gcTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGt                        >  1:1041195/1‑150 (MQ=255)
                                                                                                                                cTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGtt                       >  1:233764/1‑150 (MQ=255)
                                                                                                                                  cGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTAc                     <  1:3931977/150‑1 (MQ=255)
                                                                                                                                    aCTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGc                   >  2:3815724/1‑150 (MQ=255)
                                                                                                                                    aCTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGc                   >  1:3891194/1‑150 (MQ=255)
                                                                                                                                      ttcttcCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCa                 >  2:4373295/1‑150 (MQ=255)
                                                                                                                                      ttctccCGGCTGAATCCCCTCCCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCa                 <  2:3390306/150‑1 (MQ=255)
                                                                                                                                        cttcCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGt               >  1:849670/1‑150 (MQ=255)
                                                                                                                                         ttcCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCTGTTACGCCAGTa              >  2:3517134/1‑150 (MQ=255)
                                                                                                                                            cGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAAcc           >  1:471513/1‑150 (MQ=255)
                                                                                                                                            cGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAAcc           >  1:3263927/1‑150 (MQ=255)
                                                                                                                                               cTGAATCCCCTCTCCGCCACTGAAAAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTc        <  1:2707283/150‑1 (MQ=255)
                                                                                                                                               cTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTc        >  2:4263873/1‑150 (MQ=255)
                                                                                                                                                tGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCt       <  1:559448/150‑1 (MQ=255)
                                                                                                                                                 gAATCCCCCTTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCTg      <  2:3322328/150‑1 (MQ=255)
                                                                                                                                                  aaTACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCTGc     >  2:1338740/1‑150 (MQ=255)
                                                                                                                                                  aaTACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCTGc     >  2:1091053/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCCATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCTGCCAT  >  AP012306/1944139‑1944437

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