Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,435,0650CG20.0% 200.2 / 10.9 76A141G (GCT→GGT) dapH2,3,4,5‑tetrahydropyridine‑2,6‑dicarboxylate N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (21/39);  new base G (15/0);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 2.44e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTGAACCAGGTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTGGTAAAG  >  USA300TCH1516_ALE/1434926‑1435211
                                                                                                                                           |                                                                                                                                                  
aTTGAACCAGGTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATGGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTgg                                                                                                                                           >  1:897611/1‑149 (MQ=255)
 ttcaaacAGGTGCGCTTATTAGAGAACAAGCCATTAGTGAAGTTGGGGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCGGTGGTAGTCGAGGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGt                                                                                                                                          <  2:401770/143‑1 (MQ=255)
    aaCCCGGTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAGGGTACAATGTTTGATATGAATGATACTCTCGGTGGTCGTGCTACAAATGGTaaa                                                                                                                                       >  1:520997/1‑149 (MQ=255)
    aaCCAGGTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACCACTGGTaaa                                                                                                                                       >  1:191083/1‑149 (MQ=255)
    aaCCAGGTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGGCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaa                                                                                                                                       <  1:550280/149‑1 (MQ=255)
     aCCAGGTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaaa                                                                                                                                      <  1:338186/149‑1 (MQ=255)
     aCCAGGTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCCACTCTCGGTGGTCGTGCTACAACTGGTaaaa                                                                                                                                      <  1:47898/149‑1 (MQ=255)
        aGGTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTGGTGCTACAACTGGTAAAAATg                                                                                                                                   >  1:925293/1‑149 (MQ=255)
         ggggCGTTGATTAGAGAACAAGCCATTATTGACGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAGGGCACAATGGATGATAGGAATGCTACTCTCGGTGGTCGTGGTATCACTGGCaaaaaaga                                                                                                                                  >  2:458683/4‑146 (MQ=255)
         ggggCGGTTATCAGGGAACAAGCCATGATTGGGGAGGGTGCTGTGGGTATGATGGGCGCAACAATTAATATTGGCGGAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGt                                                                                                                                  <  2:28481/146‑1 (MQ=255)
         ggTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGt                                                                                                                                  >  1:486375/1‑149 (MQ=255)
           tGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGGTACAACTGGTAACAATGTaa                                                                                                                                >  1:853126/1‑148 (MQ=255)
           tGCGTTTATTAGAGAACAAGCCATTATCGAGGCTGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAACATGTAc                                                                                                                                >  1:1102733/1‑149 (MQ=255)
             cGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCGGGCTTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACAt                                                                                                                              <  1:483069/149‑1 (MQ=255)
                  attaGAGAACAAGCCATTATTGAAGATGGTGCTGGCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACAACTGGTAAAAATGTACATGTAgg                                                                                                                         >  2:693791/1‑149 (MQ=255)
                      gagaACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCATTGGTAAAAAAGTACATGTAGGGGCt                                                                                                                     >  2:184158/1‑149 (MQ=255)
                       agaACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTATCATTAATAAAAAAAAACATGTAGGGGCTg                                                                                                                    >  2:365576/1‑149 (MQ=255)
                         aaCAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgc                                                                                                                  >  1:1010377/1‑149 (MQ=255)
                         aaCAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTACCATTAAAAAAAAAGTATATATAAGGGCTGgc                                                                                                                  >  2:733563/1‑149 (MQ=255)
                          aCAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcg                                                                                                                 <  1:432566/149‑1 (MQ=255)
                          aCAAGCCATTATTGAAGATGCTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTACAAATGTACATGTAGGGGCTGgcg                                                                                                                 <  1:825130/149‑1 (MQ=255)
                           aaacatattattGAAAATGGTGATGCAGTTATGATGGGCACAAATATTAATATTGGCGCAGTCGTTGGCGAAGGTAGAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcgc                                                                                                                <  2:1045951/143‑1 (MQ=255)
                            caGCCACTATTGAACAACGTGCTCTCGTTATGATGGGCGACATATTTAATCTTTTCGCTGTCGTTGTCGAGGGTATAATGATCGATATGAATGCTACCCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCa                                                                                                               <  2:349954/148‑1 (MQ=255)
                            aaGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCa                                                                                                               <  1:372861/149‑1 (MQ=255)
                             aGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAg                                                                                                              <  1:33577/149‑1 (MQ=255)
                              ggcATTATTGAAGTTGGTGCAGTCGTGTGGATGGGCGAAACAATTTATATTGGCGCAGGCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGt                                                                                                             <  2:17389/147‑1 (MQ=255)
                              gCCATTAATAAAATTGGCGCTGTTGTAATGATGGGCGGAACAAGTAATATTTGTGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTAAAACTGGTAAAACTGTACACGTAGGGGCTGGCGCatt                                                                                                             <  2:380207/149‑3 (MQ=255)
                               ccATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTa                                                                                                            >  1:972679/1‑149 (MQ=255)
                               ccATTAGTGAAGATGGGGCGGTCGGTATGATGGAGGCAAGAGTTTTCAGTGGCGCAGTCGTTGGGGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTa                                                                                                            <  2:235268/149‑1 (MQ=255)
                                        aacaTGGTGCTGTCGTTACGATGGCCGCAACAAGTAATATTGGCGCAGTCCTTGGCGAAGGTACAATGATTGATATGACTGCTACTCTCGGTGGTCGTGCTACAACTGGTAACAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgt                                                                                                   <  1:112362/146‑1 (MQ=255)
                                        aaGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCATTGAAAAAAAAGTACATGTAGGGGCTGGCGCAGTATTAGCAggg                                                                                                   >  2:884882/1‑148 (MQ=255)
                                        aaGATGGTGCTGCCGTTATGAGGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgt                                                                                                   <  1:13634/149‑1 (MQ=255)
                                         aGATGGGGCTGACGTAAAGAAGAGAGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtg                                                                                                  <  2:432792/149‑1 (MQ=255)
                                              gTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTACCATTGAAAAAAAAGTACATGTAGGGGCTGGCGGAGTATTAGCAGGTGTGAtg                                                                                             >  2:409750/1‑148 (MQ=255)
                                               tGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCCTACCACTAGAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTg                                                                                            >  2:743054/1‑149 (MQ=255)
                                                gCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGa                                                                                           >  1:733012/1‑149 (MQ=255)
                                                 ttgacGTTATGATGGGCGCAACAATTAATATTGGAGCAGTCGTTGGCGAAGGTACAATGATTGACATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGaa                                                                                          <  2:808003/145‑1 (MQ=255)
                                                 tctgcGTAATGAGGGGCGCAACAATTAATATTTGCGCAGTCGTTGGGGAAGGTTCAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTATAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTGGCAGGTGTGATTGaa                                                                                          <  2:518790/145‑1 (MQ=255)
                                                  tGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAc                                                                                         >  1:512053/1‑149 (MQ=255)
                                                     ctctattatGGGTGCAATAATTAAGATTGGCCCAGTCGTTGGCGAATGTACAATGATTGATATGAATGCTACTCTGGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                                                                                      <  1:508444/146‑1 (MQ=255)
                                                     cGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                                                                                      <  1:908804/149‑1 (MQ=255)
                                                          tgatgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGGATTAACAGGTGAGATTGAACCCCCtcgg                                                                                 >  2:929789/1‑146 (MQ=255)
                                                                  gccacaaATAATATTGGCGCAGTCGTTTGCGCAGGTACAACGATTGATATGAATGCTACCCACGGTGGTGGTGAAAAAAATAGAAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCAcc                                                                         <  2:673980/146‑1 (MQ=255)
                                                                  gcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACCCTCGGTGGTCGTGCTACAACTTGTAAAAATGTACAAGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCAcc                                                                         <  1:235546/149‑1 (MQ=255)
                                                                  ccaataaTTAAGATTGGCGCAGTCGTTGGCGAAGGTACAATCATTGATATGATTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCAcc                                                                         <  2:701879/148‑1 (MQ=255)
                                                                     acaaTTAATGTTGGCGCAGTGGTTGGCGAAGGTACAATGATTGATATGAATGCTACGCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGt                                                                      <  2:188994/149‑1 (MQ=255)
                                                                           aaTATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGAAGGTGGGATTGAACCACCTAGTGCTTCACCGGTTATAAt                                                                >  2:125737/1‑149 (MQ=255)
                                                                            atatTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTCGTATCATTAAAAAAAAAGTACATGTAAGGGCTGGCACCGTATTAACCGGTGAGATGGAACCCCCT‑GTGCTTCACCGGTCATAAtgg                                                              >  2:487600/1‑147 (MQ=255)
                                                                              atTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCTTAGTGCTTCACCGGTTATAATCGa                                                             <  1:560471/149‑1 (MQ=255)
                                                                              atTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCGCAGTGCTTCACCGGTTATAATCGa                                                             <  1:560514/149‑1 (MQ=255)
                                                                                tGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAgg                                                           >  1:257297/1‑149 (MQ=255)
                                                                                tGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAgg                                                           >  1:173076/1‑149 (MQ=255)
                                                                                  gcgcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGgat                                                         >  1:68885/1‑149 (MQ=255)
                                                                                   cgcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAAAACATGTCGGGGCTGGACCAGTATTAGCAGGTGTGATTGAAACCCCTAGTGCTTCACCGGTTCTAATCGTGgaag                                                        >  2:613427/1‑147 (MQ=255)
                                                                                     cgcgcagtGGGGAAGATAAAATGAGTGATATGAATGCGACGGTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATGAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGgatgat                                                      <  2:473085/142‑1 (MQ=255)
                                                                                        tccttGTCGAATGTACAATGATTGATATGAAGGCTTCTCGCGGTGGTCGTGCTACAACTGGGAAAAATGTACAGGTAGGGGCTGGCGCAGTATTAGCAGGTGGGATTGATCCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTa                                                   <  2:873537/146‑1 (MQ=255)
                                                                                           ttGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTa                                                <  1:184158/149‑1 (MQ=255)
                                                                                              gccaaGGGACAATGATTGATAGGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGTAGGTGGGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATc                                             <  2:216993/146‑1 (MQ=255)
                                                                                               cGAAGGTACAAGGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCg                                            >  1:801674/1‑149 (MQ=255)
                                                                                                  aGGTACAATTATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTTGTAAAAAGGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTg                                         <  2:637229/149‑1 (MQ=255)
                                                                                                          tgtttgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGGGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCa                                 <  1:985109/146‑1 (MQ=255)
                                                                                                           tattgatATGAATGCTCCTCTCGGGGGGCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAg                                <  2:897611/148‑1 (MQ=255)
                                                                                                            attgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTGATCGGTGCAAATGCAGt                               <  1:818842/149‑1 (MQ=255)
                                                                                                               gatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTAt                            <  2:867890/149‑1 (MQ=255)
                                                                                                                ataaGAATGCTACTCTGGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGAAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTGAtt                           <  2:299491/149‑1 (MQ=255)
                                                                                                                 tatGAATGCTACTCTCGGTGGTCGCGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTAttt                          >  1:556315/1‑149 (MQ=255)
                                                                                                                 tatGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAAAACATGTGGGGGGTAGCGCAGAATTGGCAGATAAGATTGAACGCCCTAGTGCTGCACCGGTTATAATCGAGGATGACGTATCAATCAGTGCCAATTCAGTTAttt                          >  2:175817/1‑149 (MQ=255)
                                                                                                                    gAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCACAGTATTGGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAGTCGAGGATGATGTATTAACAGGTGCAACAGCCATTACTTTAg                       >  2:1096910/1‑149 (MQ=255)
                                                                                                                      aTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGaa                     >  1:15024/1‑149 (MQ=255)
                                                                                                                      aTGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAAAACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGACCCACGCGGTGATGCAGCGGTTATAATCGAGGATGATGTATTAATCGGTGCCACTGCAGTTTTTTTAGaa                     >  2:951958/1‑149 (MQ=255)
                                                                                                                           aCTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAAGGCAGTTATTTTAGAAGgtgt                >  1:313239/1‑149 (MQ=255)
                                                                                                                              ctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAAGCGGTGCAAATGCAGTTATTTTAGAAGGTGTACg             >  1:234686/1‑149 (MQ=255)
                                                                                                                                  gtggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGtt         <  1:693791/149‑1 (MQ=255)
                                                                                                                                   tggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTg        <  1:450005/149‑1 (MQ=255)
                                                                                                                                     gtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTGGt      >  2:703379/1‑149 (MQ=255)
                                                                                                                                         tGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATGATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTGGTAAAg  >  2:140423/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
ATTGAACCAGGTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTGGTAAAG  >  USA300TCH1516_ALE/1434926‑1435211

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑