Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,435,0650CG24.5% 244.3 / 25.1 107A141G (GCT→GGT) dapH2,3,4,5‑tetrahydropyridine‑2,6‑dicarboxylate N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (32/48);  new base G (26/0);  total (58/49)
Fisher's exact test for biased strand distribution p-value = 6.04e-09
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 107 total aligned reads displayed.

CGAATTGAACCAGGTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTGGTAA  >  USA300TCH1516_ALE/1434923‑1435209
                                                                                                                                              |                                                                                                                                                
cGAATTGAACCAGGTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGGTACaa                                                                                                                                             >  2:1162886/1‑148 (MQ=255)
cGAATTGAACCAGGTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACaa                                                                                                                                             >  1:268185/1‑148 (MQ=255)
  aaTTGAACGAGGTGCCTTTATTAGAGAACAAGACATTATCGAAGATGGTTCTGTGGTTATGATGGGCGCAAAGATGAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACt                                                                                                                                           <  2:250464/148‑1 (MQ=255)
      gAACCAGGTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTa                                                                                                                                       <  1:589220/148‑1 (MQ=255)
         ccAGGAGCGTTTATTAGAGAAAAAGACATTATTGAAGATGGTGCTGTCGTTATTATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaaa                                                                                                                                    <  2:480357/148‑1 (MQ=255)
             gTGCGTTTATTAGAGAACCAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGt                                                                                                                                <  1:1071505/148‑1 (MQ=255)
              tGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTa                                                                                                                               <  1:582662/148‑1 (MQ=255)
                cGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACa                                                                                                                             <  1:1142718/148‑1 (MQ=255)
                 gTTTATTAGATAACAAGCCATTATTGAAGATGGTGCTGGCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACAt                                                                                                                            <  1:1196064/148‑1 (MQ=255)
                      ttaGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGCTATAACTAGAAAAAAGATACATgcaca                                                                                                                       >  2:646412/1‑144 (MQ=255)
                       taGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGAGATGAATGCTACTCTCGGTGGTCGTGCTGCAACTGGTAAAAATGTACATGTAggg                                                                                                                      >  1:1335402/1‑148 (MQ=255)
                        agagAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAgggg                                                                                                                     >  1:1216444/1‑148 (MQ=255)
                        agagAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTCGTATCATTAGAAAAAAAGTACATGTAgggg                                                                                                                     >  2:694331/1‑148 (MQ=255)
                        agagAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTCCAATTAAAAAAAAAATACATTTAgggg                                                                                                                     >  2:1200304/1‑148 (MQ=255)
                         gagaACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAAGGTACATGTAGGGGc                                                                                                                    >  1:1235550/1‑148 (MQ=255)
                          agaACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCt                                                                                                                   >  1:791070/1‑148 (MQ=255)
                            aaaaaGCCATTATTGAAGATGGTGCTGTCGTTTTGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACACTGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTgg                                                                                                                 <  2:786073/145‑1 (MQ=255)
                            aaCGAGCCATTATTTAAGATGGTGCTTGCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTgg                                                                                                                 <  2:803735/148‑1 (MQ=255)
                               aaGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACAACTAAAAAAAAAATACATGAAGGGGCTGgcgc                                                                                                              >  2:1149/1‑148 (MQ=255)
                                 gcaatttttGAAGATGGTGGTGTCGTTAGGATGGGCGCAACAATTAATATTGGGGCAGTCGTTGGCGAAGGTAAAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAg                                                                                                            <  1:37915/145‑1 (MQ=255)
                                 gCCATTATTGAAGATGGTGCTGTCGTTATGAGGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAg                                                                                                            <  1:582414/148‑1 (MQ=255)
                                 gCCATTAGTGAATAAGGTGGTGTCGTTATGAAGGGCGCAACAATTAATATTGGCACAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAg                                                                                                            <  2:135826/148‑1 (MQ=255)
                                   caattatGAAGATGGGGCTATCGGTATGATGGGCGTAACAATTAATGTTGGCGCAGTCGTTGTCGAAGGTACAATGATGGATATGAATGCTACTCGCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTa                                                                                                          <  2:707239/142‑1 (MQ=255)
                                   cATAATAGAAGATGGTGCTGTCGTAATGCAGGGCGCAACAATGGATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGACTGCTTCTGTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTa                                                                                                          <  2:1285435/148‑1 (MQ=255)
                                    attattGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGTAGTACCGATGAAAGGTAACGTGCATTACACGAAGGCTACTCTCGTTGATCGTGCTACGACTGATAGAAATGTACATGTGGGGGCTGGCGCAGTAt                                                                                                         >  1:1392300/1‑148 (MQ=255)
                                       attTAAGATGGTTCTAGCGGTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGGGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAg                                                                                                      <  2:1184297/148‑1 (MQ=255)
                                       attGAAGATGGGGCATTCGTTATGATGGGTGCAACAATTAATGTTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAg                                                                                                      <  2:1338680/148‑1 (MQ=255)
                                        ttGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTCGGGGCTGGCGCAGTATTAGc                                                                                                     >  1:179423/1‑148 (MQ=255)
                                        ttGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGc                                                                                                     >  1:1164182/1‑148 (MQ=255)
                                        ttGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGc                                                                                                     >  1:806476/1‑148 (MQ=255)
                                         tGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                                                                    >  1:1167532/1‑148 (MQ=255)
                                         tGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                                                                    >  1:420433/1‑148 (MQ=255)
                                         tGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                                                                    >  1:1329508/1‑148 (MQ=255)
                                         tGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAAAACATATAGGGGCTGGCGCAGTATTAGcg                                                                                                    >  2:25179/1‑147 (MQ=255)
                                          gAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCATTGAAAAAAATGTACATGTAGGGGCTGGCACAGTATTAGCAg                                                                                                   >  2:169231/1‑148 (MQ=255)
                                           aaGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAAGTAAAAAAAAAATACATGTAGGGGCTGGCGAAGTAAGAGCacg                                                                                                  >  2:1227963/1‑146 (MQ=255)
                                            aGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACACGTAGGGGCTGGCGCAGTATTAGCAGgt                                                                                                 >  1:661539/1‑148 (MQ=255)
                                            aGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGgt                                                                                                 >  2:655407/1‑148 (MQ=255)
                                            aGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAggg                                                                                                 >  2:894240/1‑147 (MQ=255)
                                            aGATGGGGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATAATTGATATGAATGCGACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgt                                                                                                 <  2:1412530/148‑1 (MQ=255)
                                               tGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGtgtg                                                                                              <  1:1174347/148‑1 (MQ=255)
                                                 gtccTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATTAATGCGACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCATTATTAGCAGGTGTGAt                                                                                            <  1:695697/145‑1 (MQ=255)
                                                 gTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                                                                            <  1:89040/148‑1 (MQ=255)
                                                   gCTGTCGTTATGATGGGCGCAACAATTAATATGGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGGATTAGCAGGTGTGAGTg                                                                                          >  1:146675/1‑148 (MQ=255)
                                                      gTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAc                                                                                       >  1:549060/1‑148 (MQ=255)
                                                      gTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGAGGCTGGCGCAGTATTAGCAGGTGTGATTGAAc                                                                                       >  1:243900/1‑148 (MQ=255)
                                                        cGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAccc                                                                                     >  1:600991/1‑148 (MQ=255)
                                                         gTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGTTATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                                                                                    <  2:755301/148‑1 (MQ=255)
                                                         gTTATGATGGGCGCAACAATTAACATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                                                                                    <  1:55719/148‑1 (MQ=255)
                                                         gTTATAAGGGGCGTAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                                                                                    <  2:840938/148‑1 (MQ=255)
                                                              gatgGGCGCAACAATTAATATTGGCGCAGTCGTTGGAGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGt                                                                               >  1:894965/1‑148 (MQ=255)
                                                                    cgcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCATTAAAAAAAAAAAACATGTAGGGGCCGGCGCAGTACTAGCAAGATTGATTGAACCCCCTACTGcagca                                                                         >  2:1087074/1‑144 (MQ=255)
                                                                     gaaacaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCAc                                                                        <  1:591901/146‑1 (MQ=255)
                                                                     acaacaaATAAAACTAGCGCAGTCGTTAGCAAAGAAACAACGATTAATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCAc                                                                        <  1:375753/147‑1 (MQ=255)
                                                                         caaTTAATATTGGTGCGGTGGTTGGCGAAGGTACAAGGATTGATATGAATGCTGCCCTCGGGGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGTTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGt                                                                    <  1:46820/148‑1 (MQ=255)
                                                                          caTTAATATGTGCCAAGTCGTTGGCGAAGGTAAAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTCGCAGGCGCGATTGCACCCCCTAGTGCTTCACCGGtt                                                                   <  2:1398803/147‑1 (MQ=255)
                                                                          aaTTAATATGGGCGCAGTCGTTGGCGAGGGTGCAATGATTGAGATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGtt                                                                   <  2:220957/148‑1 (MQ=255)
                                                                               atatTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAAt                                                              <  1:156017/148‑1 (MQ=255)
                                                                               atatTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCCTATCACTAATAAAAATATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAAt                                                              >  2:418411/1‑148 (MQ=255)
                                                                                   tGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTCGTACCATTAAAAAAAAAGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAg                                                          >  2:1245472/1‑148 (MQ=255)
                                                                                   tGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCACCTAGTGCTTCACCGGTAATAAGCGAg                                                          >  2:1069349/1‑148 (MQ=255)
                                                                                     gcgcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGAATTATCAGGTGGGATTGAAGCCCCTAGTGCTTCAACGGTTATAAACGAGga                                                        >  2:211889/1‑148 (MQ=255)
                                                                                         aGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTCTAATCGAGgatgat                                                    >  1:113084/1‑148 (MQ=255)
                                                                                         aGTCGTTGGCGAAGGTACAATGATTGAAATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACTTGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGgatgat                                                    <  1:979172/148‑1 (MQ=255)
                                                                                            cGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGAGGATGTa                                                 >  1:766570/1‑148 (MQ=255)
                                                                                               tGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAAGCGAGGATGATGTATTa                                              <  2:1005977/148‑1 (MQ=255)
                                                                                                 gCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAAGAGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTaaa                                            >  2:254310/1‑147 (MQ=255)
                                                                                                  cGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCGGTATTAGCAGGTGTGATTGAACCCCCTAGAGCTGCACCGGTTATAAACGAGGATGATGTATTAATc                                           >  2:187177/1‑148 (MQ=255)
                                                                                                   gAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGTTGTGATTGAACCCCCTAGTGCTTCCCCGGTTATACTCGAGGATGATGTATTAATCg                                          >  2:1012947/1‑148 (MQ=255)
                                                                                                   gAAGGTACAATGATTGATATGAATGCTACTCTCAGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGGTGATGTATTAATCg                                          >  1:703971/1‑148 (MQ=255)
                                                                                                    aaGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCgg                                         >  1:1236822/1‑148 (MQ=255)
                                                                                                    aaGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGAAGTATTAATCgg                                         >  1:593190/1‑148 (MQ=255)
                                                                                                     aGGTAAAATGATTGAGAGGAATGCGACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGt                                        <  2:456353/148‑1 (MQ=255)
                                                                                                        tACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTGAAAAAAAAATACATATAGGGGCTGGAGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCTGTTATAATAGAGAACGATGGATTAATCgggtca                                     >  2:1111803/1‑144 (MQ=255)
                                                                                                         aaaatgaatgATATAAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCaa                                    <  2:1039641/146‑1 (MQ=255)
                                                                                                          cAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAGTACATGTAGGGGCGGGCGCAGTATTAGCAGGTGTGATTGAACCCGCTAGTGCTAAAGCGATAATAATCGAGGATGATGTATTAATCGGTGCaac                                   >  2:1053515/1‑147 (MQ=255)
                                                                                                              gattaacATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCa                               <  1:599358/148‑1 (MQ=255)
                                                                                                              gaatgaaATAAATGCAACACGCGGTGGTCGGGCTACAACTGGTAAAAATGTACATGAAGGGGCTCGCTCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCa                               <  2:670693/145‑1 (MQ=255)
                                                                                                                ttgatATGAATGCTCCTCTCGGTGGTCGTGCTACAACTCGTAAAAATGTACATGTAGAGACTGTCGCATTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGt                             <  1:1386660/148‑1 (MQ=255)
                                                                                                                     atGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGGAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTAttt                        <  1:1149/148‑1 (MQ=255)
                                                                                                                     atGAAGGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTAttt                        <  1:981270/148‑1 (MQ=255)
                                                                                                                      tGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTAtttt                       >  1:1413360/1‑148 (MQ=255)
                                                                                                                       gAATGCTACTCTCGGTGGTCGTCGTATCATTGAAAAAAAAGTACATGTAGGGACTGGCGCAGCATTCGAAGGTGTGCTTGAACCCCCTAATGGTTGCACGGTTGTAATGGAGGATGATGTATTAATCGTTGCAAATGAAGATATTTTa                      >  2:616418/1‑148 (MQ=255)
                                                                                                                       gAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCACTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCGAATGCAGTTATTTTa                      >  2:937542/1‑148 (MQ=255)
                                                                                                                       gAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGAAAGTGAGATTGAACCCCCTATTGCTTCTCCGGTTATAATCGAGGGTGATGTCTTATTAGGTGCAAATGGAGTTAGTTTa                      >  2:937691/1‑148 (MQ=255)
                                                                                                                        aaTGTTACTCTCGGTGGTCGTGCTACAACTTGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGGACCGCCTAGTGCTTCACCGGTTATAATCGAGGATGGTGTATTAATCGGTGCAAATGCAGTTGtttggg                     <  2:870104/148‑4 (MQ=255)
                                                                                                                         aTGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTGCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGa                    >  2:251698/1‑148 (MQ=255)
                                                                                                                           gatACTCGCGGGGGTCGTGCTACCACTGGTAAACATGTACTTGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAg                  <  1:934387/146‑1 (MQ=255)
                                                                                                                           gCTACTCTCGGTGGTCGCCGTACCATTAATAAAAACATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAg                  >  2:1360233/1‑148 (MQ=255)
                                                                                                                           gCTACTCACGGTGGTAGTGATAAAACTGGTAAAAATGTACATGTTGGGGCTGGCGCAGTATTAGCAGGTGTGATTGATGCCCCTAGTGCTTCACCGGTTATAATAGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAg                  <  2:1056264/148‑1 (MQ=255)
                                                                                                                            cTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAgg                 >  1:1149075/1‑148 (MQ=255)
                                                                                                                            cTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGGTGATGTATTACTCGGTGCCAATGAAGTTATTTTAGAAgg                 >  2:373020/1‑148 (MQ=255)
                                                                                                                                 ctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTAc            >  1:1376223/1‑148 (MQ=255)
                                                                                                                                      tggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACgtgtg       >  1:616188/1‑147 (MQ=255)
                                                                                                                                        gtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCTCCTAGTGCTTCACCGGTTATAATAGAGGATGATGTATTAATCTGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTgg     <  2:682586/148‑1 (MQ=255)
                                                                                                                                        gtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTgg     <  1:1238844/148‑1 (MQ=255)
                                                                                                                                        gtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTgg     <  2:1164182/148‑1 (MQ=255)
                                                                                                                                        gtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTgg     <  1:1028025/148‑1 (MQ=255)
                                                                                                                                        gacGTGCTACAACTGGTAAAAATGAACAGGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGttag     <  2:370597/146‑3 (MQ=255)
                                                                                                                                           gTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATTTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACgtgtcgccca  >  2:961242/1‑142 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CGAATTGAACCAGGTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTGGTAA  >  USA300TCH1516_ALE/1434923‑1435209

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: ‑