Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP012306472,8680TG24.2% 562.8 / 19.7 211D435A (GAC→GCC) cusSsensory histidine kinase in two‑component regulatory system with CusR, senses copper ions
Reads supporting (aligned to +/- strand):  ref base T (34/126);  new base G (51/0);  total (85/126)
Fisher's exact test for biased strand distribution p-value = 2.01e-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 282 total aligned reads displayed.

TGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGTTC  >  AP012306/472719‑473017
                                                                                                                                                     |                                                                                                                                                      
   tGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCaa                                                                                                                                                     <  2:296242/150‑1 (MQ=255)
     aaGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAAcg                                                                                                                                                   <  1:612514/150‑1 (MQ=255)
     aaGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAAcg                                                                                                                                                   >  1:4041100/1‑150 (MQ=255)
        gATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGa                                                                                                                                                >  2:1668286/1‑150 (MQ=255)
          ttAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATa                                                                                                                                              <  1:2515434/150‑1 (MQ=255)
          ttAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATa                                                                                                                                              >  1:836849/1‑150 (MQ=255)
           tAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGACCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGAGGGGCAACGCGagga                                                                                                                                             >  2:1000385/1‑147 (MQ=255)
              gCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGaaa                                                                                                                                          <  2:1876768/150‑1 (MQ=255)
              gCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGaaa                                                                                                                                          <  1:1151440/150‑1 (MQ=255)
               cGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAAc                                                                                                                                         >  1:2106330/1‑150 (MQ=255)
                gggTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCACAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACg                                                                                                                                        >  1:2538954/1‑150 (MQ=255)
                  gTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGt                                                                                                                                      >  1:2480375/1‑150 (MQ=255)
                    aaTGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCa                                                                                                                                    <  1:3503211/150‑1 (MQ=255)
                     aTGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCaa                                                                                                                                   >  1:2405282/1‑150 (MQ=255)
                      tgtgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCaaa                                                                                                                                  <  1:2840324/150‑1 (MQ=255)
                      tgtgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCaaa                                                                                                                                  <  2:3731321/150‑1 (MQ=255)
                      tgtgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAAAGGTCaaa                                                                                                                                  >  1:224936/1‑150 (MQ=255)
                       gtgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGTGGGGGGAACGCGTTAAAAGCAGGCAAAc                                                                                                                                 >  2:1576701/1‑150 (MQ=255)
                       gtgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAAc                                                                                                                                 <  2:1484405/150‑1 (MQ=255)
                       gtgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAAAGGTCAAAc                                                                                                                                 >  2:1745310/1‑150 (MQ=255)
                        tgATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACCCGATAGAAACAGTcacaca                                                                                                                                >  1:3219755/1‑150 (MQ=255)
                           tAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATc                                                                                                                             <  1:1538268/150‑1 (MQ=255)
                             aCAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCCGTAGAAACCGTCcaacaaaagc                                                                                                                           >  2:1394972/1‑146 (MQ=255)
                               aaaCCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAAAGGGCCAACAAACCCgg                                                                                                                         >  2:2226677/1‑150 (MQ=255)
                                   cTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGTCAAACAATCCCGGTAAg                                                                                                                    >  1:673829/1‑150 (MQ=255)
                                   cTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAAAGCGGTAGAAACCGGCAAACAATAGCGGttat                                                                                                                    >  1:2919647/1‑147 (MQ=255)
                                    ttGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGt                                                                                                                   <  1:2157631/150‑1 (MQ=255)
                                      gTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGGTAGAAAAGGGCAAAAAAACGCGGGAAGGGc                                                                                                                 >  2:2685417/1‑150 (MQ=255)
                                       tCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCt                                                                                                                >  1:1599192/1‑150 (MQ=255)
                                         ccccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCg                                                                                                              <  2:515889/150‑1 (MQ=255)
                                          cccGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCCATAAAAACAGTCAAAAAAAAGCGGCAAGGGCTCgg                                                                                                             >  2:3270969/1‑150 (MQ=255)
                                             gcgcATCCGACGTTACCGAAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATTCCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGgcg                                                                                                          <  1:1014962/150‑1 (MQ=255)
                                               gcATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCa                                                                                                        <  1:1892429/150‑1 (MQ=255)
                                               gcATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCCACCCGATAGAAACGGGCAAACAATCGCGGGAAGTGGCCCGGCGCa                                                                                                        >  1:3192810/1‑150 (MQ=255)
                                               gcATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACACGAAAGAAACGGTCAAACAAACGCGGTAAGTGGCCGGGCGCa                                                                                                        >  2:261265/1‑150 (MQ=255)
                                                cATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCaa                                                                                                       >  1:1991577/1‑150 (MQ=255)
                                                 aTCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAAt                                                                                                      <  1:3596981/150‑1 (MQ=255)
                                                  tCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGGCAAACAACCGCGGGAAGTGCACCGGCGCAATg                                                                                                     >  1:3451868/1‑150 (MQ=255)
                                                   ccGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATgg                                                                                                    <  1:232014/150‑1 (MQ=255)
                                                    cGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCCACGCGATAGAAACGGTCCAACAATCGCGGGAAGTGCTCGGGCGCAATggg                                                                                                   >  2:3200469/1‑150 (MQ=255)
                                                     gACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAAAAACCGTCCAACAATCGCCGGAAGTGCCCCGGGCCAATGGGc                                                                                                  >  1:3931642/1‑150 (MQ=255)
                                                      aCGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGCAGGGTCAACCCGATAAAAACGGGCAAACAATCGCGGGTAGTGACCGGGCGCAACggggg                                                                                                 >  2:2337558/1‑148 (MQ=255)
                                                        gTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTa                                                                                               >  1:895229/1‑150 (MQ=255)
                                                            ccGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGGACCgg                                                                                           >  2:1267047/1‑150 (MQ=255)
                                                                aaCCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGtt                                                                                       <  1:1356188/150‑1 (MQ=255)
                                                                aaCCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGGTAGAAAACGGCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCgggggt                                                                                       >  2:4081338/1‑148 (MQ=255)
                                                                 aCCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGttt                                                                                      <  1:4011032/150‑1 (MQ=255)
                                                                 aCCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGGACCGGGGgtt                                                                                      >  1:2198540/1‑150 (MQ=255)
                                                                    gTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATTGGCGTACCGGGGGTTTcg                                                                                   >  1:365672/1‑150 (MQ=255)
                                                                      gCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAAAGGTCAAACAAACGCGGTAAGTGCTCGGGGGCAATGGGCGTACCGGGGGTTTtgac                                                                                 >  2:1987435/1‑150 (MQ=255)
                                                                       ccTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGGGCAATGGGCGTACCGGGGGTTTcgacg                                                                                >  2:3249129/1‑150 (MQ=255)
                                                                        cTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGTCAACGCGATAGAAACGGTCAAACAAACGCGGTAAGGGCTCGGGCGGAAAAGGGGGAACGGGGGTTTtgacga                                                                               >  2:1920509/1‑150 (MQ=255)
                                                                        cTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGGGCTCGGGCGCAATGGGCGTACCCGGGGTTTtgacga                                                                               >  1:1646896/1‑150 (MQ=255)
                                                                           tATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATaa                                                                            <  1:4147626/150‑1 (MQ=255)
                                                                           tATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATaa                                                                            <  1:2906118/150‑1 (MQ=255)
                                                                             tGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGGTAGAAACGGGCAAACAATCGCGGTAAGTGCTCGGGCGCCATGGGCGTTCCCGGGGTTTTGACGATaaat                                                                          >  2:2504087/1‑148 (MQ=255)
                                                                              gcgcGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACtt                                                                         <  2:1116396/150‑1 (MQ=255)
                                                                                      acaaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAgg                                                                 <  1:1294086/150‑1 (MQ=255)
                                                                                          tCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGAt                                                             <  2:1674078/150‑1 (MQ=255)
                                                                                           cGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATc                                                            <  1:1487075/150‑1 (MQ=255)
                                                                                           cGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATc                                                            <  1:2968490/150‑1 (MQ=255)
                                                                                            gATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGGAAAGGGGGAGAAAAGGGGAAAAAAAAGGGGGAAGGGGGGGGGGGGGATGGGGGGAGCGGGGGTTTTGACGAGAAATTGGACCAgggggagg                                                           >  2:425285/1‑144 (MQ=255)
                                                                                             aTTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGCGGGGTCAACGCGATAAAAACGGGTAAACAAACGCGGCAAAGTGCCGGGGCGCATGGGCGTTCCCGGGATGTTGAGGCGAACCTTCACCAGGAGGTCGa                                                          >  2:2768169/1‑150 (MQ=255)
                                                                                             aTTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGa                                                          >  1:2805645/1‑150 (MQ=255)
                                                                                                   cTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCt                                                    <  2:902566/150‑1 (MQ=255)
                                                                                                    tATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGGTAGAAAAGGTCAAACAATTGCGGTAAGTGCTCGGGCGCAATTGGCGTAACGGGGGTTTCGACGATAACTTGCACCAGGTGGTCGACCGTCTg                                                   >  2:2284158/1‑150 (MQ=255)
                                                                                                    tATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGTCAAAAAATCCCGGTAAATGCTCGGGCGCAAAGGGGGGACCGGGGGTTTTGACGGTAACTTTCACCAAGGGATCGACCGTCTg                                                   >  2:3542071/1‑150 (MQ=255)
                                                                                                       cGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACa                                                <  1:3567590/150‑1 (MQ=255)
                                                                                                        gCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAg                                               >  1:1077177/1‑150 (MQ=255)
                                                                                                           aGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCCATGGGCGTACCCGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAgcgc                                            >  1:3764410/1‑150 (MQ=255)
                                                                                                           aGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGACCGTCTGACAgcgc                                            >  1:3041360/1‑150 (MQ=255)
                                                                                                               cAATACCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGGGTACCGGGGGTTTCGACGATAACTTGCACCAGGGGATCGACCGTCTGACAGCGCACTa                                        >  2:477982/1‑150 (MQ=255)
                                                                                                               cAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTa                                        <  2:657663/150‑1 (MQ=255)
                                                                                                                aaTACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCCACGCGATAGAAACGGGCAAACAATCGCGGGAAGTGGTCGGGCGCAATGGGGGTACCGGGGGTTTTGACGATAAATTTCACCAGGGGATCGACCCTCTGACAGCGCCCTAc                                       >  2:3600866/1‑150 (MQ=255)
                                                                                                                aaTACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTAc                                       <  1:1067275/150‑1 (MQ=255)
                                                                                                                aaTACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGGCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTAc                                       >  1:261829/1‑150 (MQ=255)
                                                                                                                aaTACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTAc                                       >  2:580650/1‑150 (MQ=255)
                                                                                                                  tACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACaa                                     <  2:145581/150‑1 (MQ=255)
                                                                                                                   aCCGCTACCTTCACCTTTTCGCTGGCGGGGGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTTGACGATAACTTGCACCAGGTGGTCGACCGTCTGACAGCGCACTACAAt                                    >  2:3447037/1‑150 (MQ=255)
                                                                                                                    ccGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGGCAAACAATCACGGTAAAAGGCCGGGGGCACTGGGCGTACCGGGGGTTTCGACGATAAATTGCCCCAGGGGATCGACCGTCTTGCAGCGCACTCCAAtt                                   >  2:1962548/1‑150 (MQ=255)
                                                                                                                       cTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGGTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGtc                                >  1:3881352/1‑150 (MQ=255)
                                                                                                                          ccTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGtctct                             <  1:3793982/150‑1 (MQ=255)
                                                                                                                           cTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGtctctc                            <  2:2854395/150‑1 (MQ=255)
                                                                                                                           cTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGtctctc                            <  1:4108274/150‑1 (MQ=255)
                                                                                                                            ttCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAAAGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTAACGGGGGTTTCGACGATAACTTGCACCAGGGGATCGACCGTCTGACAGCGCACTACAATTGTCTCTcc                           >  2:3860934/1‑150 (MQ=255)
                                                                                                                             tCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCg                          <  2:3796846/150‑1 (MQ=255)
                                                                                                                                 cTTTTCGCTGGCGGGAGGGGGCAACGCGGGAAAAACGGTGAAAAAATCGCCGGAAGGGGCCCGGTGCTAGGGGCCTGCCGGGGTTTTTGGCGATAACTTGCCCCGGGTGAGCGCCCGTCCGAAAGGCCACTACAATTTTACCTCCgcatg                      >  2:21573/1‑146 (MQ=255)
                                                                                                                                  ttttCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTgg                     <  2:1045456/150‑1 (MQ=255)
                                                                                                                                   tttCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGc                    >  1:3349219/1‑150 (MQ=255)
                                                                                                                                   tttCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGc                    <  2:2378224/150‑1 (MQ=255)
                                                                                                                                      cGCTGGCGGGAGGGGGCAACGCGATAGGAACGGGCAAACAATCCCGGTAAGTGCTCGGGCGCGATGGGCGGGCCGGGGTTTTCGACGATAACTTGCACCAGGTGATAGACCGTCTGACAGCGCACTCCAATTTTCTCTCCAGTTGGCGta                 >  2:1150566/1‑150 (MQ=255)
                                                                                                                                       gCTGGCGGGAGGTGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGtat                <  2:2077552/150‑1 (MQ=255)
                                                                                                                                       gCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGtat                <  2:3116830/150‑1 (MQ=255)
                                                                                                                                        cTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGGAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGtata               >  2:1245112/1‑150 (MQ=255)
                                                                                                                                          ggCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAAc             <  1:1207246/150‑1 (MQ=255)
                                                                                                                                             gggAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCa          >  1:1799228/1‑150 (MQ=255)
                                                                                                                                              ggAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAg         <  1:2076201/150‑1 (MQ=255)
                                                                                                                                                aGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAggg       >  2:3032243/1‑150 (MQ=255)
                                                                                                                                                  gggTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCg     <  2:3318819/150‑1 (MQ=255)
                                                                                                                                                    gTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGGTTGGCGTATAACGCAGGGCGtt   <  2:756597/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 < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: