BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	638,614	G→T	T14K (ACA→AAA)	pdxK ←	Pyridoxine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	638,614	0	G	T	98.1%	351.7 / ‑0.1	107	T14K (ACA→AAA)	pdxK	Pyridoxine kinase
Reads supporting (aligned to +/- strand): ref base G (2/0); new base T (52/53); total (54/53)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

Only 100 of 108 total aligned reads displayed.

CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATG  >  USA300TCH1516_ALE/638477‑638751

cATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTCTCCATGGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGATTTCAAATCTGCTTGCATACCTGCGCCAGCACTTttgt                                                                                                                                         >  2:521445/1‑141 (MQ=255)

 aTCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTc                                                                                                                                        <  2:420988/141‑1 (MQ=255)

     aTTGGTAATGGTGCAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGaa                                                                                                                                    <  1:532992/141‑1 (MQ=255)

       tGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAcc                                                                                                                                  <  1:84048/141‑1 (MQ=255)

       tGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAATTCACGGGACCGTGCCCGACGTATCTAATTCTTGGAACTTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAcc                                                                                                                                  <  1:302173/141‑1 (MQ=255)

        ggTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCg                                                                                                                                 <  3:186926/141‑1 (MQ=255)

         gTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCgg                                                                                                                                <  2:22330/141‑1 (MQ=255)

          tAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGc                                                                                                                               >  2:160440/1‑141 (MQ=255)

            aTGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTTTCAGAACCGACaa                                                                                                                             >  4:203457/1‑141 (MQ=255)

            aTGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCaa                                                                                                                             >  1:157100/1‑141 (MQ=255)

            aTGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCaa                                                                                                                             >  4:335731/1‑141 (MQ=255)

             tGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAGCGGCGAt                                                                                                                            >  4:230780/1‑141 (MQ=255)

                tgtAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGt                                                                                                                         <  2:183542/141‑1 (MQ=255)

                 gtAGCATCGTGTGACCTTGTATCTTTATCGATAGTAGCGATGGCAGTTAAAGGGACCATGGTATACGTATCGAATTCTTGGAACGTTTTCAAATCTGCTTGCAGACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGtt                                                                                                                        <  4:328178/141‑1 (MQ=255)

                 gtAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGtt                                                                                                                        >  1:373514/1‑141 (MQ=255)

                 gtAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGtt                                                                                                                        <  4:141002/141‑1 (MQ=255)

                   aaCATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑aa                                                                                                                     <  2:800194/141‑1 (MQ=255)

                    aCATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑aaa                                                                                                                    <  2:240224/141‑1 (MQ=255)

                      aTCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAAc                                                                                                                  <  4:180234/141‑1 (MQ=255)

                      aTCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAAc                                                                                                                  <  4:347999/141‑1 (MQ=255)

                      aTCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAAc                                                                                                                  <  4:139762/141‑1 (MQ=255)

                         gtgtGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAActtt                                                                                                               >  2:260768/1‑141 (MQ=255)

                         gtgtGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAActtt                                                                                                               >  2:844951/1‑141 (MQ=255)

                          tgtgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAActttc                                                                                                              <  2:367557/141‑1 (MQ=255)

                           gtgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAActttct                                                                                                             >  1:856999/1‑141 (MQ=255)

                            tgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAActttctt                                                                                                            >  4:351465/1‑141 (MQ=255)

                            tgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAActttctt                                                                                                            >  2:362606/1‑141 (MQ=255)

                               ctaTGTATCTTTATCCATAGTAACGATGGCCGTTAAAGCGACCATGCCATACGTATCTATTTCTTGGAACTTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTaa                                                                                                         <  2:402990/139‑1 (MQ=255)

                               ccATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTaa                                                                                                         >  4:17088/1‑141 (MQ=255)

                                   gTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCACAACCGGCAATTGTT‑AAAACTTTCTTTAAAGcc                                                                                                     >  1:595238/1‑141 (MQ=255)

                                     aTCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGCACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCATAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCAt                                                                                                   >  2:711927/1‑141 (MQ=255)

                                     aTCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCAt                                                                                                   <  2:263181/141‑1 (MQ=255)

                                     aTCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCAt                                                                                                   >  2:804317/1‑141 (MQ=255)

                                      tCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTNTCTAATTNTTGGNACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACNTTCTTTAAANCCAtt                                                                                                  >  4:288226/1‑141 (MQ=255)

                                        tttATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGa                                                                                                <  1:863898/141‑1 (MQ=255)

                                        tttATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGa                                                                                                >  2:157794/1‑141 (MQ=255)

                                         ttATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGGACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAg                                                                                               >  4:101038/1‑141 (MQ=255)

                                          tATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGc                                                                                              <  4:219666/141‑1 (MQ=255)

                                          tATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCATGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTACAGCCATTGAGc                                                                                              >  4:166215/1‑141 (MQ=255)

                                            tCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCtt                                                                                            >  3:175762/1‑141 (MQ=255)

                                              cATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCa                                                                                          >  1:344488/1‑141 (MQ=255)

                                               aTAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCAc                                                                                         <  4:115681/141‑1 (MQ=255)

                                                   tAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCt                                                                                     >  3:40283/1‑141 (MQ=255)

                                                      cGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACa                                                                                  >  2:109061/1‑141 (MQ=255)

                                                        aTGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACata                                                                                >  1:129248/1‑141 (MQ=255)

                                                         tGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGGTTTCAAATCTGCTTGCATACCTGAGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataa                                                                               >  2:296230/1‑141 (MQ=255)

                                                           gCAGTTAAAGCGACCATGCCATGCGAAACGAATACTTGGAACGCCATCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGCACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCAATGAGCTTCACTCCTACataata                                                                             >  4:303129/1‑141 (MQ=255)

                                                            cAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataa                                                                            >  2:27324/1‑141 (MQ=255)

                                                            cAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataa                                                                            >  1:832602/1‑141 (MQ=255)

                                                               ttAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAatat                                                                         <  2:597160/141‑1 (MQ=255)

                                                               ttAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAatat                                                                         >  2:261044/1‑141 (MQ=255)

                                                                tAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATAtt                                                                        <  4:304288/141‑1 (MQ=255)

                                                                tAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATAtt                                                                        <  2:583312/141‑1 (MQ=255)

                                                                tAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATAtt                                                                        <  2:583342/141‑1 (MQ=255)

                                                                   aGCGACCATGCCATACGTATCCAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTa                                                                     >  2:160881/1‑141 (MQ=255)

                                                                     cGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTAtt                                                                   <  2:703008/141‑1 (MQ=255)

                                                                      gACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtc                                                                  <  1:504672/141‑1 (MQ=255)

                                                                      gACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACGTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtc                                                                  >  2:733215/1‑141 (MQ=255)

                                                                        ccATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAGCCGGCAATTGTG‑AAAACTTTCTGTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtcat                                                                >  1:736420/1‑141 (MQ=255)

                                                                        ccATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTaaaaca                                                                >  1:91101/1‑136 (MQ=255)

                                                                        ccATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtcat                                                                >  1:11065/1‑141 (MQ=255)

                                                                        ccATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAAGTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtcat                                                                >  1:86333/1‑141 (MQ=255)

                                                                            gCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcata                                                            <  2:4753/141‑1 (MQ=255)

                                                                              cATACGTATCTAATTCTTGGGACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTATAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatatc                                                          <  2:544948/141‑1 (MQ=255)

                                                                              cATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatatc                                                          <  1:109061/141‑1 (MQ=255)

                                                                              cATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatatc                                                          <  2:12021/141‑1 (MQ=255)

                                                                               aTACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGCACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATCGAGCTTCACTCCTACATAATAATATTGTATTCatcatatca                                                         <  4:183496/141‑1 (MQ=255)

                                                                               aTACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATCGAGCTTCACTCCTACATAATAATATTGTATTCatcatatca                                                         <  4:183510/141‑1 (MQ=255)

                                                                                     aTCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCAttttt                                                    <  2:75572/141‑1 (MQ=255)

                                                                                      tCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTaa                                                  >  2:491782/1‑141 (MQ=255)

                                                                                       cTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAAc                                                 <  2:679943/141‑1 (MQ=255)

                                                                                       cTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAAc                                                 >  1:290595/1‑141 (MQ=255)

                                                                                         aaTTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCt                                               <  4:250827/141‑1 (MQ=255)

                                                                                               tGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGa                                         <  1:469116/141‑1 (MQ=255)

                                                                                               tGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGa                                         <  1:469149/141‑1 (MQ=255)

                                                                                                 gAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGaaa                                       >  4:339180/1‑141 (MQ=255)

                                                                                                 gAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGaaa                                       <  2:212122/141‑1 (MQ=255)

                                                                                                 gAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGaaa                                       <  1:14342/141‑1 (MQ=255)

                                                                                                       tttCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATAtt                                 >  4:104792/1‑141 (MQ=255)

                                                                                                         tCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTaa                               <  4:372161/141‑1 (MQ=255)

                                                                                                               cTGCTTGCATACCTGCTCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTc                         <  1:355985/141‑1 (MQ=255)

                                                                                                                 gCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCaa                       >  1:619242/1‑141 (MQ=255)

                                                                                                                 gCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCaa                       <  2:680942/141‑1 (MQ=255)

                                                                                                                    tGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAatat                    <  1:416308/141‑1 (MQ=255)

                                                                                                                    tGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAatat                    <  1:868243/141‑1 (MQ=255)

                                                                                                                          ccTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTtgatg              <  2:659531/141‑1 (MQ=255)

                                                                                                                           cTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTtgatga             <  1:453676/141‑1 (MQ=255)

                                                                                                                           cTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTtgatga             >  2:625313/1‑141 (MQ=255)

                                                                                                                            tacgcCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTgatgat            >  1:196561/3‑141 (MQ=255)

                                                                                                                            tGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTgatgat            >  2:448278/1‑141 (MQ=255)

                                                                                                                              cgcCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttg          <  1:358508/141‑1 (MQ=255)

                                                                                                                              cgcCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttg          >  2:249265/1‑141 (MQ=255)

                                                                                                                              cgcCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttg          >  1:289998/1‑141 (MQ=255)

                                                                                                                              cgcCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttg          <  2:518804/141‑1 (MQ=255)

                                                                                                                               gcCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgt         <  1:491686/141‑1 (MQ=255)

                                                                                                                                ctatcACTTTTGTCCGGACAGGTACTTGTC‑AAACCTGTCTTTATAGCCCTGGAGCTTCACTCCTGCATAGAAATATTGTATTCATCATATCATTTTTAACCTAATTGTAAAATATTAAGCATTCCATATTATATGAatgtt        >  4:160753/5‑141 (MQ=255)

                                                                                                                                ccAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgtt        >  3:170352/1‑141 (MQ=255)

                                                                                                                                    cACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaaa    >  2:290413/1‑141 (MQ=255)

                                                                                                                                    cACTTGTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaaa    >  1:744246/1‑141 (MQ=255)

                                                                                                                                      cTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATg  <  1:739034/141‑1 (MQ=255)

CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATG  >  USA300TCH1516_ALE/638477‑638751

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