BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,172,088	Δ1 bp	8.7%	coding (137/1047 nt)	leuB →	3‑isopropylmalate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,172,088	0	C	.	8.7%	575.2 / 46.9	173	coding (137/1047 nt)	leuB	3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand): ref base C (83/75); new base . (11/4); total (94/79)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

AAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTG  >  USA300TCH1516_ALE/2171998‑2172177

aaaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGc                                                                                            <  2:2928994/91‑1 (MQ=255)

 aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTCGTGGTGc                                                                                            <  1:175118/91‑1 (MQ=37)

 aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGcc                                                                                           <  2:3279592/91‑1 (MQ=255)

 aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGcc                                                                                           >  1:2939986/1‑91 (MQ=255)

 aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGcc                                                                                           >  2:193617/1‑91 (MQ=255)

  aTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCct                                                                                          <  2:1474128/91‑1 (MQ=255)

   ttttGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCctc                                                                                         >  2:2569256/1‑91 (MQ=255)

   ttttGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCctc                                                                                         >  1:3378616/1‑91 (MQ=255)

      tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑t                                                                                     <  2:1116185/91‑1 (MQ=255)

      tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑t                                                                                     >  1:1177521/1‑91 (MQ=255)

       gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTCGTGGTG‑CTCTA‑tt                                                                                    <  2:3280407/91‑1 (MQ=25)

       gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑tt                                                                                    >  1:2581322/1‑91 (MQ=255)

        aaCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTg                                                                                   >  1:3508861/1‑91 (MQ=255)

        aaCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTg                                                                                   >  2:2197157/1‑91 (MQ=255)

         aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGa                                                                                  <  1:53549/91‑1 (MQ=255)

         aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGa                                                                                  >  2:2717917/1‑91 (MQ=255)

         aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTCGTGGTGCCTCTA‑TTGa                                                                                  >  2:3784080/1‑91 (MQ=255)

         aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATAGAACTTTAATTATCAAATAGAGGCACCACCAAATTCGTGGTG‑CTCTA‑tttg                                                                                  <  2:232632/91‑3 (MQ=21)

         aCGGATCTCTATCACTGCTTGAAGCTATAAGTAATAAATATAAGTTTAATTATCGAATAGA‑GCACCACCAAATTCGTGGTG‑CTCTA‑tttga                                                                                 <  2:3677822/91‑4 (MQ=16)

            gATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATAc                                                                               >  1:125484/1‑91 (MQ=255)

             aTCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACa                                                                              >  2:1359630/1‑91 (MQ=255)

              tctcTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACAt                                                                             >  1:3860170/1‑91 (MQ=255)

              tctcTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACAt                                                                             <  1:996770/91‑1 (MQ=255)

                 ctATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCg                                                                          <  1:801916/91‑1 (MQ=255)

                  tATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCgg                                                                         >  1:2523004/1‑91 (MQ=255)

                    tCATTGCTTGAAATTATAAGTAATAAATATATCTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCg                                                                       <  2:3051720/91‑1 (MQ=255)

                      aTTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAg                                                                     >  2:3387012/1‑91 (MQ=255)

                       ttGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGc                                                                    >  2:2311315/1‑91 (MQ=255)

                       ttGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGc                                                                    >  2:2309741/1‑91 (MQ=255)

                        tGCTTGAAATTATAAGTAATAACTATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGcc                                                                   >  2:1102030/1‑91 (MQ=255)

                         gCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCt                                                                  <  1:2577048/91‑1 (MQ=255)

                          cTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCtt                                                                 <  1:1489028/91‑1 (MQ=255)

                            tGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTa                                                               >  1:2305761/1‑91 (MQ=255)

                             gAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTaa                                                              <  1:461694/91‑1 (MQ=255)

                              aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAAc                                                             <  1:1441979/91‑1 (MQ=255)

                              aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAAc                                                             <  1:925849/91‑1 (MQ=255)

                              aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAAc                                                             >  2:1749700/1‑91 (MQ=255)

                               aaTTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACt                                                            <  2:862877/91‑1 (MQ=255)

                                 ttATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTga                                                          <  1:1657558/91‑1 (MQ=255)

                                     aaGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGaaa                                                      >  1:2160224/1‑91 (MQ=255)

                                      aGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGcaaa                                                     >  2:3848706/1‑91 (MQ=255)

                                      aGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGaaaa                                                     >  1:2492590/1‑91 (MQ=255)

                                         aataaataTAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCt                                                  <  1:2216033/91‑1 (MQ=255)

                                          ataaataTAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑GTGATACATTCGGCGAGCCTTTAACTGAGAAAACCtt                                                 >  2:3123215/1‑91 (MQ=255)

                                          ataaataTAACTTTAATTATCAAATAGA‑GCACCACGAACTTGGTGGTGCCTCTA‑TTGATACATTCGGCCAGCCTCTAACTGCGAAAACCtt                                                 >  2:243025/1‑91 (MQ=255)

                                            aaataTAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTaa                                               <  1:1060265/91‑1 (MQ=255)

                                               tataACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATg                                            <  2:621729/91‑1 (MQ=255)

                                                ataACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGc                                           <  1:1093977/91‑1 (MQ=255)

                                                 taACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCg                                          <  1:1993048/91‑1 (MQ=255)

                                                   aCTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGAGAGCCTCTAACTGAGAAAACCTTAAATGCgtg                                        >  2:1623694/1‑91 (MQ=255)

                                                    cTTTAATTATTAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCgtgt                                       >  2:1055697/1‑91 (MQ=255)

                                                      ttAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaa                                     >  1:1330236/1‑91 (MQ=255)

                                                       tAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaaa                                    >  1:3923710/1‑91 (MQ=255)

                                                        aaTTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaaaa                                   <  1:910506/91‑1 (MQ=255)

                                                         aTTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAag                                  <  1:362775/91‑1 (MQ=255)

                                                             tCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCa                              <  2:3894000/91‑1 (MQ=255)

                                                             atcaaTAGAGGCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGc                               >  1:3751531/4‑91 (MQ=17)

                                                              tcaaTAGAGGCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCa                              >  1:711013/3‑91 (MQ=17)

                                                              cAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGATCAg                             >  1:3194616/1‑91 (MQ=255)

                                                              cAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAg                             >  1:871095/1‑91 (MQ=255)

                                                              cAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAg                             >  1:49423/1‑91 (MQ=255)

                                                              cAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGACCAg                             >  1:3194800/1‑91 (MQ=255)

                                                               caaTAGAGGCACCACCAAATTCGAGGTG‑CTCTATTTGATACATCCGGCGAGCCTTTAACTGAGAAACCCTTAAATGCCTGTAAAAGAGCAg                             >  2:2145209/2‑91 (MQ=9)

                                                               aaaTAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGa                            >  1:2806684/1‑91 (MQ=255)

                                                                 aTAGAGGCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGAt                           >  1:231436/1‑91 (MQ=18)

                                                                 aTAGAGGCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGAt                           >  2:3107851/1‑91 (MQ=18)

                                                                 aTAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAATAGCatctg                          <  1:571683/91‑5 (MQ=255)

                                                                 aTAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATg                          >  1:3905014/1‑91 (MQ=255)

                                                                  tAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGc                         <  2:3771425/91‑1 (MQ=255)

                                                                  tAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGc                         >  1:2388152/1‑91 (MQ=255)

                                                                   aga‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAACCCCTTAAATGCGTGTAAAATAGCATCTGCt                        <  2:2651595/91‑1 (MQ=255)

                                                                     aggCACCACCAAATTCGTGGTG‑CTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAACCCCTTAAATGCGTGTAAAAGCGCCGATGCTAt                      >  2:3885398/3‑91 (MQ=21)

                                                                      ggCACCACCAAATTCGTGGTG‑CTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt                     >  1:1842209/2‑91 (MQ=37)

                                                                       gCACCACGAATTTGGTGGTGCCTCTA‑TTGCTACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt                     >  2:3446694/1‑91 (MQ=255)

                                                                       gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt                     <  1:3553625/91‑1 (MQ=255)

                                                                       gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt                     <  2:626659/91‑1 (MQ=255)

                                                                         accacGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtttt                   >  2:2972482/1‑91 (MQ=255)

                                                                         accacGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtttt                   <  1:1669924/91‑1 (MQ=255)

                                                                          ccacGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTa                  >  2:3024/1‑91 (MQ=255)

                                                                          ccacGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTa                  >  2:1097410/1‑91 (MQ=255)

                                                                          ccacGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTa                  >  1:2289436/1‑91 (MQ=255)

                                                                          ccacGAATTGGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTa                  >  1:1977069/1‑91 (MQ=255)

                                                                           cacGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTAc                 >  1:3548333/1‑91 (MQ=255)

                                                                            acGAATTTGGTGGTGCCTCTA‑TTGATACGTTCGGCGAGCCTTTAACTGAGATATCCTTATATGCGTGTAAAAGTGCTGATGCTATTTTACt                <  2:417599/91‑1 (MQ=255)

                                                                               aaTTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTggg             <  1:3619515/91‑1 (MQ=255)

                                                                                aTTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTGTTTTACTGGGt            >  1:3578376/1‑91 (MQ=255)

                                                                                aTTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGt            <  1:2868693/91‑1 (MQ=255)

                                                                                aTTGGGTGGTGCCTCTA‑TTTATGCATGCGGCGAGCCTTTAGCTGAGAAAGCCGTAGGTGCGTGTAACAGCGCAGATGCTATTTTACTGGGt            <  2:922219/91‑1 (MQ=255)

                                                                                 tttGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTg           <  2:1590935/91‑1 (MQ=255)

                                                                                 tttGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTg           <  2:1218261/91‑1 (MQ=255)

                                                                                 tttGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTg           <  2:1218235/91‑1 (MQ=255)

                                                                                  ttGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGc          <  1:1955270/91‑1 (MQ=255)

                                                                                    ggtggtGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAATAGCAGATGCTATTTTACTGGGTGCaa        <  2:2250906/91‑1 (MQ=255)

                                                                                    ggtggtGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCaa        <  1:987840/91‑1 (MQ=255)

                                                                                    ggtggtGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCaa        >  2:2020686/1‑91 (MQ=255)

                                                                                    ggtggtGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCaa        <  1:1312072/91‑1 (MQ=255)

                                                                                    ggtggtGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCaa        >  2:1715869/1‑91 (MQ=255)

                                                                                      tggtgCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATc      >  1:1269432/1‑91 (MQ=255)

                                                                                         tgCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGt   >  2:2794836/1‑91 (MQ=255)

                                                                                          gCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTg  <  1:127834/91‑1 (MQ=255)

AAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTG  >  USA300TCH1516_ALE/2171998‑2172177

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