BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,435,075	T→A	9.6%	G144G (GGT→GGA)	dapH →	2,3,4,5‑tetrahydropyridine‑2,6‑dicarboxylate N‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,435,075	0	T	A	9.6%	581.8 / 7.9	196	G144G (GGT→GGA)	dapH	2,3,4,5‑tetrahydropyridine‑2,6‑dicarboxylate N‑acetyltransferase
Reads supporting (aligned to +/- strand): ref base T (68/105); new base A (21/2); total (89/107)
Fisher's exact test for biased strand distribution p-value = 1.74e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 197 total aligned reads displayed.

GATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATT  >  USA300TCH1516_ALE/1434985‑1435164

 atgGGCGCAACAATTAATATTGGCGCGGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTa                                                                                          >  2:2903080/1‑91 (MQ=255)

 atgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTa                                                                                          <  1:1435643/91‑1 (MQ=255)

 atgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTa                                                                                          >  1:1437183/1‑91 (MQ=255)

  tgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaa                                                                                         >  2:3186435/1‑91 (MQ=255)

    ggCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACCACTGGAaaaa                                                                                       >  2:3634212/1‑91 (MQ=255)

    ggCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaaa                                                                                       >  1:747735/1‑91 (MQ=255)

     gcgcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaaaa                                                                                      <  1:2776796/91‑1 (MQ=255)

     gcgcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaaaa                                                                                      <  2:1633528/91‑1 (MQ=255)

       gcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATg                                                                                    <  1:1312401/91‑1 (MQ=255)

       gcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATg                                                                                    <  1:785473/91‑1 (MQ=255)

        caacaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGt                                                                                   <  2:3642777/91‑1 (MQ=255)

            aaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACAt                                                                               >  1:3459208/1‑91 (MQ=255)

            aaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGCTATAATTGGAAAAAAAGTACAt                                                                               >  2:1672117/1‑91 (MQ=255)

             tttAATATTGGCGCGGTCGTTGGCGAAGGTACAATGATGGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATg                                                                              <  2:1711142/90‑1 (MQ=255)

             aTTAATATTGGCGCAGTCGTTGGCGAGGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATg                                                                              <  2:1352272/91‑1 (MQ=255)

              ttAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGt                                                                             >  1:3129207/1‑91 (MQ=255)

              ttAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGt                                                                             <  1:2345512/91‑1 (MQ=255)

                 atatTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTCGTATCATTAGAAAAAAAGTACATGTAgg                                                                          >  2:2722047/1‑91 (MQ=255)

                  tggtGGCGCAGGCGTGGGCGAAGGGACAATGATTGATATGGATGCTACGCTCGGTGGTCGTGCTACAACTGGAAAAAATGTACATGTAggg                                                                         <  2:1209386/88‑1 (MQ=255)

                  tatTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAggg                                                                         <  1:1165463/91‑1 (MQ=255)

                  tatTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAggg                                                                         >  1:115533/1‑91 (MQ=255)

                   atTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAgggg                                                                        <  1:2077856/91‑1 (MQ=255)

                    ttGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGc                                                                       <  1:3033716/91‑1 (MQ=255)

                    ttGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGATAAAAAAAAAAAAAATGTACATGTAGGGGc                                                                       <  2:335715/91‑1 (MQ=255)

                     tGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAATTGGTAAAAATGTACATGTAGGGGCt                                                                      >  2:582490/1‑91 (MQ=255)

                      ggCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTg                                                                     <  1:3667767/91‑1 (MQ=255)

                      ggCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTg                                                                     <  2:1764616/91‑1 (MQ=255)

                      ggCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTg                                                                     >  1:3811523/1‑91 (MQ=255)

                       gcgcAGTCGTTGGCGGCGGTAGCATGATTGATATGAATGCTACCATCGGCGGACGTGCTCCAACTGGTAAAGATGTACATGTAGGGGCTgg                                                                    <  2:2191003/91‑1 (MQ=255)

                       gcgcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTgg                                                                    <  1:3346667/91‑1 (MQ=255)

                       gcgcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTgg                                                                    <  1:2960711/91‑1 (MQ=255)

                         gcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcg                                                                  <  2:2487428/91‑1 (MQ=255)

                          cAGTCGTTTGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTGGAAAAAATATACATGTAGGGGCTGgcgc                                                                 >  2:1015481/1‑91 (MQ=255)

                          cAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcgc                                                                 >  1:2858725/1‑91 (MQ=255)

                          cAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcgc                                                                 <  1:186390/91‑1 (MQ=255)

                          cAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGCTACAACTGGAAAAAATGTACATGTAGGGGCTGgcgc                                                                 >  2:3428913/1‑91 (MQ=255)

                            gTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAg                                                               >  1:2548168/1‑91 (MQ=255)

                            gTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAg                                                               >  1:3859308/1‑91 (MQ=255)

                              cGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTa                                                             <  1:2915781/91‑1 (MQ=255)

                                tttgCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTAtt                                                           <  2:1045775/88‑1 (MQ=255)

                                  ggCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGGTACCACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAg                                                         >  2:1318955/1‑91 (MQ=255)

                                   gCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGc                                                        >  1:1801686/1‑91 (MQ=255)

                                    cGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                       <  1:1406788/91‑1 (MQ=255)

                                    cGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTAGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                       >  1:2437822/1‑91 (MQ=255)

                                      aaGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTGAAAAAAAAGTACCCGTAGGGGCTGGCGCAGCATTAGCAgg                                                     >  2:1265830/1‑91 (MQ=255)

                                       aGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgt                                                    <  2:103827/91‑1 (MQ=255)

                                        ggTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtg                                                   <  1:2582394/91‑1 (MQ=255)

                                         gTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtgt                                                  <  1:3565059/91‑1 (MQ=255)

                                           gcAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGa                                                <  2:3846189/90‑1 (MQ=255)

                                           aCAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGa                                                >  1:2037225/1‑91 (MQ=255)

                                           aCAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGa                                                <  2:2098668/91‑1 (MQ=255)

                                           aCAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGa                                                <  1:2509234/91‑1 (MQ=255)

                                           aCAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTACAATTAGTAAAAATATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGa                                                >  2:2404209/1‑91 (MQ=255)

                                            cAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                               <  1:1681020/91‑1 (MQ=255)

                                            cAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                               <  2:1391098/91‑1 (MQ=255)

                                            cAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                               >  1:1753616/1‑91 (MQ=255)

                                            cAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                               <  2:724199/91‑1 (MQ=255)

                                            cAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTACCACTGGTAAAAAAGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                               >  2:194770/1‑91 (MQ=255)

                                               tgattgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGa                                            <  2:2396232/91‑1 (MQ=255)

                                                gattgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGaa                                           <  1:1595205/91‑1 (MQ=255)

                                                 attgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAc                                          <  1:132613/91‑1 (MQ=255)

                                                    gatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                                       <  1:400777/91‑1 (MQ=255)

                                                    gatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                                       >  1:3263346/1‑91 (MQ=255)

                                                    gatATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAAGAGTGACTGAAcccc                                       >  2:200217/1‑91 (MQ=255)

                                                     atatGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAccccc                                      <  1:2935795/91‑1 (MQ=255)

                                                      tatGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCt                                     >  1:880487/1‑91 (MQ=255)

                                                      tatGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCt                                     <  1:2687901/91‑1 (MQ=255)

                                                      tatGAATGCTACTCTCGGTGGTCGCGCTACCACTAGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCt                                     >  2:1657021/1‑91 (MQ=255)

                                                      tatGAATGCTACTCTCGGTGGTCGCCGTACCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCt                                     >  2:3353999/1‑91 (MQ=255)

                                                       atGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTa                                    <  1:3445693/91‑1 (MQ=255)

                                                        tGAATGCTACTCTCGGTGGTCGCCGTACCATTAGTAAAAAAGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAg                                   >  2:509107/1‑91 (MQ=255)

                                                          aaTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTg                                 <  1:2910274/91‑1 (MQ=255)

                                                          aaTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTg                                 <  1:2225890/91‑1 (MQ=255)

                                                            tGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCt                               <  2:809656/91‑1 (MQ=255)

                                                            tGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCt                               <  2:2060324/91‑1 (MQ=255)

                                                              cTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTc                             >  1:2043885/1‑91 (MQ=255)

                                                               tACTCTCGGTGGTCGTGGTACCACTAGAAAAAAAGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTcc                            >  2:2807988/1‑90 (MQ=255)

                                                                  tctcGGTGGTCGTGCTACCACTAGAAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCg                         >  2:1438476/1‑91 (MQ=255)

                                                                   ctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCtt                              >  1:3151453/1‑85 (MQ=255)

                                                                   ctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCtt                              <  2:3151453/85‑1 (MQ=255)

                                                                    tcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGt                       <  1:1829481/91‑1 (MQ=255)

                                                                    tcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGt                       <  1:2277182/91‑1 (MQ=255)

                                                                    tcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGt                       <  2:2676534/91‑1 (MQ=255)

                                                                     cGGTGGTCGCCGTACCATTGGAAAAAATATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGtt                      >  2:2448447/1‑91 (MQ=255)

                                                                      ggtggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTta                     >  1:2510795/1‑91 (MQ=255)

                                                                      ggtggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGGGTGATTGAACCCCCTAGTGCTTCACCGGTta                     >  1:3500877/1‑91 (MQ=255)

                                                                      ggtgggCGTGCTACGACTGGTTAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCGCCGGTta                     <  2:1588099/91‑1 (MQ=255)

                                                                       gtggtCGTGCTACAGCTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTtat                    <  2:2577438/91‑1 (MQ=255)

                                                                          gtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAAt                 >  1:1519621/1‑91 (MQ=255)

                                                                           tCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATc                >  1:3422244/1‑91 (MQ=255)

                                                                            cGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCg               <  1:521768/91‑1 (MQ=255)

                                                                             gTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGa              <  1:98192/91‑1 (MQ=255)

                                                                             gTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGa              >  1:69641/1‑91 (MQ=255)

                                                                             gTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGa              >  1:1248883/1‑91 (MQ=255)

                                                                                cTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGga           >  1:3790581/1‑91 (MQ=255)

                                                                                  aCAACTGGTAACAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGgatg         <  1:1983601/91‑1 (MQ=255)

                                                                                     aCTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGatgatg      <  1:1065265/91‑1 (MQ=255)

                                                                                      cTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGt     >  1:2056085/1‑91 (MQ=255)

                                                                                        ggTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTAt   <  2:1834357/91‑1 (MQ=255)

                                                                                        ggTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTAt   <  2:945105/91‑1 (MQ=255)

GATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATT  >  USA300TCH1516_ALE/1434985‑1435164

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