Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,944,2880AC35.7% 398.7 / 33.9 196G361G (GGT→GGGmenEo‑succinylbenzoate‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (103/23);  new base C (0/70);  total (103/93)
Fisher's exact test for biased strand distribution p-value = 2.09e-33
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 242 total aligned reads displayed.

ACCCATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCTGCCAT  >  AP012306/1944139‑1944437
                                                                                                                                                     |                                                                                                                                                     
aCCCATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATa                                                                                                                                                       <  1:2402240/150‑1 (MQ=255)
   cATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAccc                                                                                                                                                    >  1:2612026/1‑150 (MQ=255)
        aCTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTcc                                                                                                                                               <  2:512776/150‑1 (MQ=255)
         cTAAGATCGACGCTCTCGTGGTCATCCTCCACCCCCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCg                                                                                                                                              <  2:1238824/150‑1 (MQ=255)
          tAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGc                                                                                                                                             >  2:915586/1‑150 (MQ=255)
          tAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGc                                                                                                                                             >  2:1754033/1‑150 (MQ=255)
          tAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGc                                                                                                                                             >  2:1543089/1‑150 (MQ=255)
          tAAGATCAACGCTCTCGTGGTCAGCCTCCATCACCGCCACCGGTCGATGACCAAACGCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGc                                                                                                                                             <  2:3522382/150‑1 (MQ=255)
                cacccctctCGTGGTCATACTCCACCCCCGCCCCCGGCGGATGACCAACCCCCTTGCCGGCAGCGGGGACGATAAACCCCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGCCTTCTTCCGGCTAAATCCCCTCTCCGCCACTga                                                                                                                                       <  2:650131/145‑1 (MQ=255)
                   cGCTCTGGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTCCACCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTgaaga                                                                                                                                    <  2:2193949/150‑1 (MQ=255)
                       ctggTGGTCATCCTCCACCCCCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAg                                                                                                                                <  2:2566887/147‑1 (MQ=255)
                           tGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTg                                                                                                                            >  2:2310450/1‑150 (MQ=255)
                           tGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTg                                                                                                                            >  2:472494/1‑150 (MQ=255)
                             gTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTc                                                                                                                          >  2:243889/1‑150 (MQ=255)
                               cATACTCCATCACCTCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTa                                                                                                                        >  1:3528214/1‑150 (MQ=255)
                               cATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTa                                                                                                                        <  1:1466749/150‑1 (MQ=255)
                               cATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTa                                                                                                                        >  1:3514379/1‑150 (MQ=255)
                               cATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTa                                                                                                                        >  2:1125294/1‑150 (MQ=255)
                                aTACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTaa                                                                                                                       >  1:870698/1‑150 (MQ=255)
                                   ccccATCCCCGCCCCCGGTCGATGACCAACCTCCTTGCCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTCCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACg                                                                                                                    <  2:3761800/148‑1 (MQ=255)
                                     ccATCCCCGCCACCGGTCGATGACCAACCTCCTTGTCGGAAACGGGGACGATAAACACCTGAAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCGGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTc                                                                                                                  <  2:2655213/150‑1 (MQ=255)
                                     ccATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTc                                                                                                                  >  1:3805238/1‑150 (MQ=255)
                                      cATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc                                                                                                                 <  1:2025335/150‑1 (MQ=255)
                                          aCCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACa                                                                                                             >  1:2688900/1‑150 (MQ=255)
                                             gCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATg                                                                                                          >  1:2778219/1‑150 (MQ=255)
                                                 ccGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCa                                                                                                      >  1:154566/1‑150 (MQ=255)
                                                   ggTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGc                                                                                                    <  1:3613058/150‑1 (MQ=255)
                                                     tggATGACCAATCTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAACCCGCAGGATGTGCACCAATTACGCGCTCGACTTCTTCCGGCTGAACACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCtt                                                                                                  <  2:1813168/148‑1 (MQ=255)
                                                     tCGATGACCAAACTCCTTGTCGGCAACGGGGATGATAAACACCTGAAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCtt                                                                                                  <  1:2880203/150‑1 (MQ=255)
                                                     tCGATGACCAAACCCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCtt                                                                                                  <  2:3104918/150‑1 (MQ=255)
                                                      cGATGACCAACCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTg                                                                                                 <  1:176635/150‑1 (MQ=255)
                                                      cGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTg                                                                                                 <  2:3365123/150‑1 (MQ=255)
                                                      cGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTg                                                                                                 <  1:588667/150‑1 (MQ=255)
                                                       gATGACCAAACTCCTTGCCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGc                                                                                                <  2:2536839/150‑1 (MQ=255)
                                                        aTGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGcc                                                                                               >  1:1231003/1‑150 (MQ=255)
                                                              aatcTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGAAATTACGCGCTCGACTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATg                                                                                         <  2:382216/147‑1 (MQ=255)
                                                              aaccTCCTTGTCGGCAACGGGGACGATAAACCCCTGAAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATg                                                                                         <  2:3871062/147‑1 (MQ=255)
                                                              aaaCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATg                                                                                         >  1:510262/1‑150 (MQ=255)
                                                                aCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCa                                                                                       <  2:2401585/150‑1 (MQ=255)
                                                                    cTTGTCGGCAACGGGGACGATAAACACCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTAAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtctc                                                                                   <  2:819426/150‑1 (MQ=255)
                                                                      tGTCGGCAACGGGGACGATAAACGCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCAc                                                                                 <  1:1509444/150‑1 (MQ=255)
                                                                      tGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCAc                                                                                 <  1:227547/150‑1 (MQ=255)
                                                                      tGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCAc                                                                                 >  2:126614/1‑150 (MQ=255)
                                                                      tGTCGGCAACGGGGACGAAAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCAc                                                                                 <  2:1417819/150‑1 (MQ=255)
                                                                         cgtcAACGGGCGCGATAATCACGCGGAAACCCGAGGGATGTGAAGTATTTACCCGCCCGCCTTTTTCGGTCGGACTACCTTTCCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACcgc                                                                              <  2:3279182/147‑1 (MQ=255)
                                                                            cAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTCCGCGCTCGACTTCTTCCGGTTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGAt                                                                           <  1:3888377/150‑1 (MQ=255)
                                                                              aCGGGGACGATAACCCCCTGCAAACCCGCAGGATGTCCAGCAATTCCGCGCTCGCCTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATcg                                                                         <  2:2763347/150‑1 (MQ=255)
                                                                              aCGGGGACGATAAACACCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATcg                                                                         <  1:2332365/150‑1 (MQ=255)
                                                                              aCGGGGACGATAAACACCTGCAAACCCGAAGGATGTGCAGCAATTACGCGCTCGACTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATcg                                                                         <  2:4000146/150‑1 (MQ=255)
                                                                                 gggACGATAAACACCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATcgcgc                                                                      <  2:2598154/150‑1 (MQ=255)
                                                                                   gACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgt                                                                    >  1:2178858/1‑150 (MQ=255)
                                                                                       aTAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtagcg                                                                <  1:1657328/150‑1 (MQ=255)
                                                                                         aaaCACCTGCAAAACCGCAGGTTGTGCAGCAATTACGCGCTCGTCTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtagcgta                                                              <  2:1898603/150‑1 (MQ=255)
                                                                                           acacCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTAcc                                                            >  1:313167/1‑150 (MQ=255)
                                                                                           acacCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTAcc                                                            >  1:237331/1‑150 (MQ=255)
                                                                                             acCTGAAAAACCGCAGGATGTGCAGCAATTACGCGCTCGCCTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAg                                                          <  2:3152069/150‑1 (MQ=255)
                                                                                               cTGCAAAACCGCAGGATGTGCAGCAATTCCGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGcc                                                        <  1:3278122/150‑1 (MQ=255)
                                                                                                tccAAAACCTGAGGATGTGCAGCAATTACGCGCTCGACTTCTGCCGGCGAAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCt                                                       <  1:453562/148‑1 (MQ=255)
                                                                                                tGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCt                                                       >  1:131711/1‑150 (MQ=255)
                                                                                                   aaacCCGCAGGATGTGCAGCAATTCCGCGCTCGCCTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCg                                                    <  2:145583/150‑1 (MQ=255)
                                                                                                   aaaaCCGCAGGATGTGCACCAATTCCGCGCTCGACTTCTCCCGGCTGAACACCCCCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCg                                                    <  2:1053349/150‑1 (MQ=255)
                                                                                                    aacccGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGt                                                   <  2:671962/147‑1 (MQ=255)
                                                                                                    aaaCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGt                                                   >  1:981960/1‑150 (MQ=255)
                                                                                                    aaaCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGt                                                   >  2:2038126/1‑150 (MQ=255)
                                                                                                      aCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCa                                                 >  2:4104520/1‑150 (MQ=255)
                                                                                                      aCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCa                                                 >  1:1946977/1‑150 (MQ=255)
                                                                                                      aCCGCAGGATGTGCAGAAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCa                                                 <  1:1813442/150‑1 (MQ=255)
                                                                                                         gCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTa                                              >  1:2044278/1‑150 (MQ=255)
                                                                                                           aGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAAc                                            <  1:3618342/150‑1 (MQ=255)
                                                                                                           aGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAAc                                            >  1:2499955/1‑150 (MQ=255)
                                                                                                             gATGTGCATCAATTACGCGCTCCTCTTCTTCCGGCTAAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCa                                          <  2:2351329/150‑1 (MQ=255)
                                                                                                             gATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCa                                          >  1:3782187/1‑150 (MQ=255)
                                                                                                              aTGTTCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAAGAGATTGTCTAAACTTCCGACAATGGTCAGCTTGCCATTATGCATCTCTCCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAg                                         >  2:4094335/1‑150 (MQ=255)
                                                                                                              aTGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAg                                         <  1:3343452/150‑1 (MQ=255)
                                                                                                                  gcagcaATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGaa                                     <  1:1057394/150‑1 (MQ=255)
                                                                                                                   cagcaATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGaaa                                    >  1:3165905/1‑150 (MQ=255)
                                                                                                                   cagcaATTACGCGCTCGACTTCTTCCGGCGGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGaaa                                    <  1:633703/150‑1 (MQ=255)
                                                                                                                     gcaATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAAcc                                  >  1:818659/1‑150 (MQ=255)
                                                                                                                     gcaATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAAcc                                  >  2:3208585/1‑150 (MQ=255)
                                                                                                                      caATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCa                                 >  1:3955009/1‑150 (MQ=255)
                                                                                                                      caATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCa                                 >  2:3004291/1‑150 (MQ=255)
                                                                                                                       aaTTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAg                                >  2:3309951/1‑150 (MQ=255)
                                                                                                                         ttACGGGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGtt                              <  1:2836989/150‑1 (MQ=255)
                                                                                                                         ttACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGtt                              >  2:1961250/1‑150 (MQ=255)
                                                                                                                             gcgcTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGccc                          <  1:3546984/150‑1 (MQ=255)
                                                                                                                                cTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGtt                       <  1:246494/150‑1 (MQ=255)
                                                                                                                                  cGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTAc                     <  2:1368171/150‑1 (MQ=255)
                                                                                                                                  cGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTAc                     >  1:103559/1‑150 (MQ=255)
                                                                                                                                  cGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTAc                     >  2:1678217/1‑150 (MQ=255)
                                                                                                                                  cGACTTCTTCCGGCGGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTAc                     <  1:1172901/150‑1 (MQ=255)
                                                                                                                                   gACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACg                    <  1:332483/150‑1 (MQ=255)
                                                                                                                                   gACTTCTTCCGGCTAAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACg                    <  1:2546753/150‑1 (MQ=255)
                                                                                                                                     cttcttCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTCCCGTTACGCCa                 <  2:1577520/150‑1 (MQ=255)
                                                                                                                                       tcttcCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAg                <  1:2772375/150‑1 (MQ=255)
                                                                                                                                        cttcCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGt               <  1:646813/150‑1 (MQ=255)
                                                                                                                                             ggCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCt          >  1:2250842/1‑150 (MQ=255)
                                                                                                                                             ggCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCt          >  2:2006964/1‑150 (MQ=255)
                                                                                                                                                ggAATCCCCCCCCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTCCCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCt       <  2:376039/149‑1 (MQ=255)
                                                                                                                                                  aaTCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCTGc     <  2:289874/150‑1 (MQ=255)
                                                                                                                                                     aCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCTGCCAt  >  2:46075/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 < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: