Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,944,2880AC36.5% 372.7 / 59.2 211G361G (GGT→GGGmenEo‑succinylbenzoate‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (99/35);  new base C (0/77);  total (99/112)
Fisher's exact test for biased strand distribution p-value = 2.25e-30
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 282 total aligned reads displayed.

ACCCATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCTG  >  AP012306/1944139‑1944433
                                                                                                                                                     |                                                                                                                                                 
aCCCATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATa                                                                                                                                                   >  2:888735/1‑150 (MQ=255)
aCCCATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATa                                                                                                                                                   >  1:2530456/1‑150 (MQ=255)
 cccATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAc                                                                                                                                                  <  2:4208391/150‑1 (MQ=255)
   cATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAccc                                                                                                                                                <  1:1308909/150‑1 (MQ=255)
   cATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAccc                                                                                                                                                <  2:305065/150‑1 (MQ=255)
      tCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCctct                                                                                                                                             >  1:211995/1‑150 (MQ=255)
       cACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCtctc                                                                                                                                            >  1:1351510/1‑150 (MQ=255)
       cACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCtctc                                                                                                                                            >  2:2826996/1‑150 (MQ=255)
       cACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCtctc                                                                                                                                            <  2:2547692/150‑1 (MQ=255)
            aGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCa                                                                                                                                       <  2:702569/150‑1 (MQ=255)
             gATCACCGCTCTCGTGGTCATCCTCCACCCCCCCCCCCGGTCGATGACCAACCCCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCAc                                                                                                                                      <  2:4559989/150‑1 (MQ=255)
             gATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCAc                                                                                                                                      <  1:2021975/150‑1 (MQ=255)
               tCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg                                                                                                                                    >  1:2416215/1‑150 (MQ=255)
               tCAACGCTCCCGTGGTCACACCCCACCCCCCCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCCCTCGATTTCTTCCGGCTGAATCCCCTCTCCGCCACTg                                                                                                                                    <  2:80786/150‑1 (MQ=255)
                 aaCGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTgaa                                                                                                                                  >  2:1984470/1‑150 (MQ=255)
                   cGCTCTCGTGGTCATACCCCATCCCCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGCCGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTgaaga                                                                                                                                <  2:3146384/150‑1 (MQ=255)
                    gCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTgaagaa                                                                                                                               <  1:3759944/150‑1 (MQ=255)
                     ctctCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAAt                                                                                                                              <  1:4039070/150‑1 (MQ=255)
                     ctctCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAAt                                                                                                                              >  2:4735729/1‑150 (MQ=255)
                     ctctCGTGGTCATACCCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAAt                                                                                                                              <  2:676328/150‑1 (MQ=255)
                      tctcGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATa                                                                                                                             <  1:3688734/150‑1 (MQ=255)
                       ctcGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAg                                                                                                                            >  2:3828473/1‑150 (MQ=255)
                            ggTCATACCCCATCCCCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGt                                                                                                                       <  2:324193/150‑1 (MQ=255)
                              tCATACTCCACCACCGCCACGGGTCGATGACAAAACTCCTTGCCGCCACCGGGGACGATAACCACCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGCCTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCt                                                                                                                     <  2:86000/150‑1 (MQ=255)
                               cATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTa                                                                                                                    >  1:2057259/1‑150 (MQ=255)
                                 tACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTaaa                                                                                                                  <  1:4195482/150‑1 (MQ=255)
                                 tACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTaaa                                                                                                                  <  1:975475/150‑1 (MQ=255)
                                  ccTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAAc                                                                                                                 <  2:2915792/149‑1 (MQ=255)
                                  ccTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAAc                                                                                                                 <  1:3876668/149‑1 (MQ=255)
                                    tCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCGAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGt                                                                                                               >  1:1376903/1‑150 (MQ=255)
                                    tCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGt                                                                                                               >  2:1888842/1‑150 (MQ=255)
                                    cccATCACCGCCCCGGGTCGATGACCAAACTCCTTGTCGGCAACGGGGAAGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGt                                                                                                               <  1:2763541/149‑1 (MQ=255)
                                          aCCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACa                                                                                                         >  1:858679/1‑150 (MQ=255)
                                           ccGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACaa                                                                                                        <  1:358790/150‑1 (MQ=255)
                                             gcaaCCGGCAGATAACCAAACTCCTTGCCGGAACCGGAGACAATAAACTCCTGCAAACCCGAAGCATGTGCAGCAATTACGCGCTCGTCTTCTACCGGCAAAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATg                                                                                                      <  2:2944148/147‑1 (MQ=255)
                                             gCCACCGGTCGATGACCAACCTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGTCTTCTCCCGGCTGAATACCCACTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATg                                                                                                      <  2:1625846/150‑1 (MQ=255)
                                             gCCACCGGTCGATGACCAACCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATg                                                                                                      <  1:3901031/150‑1 (MQ=255)
                                              ccACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATCTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTACACGCCCGACAATgg                                                                                                     >  2:2599273/1‑150 (MQ=255)
                                                    gacGATGACCAACCTCTTTGTCGGCAACGGGGACGATAAACACCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGTCTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCt                                                                                               <  2:3481758/148‑1 (MQ=255)
                                                    gTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCt                                                                                               <  1:1693151/150‑1 (MQ=255)
                                                      cGATGACCAACCTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTg                                                                                             <  2:1697975/150‑1 (MQ=255)
                                                       gATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGc                                                                                            >  2:2153721/1‑150 (MQ=255)
                                                          gACCAAACTCCTTGTCGGCAGCGGGGACAATAAACCCTGGCAAACCCGGAGGATGCGCAGCAATTCCGCGGCGGACTTCTCCCGGCTGAAACCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCAt                                                                                         <  2:4310631/150‑1 (MQ=255)
                                                          gACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCACCAATTACGCGCTCGTCTTTTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCAt                                                                                         <  2:1342734/150‑1 (MQ=255)
                                                           aCCAAACTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAACCCGCAGGATGTGCAGCAATTCCGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCAtt                                                                                        <  2:3606585/150‑1 (MQ=255)
                                                            ccAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTa                                                                                       <  1:787807/150‑1 (MQ=255)
                                                            ccAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTa                                                                                       <  1:4559384/150‑1 (MQ=255)
                                                                   ccTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtct                                                                                >  1:2907581/1‑150 (MQ=255)
                                                                             aaCGGGGACGATAACCCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATc                                                                      <  2:2868985/150‑1 (MQ=255)
                                                                             aaCGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATGGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATc                                                                      >  1:4449552/1‑150 (MQ=255)
                                                                             aaCGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATc                                                                      >  1:1105689/1‑150 (MQ=255)
                                                                                ggggACGAGAAACACCTGCAAAACCGCAGGATGGGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACGGCGATcggg                                                                   >  2:2761265/1‑148 (MQ=255)
                                                                                   gACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgt                                                                <  2:212649/150‑1 (MQ=255)
                                                                                    aCGATAAACACCTGCAAAACCGGAGGATGTGCAGCAATTCCGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgta                                                               <  2:4054908/150‑1 (MQ=255)
                                                                                      gaaaaaCACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGCTTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtagc                                                             <  2:1121950/147‑1 (MQ=255)
                                                                                      gATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGTTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtagc                                                             >  2:3235945/1‑150 (MQ=255)
                                                                                      gATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtagc                                                             >  2:2408340/1‑150 (MQ=255)
                                                                                       aTAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtagcg                                                            >  2:411730/1‑150 (MQ=255)
                                                                                         aaaCCCCTGAAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtagcgta                                                          <  2:3760498/150‑1 (MQ=255)
                                                                                         aaaCACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtagcgta                                                          >  1:3165476/1‑150 (MQ=255)
                                                                                           acacCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTAcc                                                        <  1:1870151/150‑1 (MQ=255)
                                                                                             cccTGCAAACCCGCAGGATGTGCAGCAATTCCGCGCTCGCCTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAg                                                      <  2:1546859/149‑1 (MQ=255)
                                                                                                 gCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCtt                                                  >  2:106784/1‑150 (MQ=255)
                                                                                                   aaaaCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCg                                                <  1:1480964/150‑1 (MQ=255)
                                                                                                        cGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCAtt                                           >  2:584950/1‑150 (MQ=255)
                                                                                                         gCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTa                                          <  1:4733787/150‑1 (MQ=255)
                                                                                                         gCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTa                                          >  1:1075807/1‑150 (MQ=255)
                                                                                                               tgtgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGt                                    >  1:4679389/1‑150 (MQ=255)
                                                                                                                 tgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGa                                  >  1:832571/1‑150 (MQ=255)
                                                                                                                  gcagcaATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGaa                                 >  2:3923306/1‑150 (MQ=255)
                                                                                                                   cagcaATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGaaa                                >  1:3984126/1‑150 (MQ=255)
                                                                                                                   caccaATTACGCCGCCGACTTCTTCCGGCTGAAAACCCTATCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGaaa                                <  2:2968571/147‑1 (MQ=255)
                                                                                                                     gcaATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAAcc                              >  1:4374849/1‑150 (MQ=255)
                                                                                                                      caATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCa                             <  1:2653660/150‑1 (MQ=255)
                                                                                                                      caATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCa                             <  1:3245120/150‑1 (MQ=255)
                                                                                                                       aaTTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAg                            <  1:4292244/150‑1 (MQ=255)
                                                                                                                       aaTTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAg                            <  1:1843514/150‑1 (MQ=255)
                                                                                                                       aaTTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAg                            >  1:164109/1‑150 (MQ=255)
                                                                                                                        aTTCCGCGCTCGACTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGt                           <  2:4366096/150‑1 (MQ=255)
                                                                                                                        aTTACGCGCTCGATTTCTCCGGGCTGAATCCCCTCGCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGt                           <  1:3267666/150‑1 (MQ=255)
                                                                                                                         ttACGCGCTGGACTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGtt                          <  2:2037197/150‑1 (MQ=255)
                                                                                                                         ttACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGtt                          >  2:1264349/1‑150 (MQ=255)
                                                                                                                          tACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTg                         <  1:3371863/150‑1 (MQ=255)
                                                                                                                           aCGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGc                        >  2:2202632/1‑150 (MQ=255)
                                                                                                                               gcTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGCCTGTCTAAACGTCCGACAATGGTCCGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGt                    >  2:1283716/1‑150 (MQ=255)
                                                                                                                               gcTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGt                    >  1:73776/1‑150 (MQ=255)
                                                                                                                                   gccttcttCCGGCGAAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACg                <  1:2166541/148‑1 (MQ=255)
                                                                                                                                   gccttctcCCGGCTGAATCCCCTCTCCGCCACTGAAGAATATATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACg                <  2:1959207/148‑1 (MQ=255)
                                                                                                                                   gACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACg                <  1:4651813/150‑1 (MQ=255)
                                                                                                                                      ttctccCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCa             <  2:1206759/150‑1 (MQ=255)
                                                                                                                                       tcttcCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAg            <  1:5657/150‑1 (MQ=255)
                                                                                                                                       tcttcCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAg            <  1:3700168/150‑1 (MQ=255)
                                                                                                                                        cttcCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGt           <  1:2771132/150‑1 (MQ=255)
                                                                                                                                         ttcCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTa          >  1:1226312/1‑150 (MQ=255)
                                                                                                                                          tcCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTaa         >  2:1771586/1‑150 (MQ=255)
                                                                                                                                             ggCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCt      >  2:3856287/1‑150 (MQ=255)
                                                                                                                                              gCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCtt     <  2:2630674/150‑1 (MQ=255)
                                                                                                                                                taaaTCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCt   <  2:771090/148‑1 (MQ=255)
                                                                                                                                                 gAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCTg  <  2:2716628/150‑1 (MQ=255)
                                                                                                                                                 gAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTa                  <  2:3778206/134‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ACCCATTCACTAAGATCAACGCTCTCGTGGTCATACTCCATCACCGCCACCGGTCGATGACCAAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGCGTACCAGCCTTCGTCATTAACCAGTGAAACCAGTTGCCCGTTACGCCAGTAACCTTCTG  >  AP012306/1944139‑1944433

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