Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP012306472,8680TG22.5% 547.6 / 11.1 200D435A (GAC→GCC) cusSsensory histidine kinase in two‑component regulatory system with CusR, senses copper ions
Reads supporting (aligned to +/- strand):  ref base T (39/116);  new base G (45/0);  total (84/116)
Fisher's exact test for biased strand distribution p-value = 9.87e-22
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 255 total aligned reads displayed.

TGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGTT  >  AP012306/472719‑473016
                                                                                                                                                      |                                                                                                                                                     
  cTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCa                                                                                                                                                     <  1:3686046/150‑1 (MQ=255)
   tGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCaa                                                                                                                                                    <  1:326121/150‑1 (MQ=255)
    gAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAAc                                                                                                                                                   <  1:246036/150‑1 (MQ=255)
     aaGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGCTTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAAcg                                                                                                                                                  <  2:368385/150‑1 (MQ=255)
     aaGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAAcg                                                                                                                                                  >  1:1602251/1‑150 (MQ=255)
      aGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAAcgc                                                                                                                                                 >  2:2353588/1‑150 (MQ=255)
      aGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAAcgc                                                                                                                                                 <  2:906494/150‑1 (MQ=255)
        gATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGa                                                                                                                                               >  2:4153417/1‑150 (MQ=255)
           tAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAg                                                                                                                                            >  2:1806227/1‑150 (MQ=255)
            aaGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGTCCACGCGATAGa                                                                                                                                           >  2:1506152/1‑150 (MQ=255)
            aaGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGa                                                                                                                                           <  1:3324999/150‑1 (MQ=255)
             aGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGGTaaaa                                                                                                                                          >  2:2044794/1‑150 (MQ=255)
             aGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGaa                                                                                                                                          <  1:934041/150‑1 (MQ=255)
              gCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGaaa                                                                                                                                         <  2:4078413/150‑1 (MQ=255)
               cGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAAc                                                                                                                                        >  2:2430798/1‑150 (MQ=255)
                gggTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACCCGATAGAAACg                                                                                                                                       >  1:3872950/1‑150 (MQ=255)
                gggTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACg                                                                                                                                       >  1:1145926/1‑150 (MQ=255)
                  gTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACCCGATAAAAACggc                                                                                                                                     >  2:3455728/1‑149 (MQ=255)
                   tAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTc                                                                                                                                    <  2:3997696/150‑1 (MQ=255)
                   tAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTc                                                                                                                                    >  1:2969167/1‑150 (MQ=255)
                    aaTGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCa                                                                                                                                   >  2:837547/1‑150 (MQ=255)
                        tgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACa                                                                                                                               <  1:1989621/150‑1 (MQ=255)
                           tAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATc                                                                                                                            <  1:494413/150‑1 (MQ=255)
                           tAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATc                                                                                                                            <  1:3846842/150‑1 (MQ=255)
                            aaCAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCTCTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATcg                                                                                                                           >  1:2427177/1‑150 (MQ=255)
                            aaCAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATcg                                                                                                                           <  2:869234/150‑1 (MQ=255)
                               aaaCCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGTCCACGCGGTAGAAACGGGCAAACAATCGCgg                                                                                                                        >  1:3115251/1‑150 (MQ=255)
                                 aCCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTa                                                                                                                      <  1:1347231/150‑1 (MQ=255)
                                  ccTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAAAAAACGCGGTaa                                                                                                                     >  2:1178441/1‑150 (MQ=255)
                                    ttGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGt                                                                                                                   >  1:2883707/1‑150 (MQ=255)
                                     tGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGTCAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTg                                                                                                                  >  1:278126/1‑150 (MQ=255)
                                        cccccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGGCACACAATCGCGGTAAGTGCTc                                                                                                               >  1:3730420/1‑150 (MQ=255)
                                         ccccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGGCAAACAAACGCGGGAAGTGCTCg                                                                                                              >  2:1874431/1‑150 (MQ=255)
                                           ccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCggg                                                                                                            <  1:305907/150‑1 (MQ=255)
                                            cgcgCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGATGTGCTCGGgc                                                                                                           >  1:4200411/1‑150 (MQ=255)
                                            cgcgCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGgc                                                                                                           >  2:1515764/1‑150 (MQ=255)
                                            cgcgCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGgc                                                                                                           <  2:2403388/150‑1 (MQ=255)
                                               gcATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCa                                                                                                        <  1:2166876/150‑1 (MQ=255)
                                               gcATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCCGGCGCa                                                                                                        >  1:2279423/1‑150 (MQ=255)
                                                  tCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATg                                                                                                     <  1:2207042/150‑1 (MQ=255)
                                                    cGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATggg                                                                                                   >  1:772036/1‑150 (MQ=255)
                                                     gACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGc                                                                                                  >  1:172360/1‑150 (MQ=255)
                                                       cGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCCACGCGGTAGAAACGGCACAAAAAACTCGGTAAGAGCTTGGGCCCAATGGGGGt                                                                                                >  2:2975571/1‑150 (MQ=255)
                                                        gTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCgga                                                                                               >  2:1137706/1‑148 (MQ=255)
                                                         ttACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTAc                                                                                             <  1:3544404/150‑1 (MQ=255)
                                                           aCCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAAACGCGGGAAGTGCTCGGGCGCAATGGGCggaacg                                                                                            >  2:3106310/1‑145 (MQ=255)
                                                             cGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCggg                                                                                         >  1:1296209/1‑150 (MQ=255)
                                                               cAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGt                                                                                       <  1:2364167/150‑1 (MQ=255)
                                                                 aCCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGttt                                                                                     <  2:2473443/150‑1 (MQ=255)
                                                                    gTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTcg                                                                                  >  1:4116590/1‑150 (MQ=255)
                                                                      gCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTTAACGCGATAGAAACGGTCAAAAAATCGCGGTAAGTGGTTGGGGGCAATGGGGGTACCGGGGTGTTcgac                                                                                >  2:683769/1‑150 (MQ=255)
                                                                      gCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGTTCGGGCGCAATGGGCGTACCGGGGTTTTcgac                                                                                <  2:2375417/150‑1 (MQ=255)
                                                                      gCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTcgac                                                                                <  1:436735/150‑1 (MQ=255)
                                                                      gCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCCGGGGTTTcgac                                                                                >  1:1306550/1‑150 (MQ=255)
                                                                           tATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATaa                                                                           <  1:876298/150‑1 (MQ=255)
                                                                            aTGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCCAACCCGTAGAAAACGTCAAACAAAACCGGGAAGTGGTCGGGCGCCATGGGGGTACCGGGGGTTTTGACCATAAc                                                                          >  2:1390078/1‑150 (MQ=255)
                                                                                   acaacaATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCAcc                                                                   >  1:1510753/1‑150 (MQ=255)
                                                                                   acaacaATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCAcc                                                                   >  1:2668779/1‑150 (MQ=255)
                                                                                          tCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAAAGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGGGAt                                                            >  2:3864627/1‑150 (MQ=255)
                                                                                           cGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATc                                                           <  1:678731/150‑1 (MQ=255)
                                                                                              tttCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGTCAAAAAAACGCGGTAAGTGCTCCGGCGCCAAGGGCGGAACGGGGGTTTCGACGATAACTTGCACCCAGGGATCGAc                                                        >  2:4051701/1‑150 (MQ=255)
                                                                                               ttCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGAcc                                                       <  1:2730055/150‑1 (MQ=255)
                                                                                                tCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAAATTGCACCAGGTGATCGACCg                                                      >  1:1521484/1‑150 (MQ=255)
                                                                                                  aCTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTc                                                    <  2:19636/150‑1 (MQ=255)
                                                                                                  aCTATCGCCAGGCCAATACCGCTACCTTCACCTTCTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCCATGGGCGGACCGGGGGTTTCGACGATAACTTGCACCAGGTGGTCGACCGTc                                                    >  2:126538/1‑150 (MQ=255)
                                                                                                    tATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAAAGGGCGTACCGGGGGTTTCGACGATAACTTTCACCAGGTGATCGACCCTCTg                                                  >  2:259523/1‑150 (MQ=255)
                                                                                                         ccAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGGCAAACAAACGCGGGAAATGCTCGGGGGCAATGGGGGGACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGACCGTCTGACAgc                                             >  2:961202/1‑150 (MQ=255)
                                                                                                            ggCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCa                                          >  1:3621793/1‑150 (MQ=255)
                                                                                                            ggCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGACCGTCTGACAGCGCa                                          >  2:3897532/1‑150 (MQ=255)
                                                                                                             gCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAAAAAAGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGGACCGGGGTTTTTGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCAc                                         >  2:1450173/1‑150 (MQ=255)
                                                                                                             gCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCAc                                         <  1:392923/150‑1 (MQ=255)
                                                                                                              ccAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGACCGTCTGACAGCGCACt                                        >  1:1274108/1‑150 (MQ=255)
                                                                                                              ccAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGGTAAAAAAGGGCAAACAAACGCGGGAAGTGGTCGGGCGCAAAGGGGGGACCGGGGGTTTCGACGATAACATTTACCCGGGGATCGACCGTCTGGCCACGCACt                                        >  2:3187147/1‑150 (MQ=255)
                                                                                                               cAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTa                                       <  1:3988510/150‑1 (MQ=255)
                                                                                                               cAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTa                                       >  2:2248160/1‑150 (MQ=255)
                                                                                                                aaTACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCCGGCGCAATGGGCGGACCGGGGGTTTCGACGGAAACTTGCACCCGGTGATCGACCGTCTGACAGCGCACTAc                                      >  2:2500450/1‑150 (MQ=255)
                                                                                                                  tACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACaa                                    >  1:1592460/1‑150 (MQ=255)
                                                                                                                    ccGCTACCTTTACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAAtt                                  <  2:286286/150‑1 (MQ=255)
                                                                                                                    ccGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGGACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAACGCACTACAAtt                                  >  1:328263/1‑150 (MQ=255)
                                                                                                                     cGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTg                                 <  1:1360902/150‑1 (MQ=255)
                                                                                                                     cGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTg                                 >  1:1838799/1‑150 (MQ=255)
                                                                                                                      gCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAATGGTCGGGCGCCAAGGGGGTAACCGGGGTTTTGGCGGTAACTTTCCCCCGGGGGTCGCCCCTCCTACAACGCACCACaaattt                                >  2:4272834/1‑146 (MQ=255)
                                                                                                                      gCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGt                                <  2:722708/150‑1 (MQ=255)
                                                                                                                        tACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGGTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGtct                              >  2:497189/1‑150 (MQ=255)
                                                                                                                          ccTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCCCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGACCGCCTGTCAGGGCACTACAATTGtctct                            >  2:3273429/1‑150 (MQ=255)
                                                                                                                           cTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGGGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAAATTGCACCAGGGGATCGACCGTCTGACAGCGCCCTACAATTGtctctc                           >  2:2083471/1‑150 (MQ=255)
                                                                                                                           cTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGtctctc                           <  2:2727379/150‑1 (MQ=255)
                                                                                                                              cACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCgg                        <  1:3637054/150‑1 (MQ=255)
                                                                                                                               aCCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAAAGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTTGACGATAACTTGCACCAGGGGATCGACCGTCTGACAGCGCACTACAAATGTCTCTCCCGt                       >  2:1230138/1‑150 (MQ=255)
                                                                                                                                ccTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGtt                      <  2:3333778/150‑1 (MQ=255)
                                                                                                                                ccTTTTCGCTGGCGGCAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGtt                      <  2:3508679/150‑1 (MQ=255)
                                                                                                                                      cGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGta                >  2:2190502/1‑150 (MQ=255)
                                                                                                                                        cTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGtata              <  1:2624291/150‑1 (MQ=255)
                                                                                                                                         tGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCTATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATaa             <  2:3542494/150‑1 (MQ=255)
                                                                                                                                         tGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATaa             >  1:2569136/1‑150 (MQ=255)
                                                                                                                                            cGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGc          <  1:18764/150‑1 (MQ=255)
                                                                                                                                             gggAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCa         >  2:3773047/1‑150 (MQ=255)
                                                                                                                                                aGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAggg      >  2:916783/1‑150 (MQ=255)
                                                                                                                                                aGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAggg      <  2:3043146/150‑1 (MQ=255)
                                                                                                                                                   gggTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCg    >  2:4018711/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGTT  >  AP012306/472719‑473016

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