Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP012306472,8680TG23.7% 568.2 / 16.2 220D435A (GAC→GCC) cusSsensory histidine kinase in two‑component regulatory system with CusR, senses copper ions
Reads supporting (aligned to +/- strand):  ref base T (38/129);  new base G (52/0);  total (91/129)
Fisher's exact test for biased strand distribution p-value = 4.11e-26
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 320 total aligned reads displayed.

TGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGTTC  >  AP012306/472719‑473017
                                                                                                                                                     |                                                                                                                                                     
  cTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCACCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCa                                                                                                                                                     <  1:3211710/150‑1 (MQ=255)
     aaGGATTAAGCGGGTAATGTGATAACAAACCTTGACCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAAcg                                                                                                                                                  <  1:1201225/150‑1 (MQ=255)
      aGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAAcgc                                                                                                                                                 <  2:2360028/150‑1 (MQ=255)
      aGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAAcgc                                                                                                                                                 >  1:2429261/1‑150 (MQ=255)
          ttAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGaca                                                                                                                                             <  1:234431/150‑3 (MQ=255)
          ttAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATa                                                                                                                                             <  1:4241435/150‑1 (MQ=255)
           tAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAg                                                                                                                                            >  1:682323/1‑150 (MQ=255)
             aGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGaa                                                                                                                                          >  1:3613884/1‑150 (MQ=255)
              gCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGTCAACGCGATAGaaa                                                                                                                                         <  1:1861878/150‑1 (MQ=255)
              gCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGaaa                                                                                                                                         <  2:1631857/150‑1 (MQ=255)
               cGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAAc                                                                                                                                        <  1:3805392/150‑1 (MQ=255)
                gggTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACg                                                                                                                                       <  1:481790/150‑1 (MQ=255)
                  gTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGt                                                                                                                                     <  1:2169905/150‑1 (MQ=255)
                  gTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGt                                                                                                                                     >  1:3745370/1‑150 (MQ=255)
                   tAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGTc                                                                                                                                    >  2:3616274/1‑150 (MQ=255)
                   tAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGGCACGCGATAGCAAAGGtg                                                                                                                                    >  2:2606899/1‑149 (MQ=255)
                        tgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACa                                                                                                                               <  2:422567/150‑1 (MQ=255)
                         gATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACaa                                                                                                                              <  2:58452/150‑1 (MQ=255)
                              cAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATcgcg                                                                                                                         >  1:2594052/1‑150 (MQ=255)
                              cAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATcccg                                                                                                                         >  2:3573802/1‑150 (MQ=255)
                               aaaCCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCgg                                                                                                                        >  2:666256/1‑150 (MQ=255)
                               aaaCCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATTGCgg                                                                                                                        >  1:91705/1‑150 (MQ=255)
                               aaaCCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCgg                                                                                                                        >  2:3737932/1‑150 (MQ=255)
                                aaCCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGt                                                                                                                       >  1:4314762/1‑150 (MQ=255)
                                 aCCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTa                                                                                                                      <  1:4126545/150‑1 (MQ=255)
                                 aCCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTa                                                                                                                      <  1:3474514/150‑1 (MQ=255)
                                  ccTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTaa                                                                                                                     <  1:3277509/150‑1 (MQ=255)
                                   cTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAAGGGTCAACGCGATAGAAACGGGCAAACAATCGCGGTAAg                                                                                                                    >  1:3062579/1‑150 (MQ=255)
                                    ttGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGt                                                                                                                   <  1:2767963/150‑1 (MQ=255)
                                    ttGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGt                                                                                                                   <  2:3144646/150‑1 (MQ=255)
                                       tCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTTCACGCGGTAGAAAAGGTCAAACAATCCCGGGAAGTGCt                                                                                                                >  1:1010170/1‑150 (MQ=255)
                                        cccccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTc                                                                                                               <  2:2663054/150‑1 (MQ=255)
                                        cccccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTc                                                                                                               <  1:2036759/150‑1 (MQ=255)
                                           ccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTCCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGCGGGGGCAACGCGCTAGTAACCGTCAAAAAATCTCGGGCAGTGCTCtgg                                                                                                            >  1:1910493/1‑150 (MQ=255)
                                            cgcgCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGgc                                                                                                           <  1:1116301/150‑1 (MQ=255)
                                            cgcgCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGgc                                                                                                           <  1:1718695/150‑1 (MQ=255)
                                             gcgcATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATTCCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGgcg                                                                                                          <  2:433047/150‑1 (MQ=255)
                                               gcATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGGCAAACAATCGCGGGAAGTGCTCGGGCGCa                                                                                                        >  2:2740172/1‑150 (MQ=255)
                                                 aTCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAAt                                                                                                      <  2:9484/150‑1 (MQ=255)
                                                  tCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATg                                                                                                     <  2:4385260/150‑1 (MQ=255)
                                                   ccGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGGCAAACAATCGCGGTAAGGGCTCGGGCGCAATgg                                                                                                    >  1:281513/1‑150 (MQ=255)
                                                     gACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGc                                                                                                  <  1:2710919/150‑1 (MQ=255)
                                                         ttACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTAc                                                                                              >  1:1444114/1‑150 (MQ=255)
                                                               cAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGt                                                                                        <  1:1539730/150‑1 (MQ=255)
                                                               cAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAAAAACGGTCAAAAAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCcgggg                                                                                        >  2:449824/1‑149 (MQ=255)
                                                                  ccGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGtttt                                                                                     <  1:3436945/150‑1 (MQ=255)
                                                                  ccGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGGGCTCGGGCGCAATGGGGGTACCGGGGgttt                                                                                     >  2:4138453/1‑150 (MQ=255)
                                                                  ccGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACCCGATAGAAAAGGGCAAAAAATCTTGGTAAGTGCGCGGGGGCAAAGGGCGTACCcgggggtt                                                                                     >  2:2852133/1‑146 (MQ=255)
                                                                   cGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTc                                                                                    >  1:293000/1‑150 (MQ=255)
                                                                     tGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCCATGGGCGTACCGGGGGTTTcga                                                                                  >  1:802562/1‑150 (MQ=255)
                                                                              gcgcGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACtt                                                                         <  1:4064184/150‑1 (MQ=255)
                                                                                 cgACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCa                                                                      <  2:120394/150‑1 (MQ=255)
                                                                                  gACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCCATGGGCGTACCGGGGTTTTCGACGATAACTTTCAc                                                                     >  2:1121684/1‑150 (MQ=255)
                                                                                    caacaaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCa                                                                   >  1:1224253/1‑150 (MQ=255)
                                                                                     aacaaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAAATTGCACCAg                                                                  >  2:3911238/1‑150 (MQ=255)
                                                                                       caaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGTCCAACAATCGCGGTAAGGGGTCGGGCGCAAAGGGCGTAACGGGGGTTTTGTCGAAAAATTGCACCAGGt                                                                >  2:700859/1‑150 (MQ=255)
                                                                                       caaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAggg                                                                >  1:1241050/1‑149 (MQ=255)
                                                                                         aTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGa                                                              >  1:3408320/1‑150 (MQ=255)
                                                                                          tCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGAt                                                             <  1:3175804/150‑1 (MQ=255)
                                                                                           cGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATc                                                            <  1:2488588/150‑1 (MQ=255)
                                                                                             aTTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGTGGGAAGGGTACGGGCGCAATGGGCGTACCGGGGATTTCGACGGTAACTTGCACCAGGTGATCGa                                                          >  2:1825231/1‑150 (MQ=255)
                                                                                             aTTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGGGTACCGGGGGTTTTGACGATAACTTTCACCAGGTGATCGa                                                          >  2:531993/1‑150 (MQ=255)
                                                                                             aTTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGGACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGa                                                          >  2:2971383/1‑150 (MQ=255)
                                                                                                tCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCg                                                       <  2:913681/150‑1 (MQ=255)
                                                                                                 cACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGt                                                      <  1:1599174/150‑1 (MQ=255)
                                                                                                    tATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTg                                                   <  2:1721370/150‑1 (MQ=255)
                                                                                                       cGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACa                                                <  1:2236733/150‑1 (MQ=255)
                                                                                                        gCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAg                                               <  2:3232814/150‑1 (MQ=255)
                                                                                                          cAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGGACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGACCGGCTGCCAgcg                                             >  2:3382037/1‑150 (MQ=255)
                                                                                                          cAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAgcg                                             >  2:1540442/1‑150 (MQ=255)
                                                                                                             gCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGGGAAGCGATAGAAAAGGGGAAACAAAAGCGGGGAGGGGTCGGGGGCAATGGGGGGAACGGGGGTTTCGGCGATAACTTGGACCAGGGGGTTGACCGGCTGAAAggggac                                          >  2:2627914/1‑145 (MQ=255)
                                                                                                              ccAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACt                                         <  1:2395784/150‑1 (MQ=255)
                                                                                                              ccAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGGCAAACAATCGCGGTAAGTGCTCGGGGGCGATGGGGGGACCGGGGGTTTTGACGATAACTTGCACCAGGTGATCGACCGGCTGACAGCGCACt                                         >  2:1407157/1‑150 (MQ=255)
                                                                                                               cAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCCACGCGAAACAAAAGATCCAACAATCGCGGTAAGAGGTCCGGCGCACTCGGCGTACCGGGGGTTTCGACGACAACTTGCACCAGGTGAGCGACCGCCTGACAGCCCGCTa                                        >  2:1988697/1‑150 (MQ=255)
                                                                                                               cAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGACCGTCTGACAGCGCACTa                                        >  2:2880644/1‑150 (MQ=255)
                                                                                                               cAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAAAGGGCAAACAATCGCGGGAAGGGCTCCGGCGCAAAGGGCGTACCCGGGGTTTTGACGATAACTTGCACCAGGGGATCGACCGTCTGACAacgcccca                                        >  2:2051063/1‑146 (MQ=255)
                                                                                                                     cGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTg                                  <  2:2540212/150‑1 (MQ=255)
                                                                                                                      gCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAAACGGAAAGGAACGGGCAAACAATCGCGGTAAGCGGTCCGGCGCCATGGGCGGACCGGGGGTTTTGACGAAAAACTCGACCAAGGGATTGACCGTCTGACAGGGCACTTCCATTGt                                 >  2:2220443/1‑150 (MQ=255)
                                                                                                                      gCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGt                                 >  1:736813/1‑150 (MQ=255)
                                                                                                                      gCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGGGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGACCGTCTGACAGCGCCCTACAAtttt                                 >  2:3464424/1‑148 (MQ=255)
                                                                                                                      gCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACCGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGACCGTCTGACAGCGCCCTACAAtttt                                 >  2:246457/1‑148 (MQ=255)
                                                                                                                       cTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGtc                                >  1:762815/1‑150 (MQ=255)
                                                                                                                        tACCTTCACCTTTTCGCTGGCGGGGGGGGTCAACGCGGTAGAAAAGGTCAAACAATCGCGGGAAGGGCTCGGGGGCAATGGGCGTACCGGGGGTTTTGCCGATAACTTGCACCAGGGGATCGACCGTCTGACAGCGCCCCACAATTTtct                               >  2:3795054/1‑150 (MQ=255)
                                                                                                                         aCCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGtctc                              >  1:4167152/1‑150 (MQ=255)
                                                                                                                             tCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCg                          <  2:4198531/150‑1 (MQ=255)
                                                                                                                             tCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCg                          <  2:4259653/150‑1 (MQ=255)
                                                                                                                             tCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCCGACAGCGCACTACAAATGTCTCTCCg                          >  2:1746376/1‑150 (MQ=255)
                                                                                                                               aCCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGt                        >  1:880086/1‑150 (MQ=255)
                                                                                                                                 cTTTTCGCTGGCGGGAGGGGTCAACGCGGTAGAAACGGTCAAAAAATCGCGGTAAGGGCTCGGGGGCAATGGGCGGACCGGGGGTTTCGACGATAACTTGCACCAGGGGGTCGACCGTCTGACAGCGCACTACAATTGTTTCTCCGGTTg                      >  2:1688800/1‑150 (MQ=255)
                                                                                                                                   tttCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCCGTTGGc                    >  2:21550/1‑150 (MQ=255)
                                                                                                                                        cTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGtata               >  2:3865310/1‑150 (MQ=255)
                                                                                                                                           gCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACg            <  1:140389/150‑1 (MQ=255)
                                                                                                                                            cGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGATGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGc           <  1:3993301/150‑1 (MQ=255)
                                                                                                                                            cGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGc           >  2:1050169/1‑150 (MQ=255)
                                                                                                                                            cGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGAAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGc           >  2:784463/1‑150 (MQ=255)
                                                                                                                                               gAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAgg        <  2:1322186/150‑1 (MQ=255)
                                                                                                                                                aGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAggg       <  2:1512328/150‑1 (MQ=255)
                                                                                                                                                aGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAggg       >  2:1953652/1‑150 (MQ=255)
                                                                                                                                                 gggtTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGc      <  1:1304082/150‑1 (MQ=255)
                                                                                                                                                     tCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGTTc  <  2:3319688/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGTTC  >  AP012306/472719‑473017

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