Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,875,8610AC34.2% 110.0 / 10.2 73T298P (ACC→CCC) ttgI_2Toluene efflux pump outer membrane protein TtgI
Reads supporting (aligned to +/- strand):  ref base A (13/35);  new base C (25/0);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 2.41e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCC  >  ACB143/3875727‑3875988
                                                                                                                                      |                                                                                                                               
cTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGGGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCAc                                                                                                                                <  6:240240/136‑1 (MQ=255)
   aCGCGCTCAGCGTAGACCTGGCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCAcc                                                                                                                               <  2:300823/134‑1 (MQ=255)
     gcgcTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCCGCGGCCACCGccc                                                                                                                           >  3:122971/1‑135 (MQ=255)
          cAGCGTAGACCTGGCGCCCACGGCCCTGCCGGCCAGCGCCACGGCGCGGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGtt                                                                                                                      <  5:198643/136‑1 (MQ=255)
           aGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCTGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCAACGTTg                                                                                                                     >  1:236693/1‑136 (MQ=255)
           aGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCCACGTTg                                                                                                                     >  3:260050/1‑136 (MQ=255)
            gCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCAACGTTgg                                                                                                                    >  6:99862/1‑136 (MQ=255)
               tAGACCTGTCGCCCACGGCCCTGCCGGCCAGCGCCAAGGCGCTGCCGGGCGGTGCCCCGGGGGAGCTGCGGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGt                                                                                                                 <  8:221279/136‑1 (MQ=255)
                         gCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCAccgcc                                                                                                       <  6:93309/136‑1 (MQ=255)
                          cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCCCCGGTGGCGTGGCCAccgccg                                                                                                      >  5:311473/1‑136 (MQ=255)
                          cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCCACGGTGGCGTGGccccccccc                                                                                                      >  3:226669/1‑135 (MQ=255)
                          cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCCACCTTGGCGTTGCCAccgcc                                                                                                       >  6:126491/1‑135 (MQ=255)
                          cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCACCCTTGGCGTGGccccccccc                                                                                                      >  1:283505/1‑135 (MQ=255)
                          cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCAACGTTGGCGTGGCCCccgccg                                                                                                      >  7:253823/1‑136 (MQ=255)
                          cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCAACGTTGGCGTGGCCAccgccg                                                                                                      >  6:101244/1‑136 (MQ=255)
                          cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCAACGTTGGCGGGGCCAccgccg                                                                                                      >  5:248587/1‑136 (MQ=255)
                          cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCAACCTTGGCGTGGCCCccgccg                                                                                                      >  5:75672/1‑136 (MQ=255)
                          cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCCCCAACGTTGGCGTGGCCAccgccg                                                                                                      >  8:97553/1‑136 (MQ=255)
                          cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCCAAGTTGGCGTGGACCcccccg                                                                                                      >  8:303536/1‑136 (MQ=255)
                          cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCACCGTTGGCGTGGCCAccgccg                                                                                                      >  5:337226/1‑136 (MQ=255)
                          cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCACCCTTGGCGTGGCCAccgccg                                                                                                      >  3:286669/1‑136 (MQ=255)
                          cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCAccgccg                                                                                                      >  5:351268/1‑136 (MQ=255)
                          cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCAccgccg                                                                                                      >  6:356075/1‑136 (MQ=255)
                          cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACCTTGGCGTGGCCCccgccg                                                                                                      >  8:336323/1‑136 (MQ=255)
                               ggCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCCACGTTGGCGTGGCCCCCGCCGCCCTg                                                                                                 >  4:130823/1‑136 (MQ=255)
                                       cGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCCCCGATGGCGTGGCCCCCCCCGCCCGGGTCCcgcg                                                                                         >  4:33716/1‑136 (MQ=255)
                                       cGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCCACGTTGGCGTGGCCACCGCCGCCCGGGTCCcgcg                                                                                         >  5:68351/1‑136 (MQ=255)
                                       cGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCcgcg                                                                                         >  3:152850/1‑136 (MQ=255)
                                        ggCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCcgcgc                                                                                        <  6:170627/136‑1 (MQ=255)
                                        ggCCATCGCCAAGCCGCTGCCGGTCGGTGACCCCGGTGAGCTGCTCCGCCGCCGCCCGGACATCAGCAGCGCCGAGCGCCAGCTGGCGGCGTCCACCCCCAACCTAGGCCTGGCCACCGCCGCCCAGTGACcgcgc                                                                                        >  1:329377/1‑136 (MQ=255)
                                        ggCCAGCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCGGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCcgcgc                                                                                        <  3:112350/136‑1 (MQ=255)
                                              cGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCCCCGTTGGCGTGGCCCCCGCCGCCCGGTTCCCGCGCGTCCGc                                                                                  >  6:88344/1‑136 (MQ=255)
                                              cGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCAACGTTGGCGTGGCCACCCCCGCCCTGTTCCCGCGCGTccga                                                                                  >  2:234301/1‑133 (MQ=255)
                                               gCCAAGGCGGTGCGGGTGGGTGACCGGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGcc                                                                                 <  6:309217/136‑1 (MQ=255)
                                               gCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGcc                                                                                 <  7:118834/136‑1 (MQ=255)
                                                       gcTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCCACGTTGGCGTGGCCCCCGCCGCCCGGTTCCCCCGCGTCCCCCGCCGTGGc                                                                         >  8:135494/1‑136 (MQ=255)
                                                       gcTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCAACGTTGGCGGGGCCCCCCCCGCCCTGGTGCCGCGCGTCAGCCGCAggggg                                                                         >  3:116957/1‑135 (MQ=255)
                                                           ccGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGTACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCtt                                                                       <  8:20258/134‑1 (MQ=255)
                                                           ccGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCCACGTTGGCGTGGCCACCGCCGACCTGGTCCCGCGCGTACGCCTCAGTGGCTTcc                                                                     >  7:255585/1‑136 (MQ=255)
                                                           ccGGTCGGGGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTc                                                                      <  5:111977/135‑1 (MQ=255)
                                                                  gtGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCCACCTTGGCGTGGCCACCCCCCCCCTGTTACCGCCCGTCCGCCGCAGTGGCGTCCTCGGcgt                                                              >  2:322070/1‑134 (MQ=255)
                                                                     aCCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGGCCCCGCCCCCCGTGGCGTGGGCCCCGCCGCCCGGTTCCCGCGCGTCCGCCGCAGGGGCGTCCGCGGCgtttc                                                           >  1:299218/1‑134 (MQ=255)
                                                                     aCCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCAACGTTGGCGTGGCCCCCCCCGACCGGTTCCCGCGCGTCAGCCTCACTGGCGTCCTCGGCGTTAc                                                           >  1:46841/1‑136 (MQ=255)
                                                                             gAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCCCCGTTGGCGTGGGCCCCCCCCCCCGGTGCCCGCGCGGCCGCCGCACGGGGGTGCTCGGCTGTTACGccgccc                                                   >  2:168633/1‑136 (MQ=255)
                                                                              aGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGccgccc                                                   <  6:186162/135‑1 (MQ=255)
                                                                                ctgctgCGCCGCCTCCCGGACATCCGCAGCGCCTCGCGCCAGCTGGCGGGGGACACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCTCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCcgcg                                                <  4:234779/136‑1 (MQ=255)
                                                                                       gccgccgccCGGACATCCGCAGCGCCGAGCGCCAGCGGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGAGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCa                                         <  3:16725/136‑1 (MQ=255)
                                                                                                       ccGCAGCGCCGAGCGCCCGCTGGCGGCGGCCACCGCCAACGTGGGGGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTcggcg                         <  4:24662/136‑1 (MQ=255)
                                                                                                       ccGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTcggcg                         <  2:153019/136‑1 (MQ=255)
                                                                                                       ccGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTcggcg                         <  3:188393/136‑1 (MQ=255)
                                                                                                       ccGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTcggcg                         <  2:295508/136‑1 (MQ=255)
                                                                                                         gCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCggcggc                       <  8:291358/136‑1 (MQ=255)
                                                                                                         gCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCggcggc                       <  8:286251/136‑1 (MQ=255)
                                                                                                         gCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCggcggc                       <  7:130100/136‑1 (MQ=255)
                                                                                                         gCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCggcggc                       <  7:71735/136‑1 (MQ=255)
                                                                                                          cAGCGCCGAGCGCCAGCTGGCGGCGGCCCCCGCCAACGTTGGCGTGGCCACCGCCGCCCTGTTCCCGCGCGTCAGCCTCAGTGGCGTCCTCGGCTGTACCGCCGCCCGCGGCGCGCAGATCGGCGCTTCGGCGGcc                      >  3:167260/1‑136 (MQ=255)
                                                                                                               ccGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGc                 <  1:283115/136‑1 (MQ=255)
                                                                                                                cGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCACTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGc                 >  8:116653/1‑135 (MQ=255)
                                                                                                                    cgcCCGCGGGCGGCGGCCACCCCCACCGTGGGGGGGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTggg            <  5:45719/136‑1 (MQ=255)
                                                                                                                    cgcCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTggg            <  2:17594/136‑1 (MQ=255)
                                                                                                                    cgcCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTggg            <  4:93456/136‑1 (MQ=255)
                                                                                                                    cgcCAGCTGGCGGCGGCCACCGCCAACGTGGGCGGGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTggg            <  7:9809/136‑1 (MQ=255)
                                                                                                                     gcCAGCTGGCGGCGGCCCCCGCCAACGTTGGCGTGGCCCCCGCCGCCCGGTTCCCGCGCGTCAGCCTCAGTGGCGTCCTCGGCTTTACCGCCGCCCGCGGCGCGCAGATCGGCGCTTAGGCGGCCCACGCCGGGGc           >  8:165562/1‑136 (MQ=255)
                                                                                                                     gcCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGCCCGGTTCCCGCGCGTCAGCCGCAGTGGCTTGCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCGTCGGCGGCCAACGCCGGGGc           >  6:164801/1‑136 (MQ=255)
                                                                                                                     gcCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGc           >  7:192740/1‑136 (MQ=255)
                                                                                                                          cgggcggcggcCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTgg      <  5:9222/134‑1 (MQ=255)
                                                                                                                          cTGGCGGCGGCCACCTCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTgg      <  1:53365/136‑1 (MQ=255)
                                                                                                                          cTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTgg      <  6:68351/136‑1 (MQ=255)
                                                                                                                          cTGGCGGCGGCCACCGCCAACGTGGGCGTGGCCACCGCCGACCTGTGCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTgg      <  3:373459/136‑1 (MQ=255)
                                                                                                                           tGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTggg     <  1:202874/136‑1 (MQ=255)
                                                                                                                           gggcggcggcCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTggg     <  4:120310/135‑1 (MQ=255)
                                                                                                                           gggcggcggcCACCGCCAACGTGGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTggg     <  5:44531/135‑1 (MQ=255)
                                                                                                                              cggcggcCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGccc  <  5:356075/136‑1 (MQ=255)
                                                                                                                                      |                                                                                                                               
CTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCC  >  ACB143/3875727‑3875988

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