BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	638,614	G→T	5.1%	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	5.1%	703.7 / 17.7	200	T14K (ACA→AAA)	pdxK	Pyridoxine kinase
Reads supporting (aligned to +/- strand): ref base G (96/92); new base T (7/3); total (105/95)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 210 total aligned reads displayed.

ATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAAC  >  USA300TCH1516_ALE/638524‑638704

aTAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTtg                                                                                            <  1:3602312/91‑1 (MQ=255)

  aGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTtgtg                                                                                          >  2:3493536/1‑91 (MQ=255)

    tAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTc                                                                                        <  1:1366354/91‑1 (MQ=255)

    tAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTc                                                                                        <  1:171779/91‑1 (MQ=255)

     aaCGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCa                                                                                       <  2:528021/91‑1 (MQ=255)

     aaCGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCa                                                                                       >  2:2167931/1‑91 (MQ=255)

       cGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGa                                                                                     >  1:3383382/1‑91 (MQ=255)

        gATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGaa                                                                                    >  1:2157382/1‑91 (MQ=255)

          tGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAAcc                                                                                  <  1:105673/91‑1 (MQ=255)

          tGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAAcc                                                                                  <  2:3385016/91‑1 (MQ=255)

           ggCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCg                                                                                 <  2:2898045/91‑1 (MQ=255)

             cAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTATCAGAACCGGc                                                                               >  1:3584497/1‑91 (MQ=255)

             cAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGc                                                                               >  2:2732610/1‑91 (MQ=255)

             cAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGc                                                                               <  1:2995760/91‑1 (MQ=255)

              aGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCa                                                                              <  2:2011076/91‑1 (MQ=255)

              aGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCa                                                                              <  1:89280/91‑1 (MQ=255)

               gTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCaa                                                                             >  2:3371854/1‑91 (MQ=255)

                 tAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAAtt                                                                           >  2:869742/1‑91 (MQ=255)

                 tAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAAtt                                                                           <  2:3747064/91‑1 (MQ=255)

                 tAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAAtt                                                                           >  2:1634578/1‑91 (MQ=255)

                    aGCGACCATGCCATACGTATCTATTTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGtt                                                                        >  2:3128056/1‑91 (MQ=255)

                     gCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTa                                                                       <  1:834755/91‑1 (MQ=255)

                     gCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTa                                                                       >  1:2835641/1‑91 (MQ=255)

                     gCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTa                                                                       >  1:2835786/1‑91 (MQ=255)

                     gCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTa                                                                       <  1:835892/91‑1 (MQ=255)

                      ggACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAGTCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTaa                                                                      <  2:1949797/90‑1 (MQ=255)

                       gACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTaaa                                                                     <  1:1701542/91‑1 (MQ=255)

                        aCCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTaaaa                                                                    >  1:1396149/1‑91 (MQ=255)

                         ccATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAAc                                                                   >  1:3346901/1‑91 (MQ=255)

                         ccATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAAc                                                                   <  2:1352521/91‑1 (MQ=255)

                            tGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActtt                                                                >  1:3840361/1‑91 (MQ=255)

                            tGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAActtt                                                                <  1:2357284/91‑1 (MQ=255)

                             gCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAActttc                                                               >  2:1749133/1‑91 (MQ=255)

                             gCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAActttc                                                               >  1:3606084/1‑91 (MQ=255)

                             gCCATACGTACCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAActttc                                                               >  2:1748701/1‑91 (MQ=255)

                              ccATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAActttct                                                              <  2:2695822/91‑1 (MQ=255)

                               cATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAActttctt                                                             <  2:220704/91‑1 (MQ=255)

                                aTACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAActttcttt                                                            <  2:359437/91‑1 (MQ=255)

                                 tACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTa                                                           >  1:753143/1‑91 (MQ=255)

                                  aCGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTaa                                                          >  2:3261470/1‑91 (MQ=255)

                                   cGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTaaa                                                         <  2:1683977/91‑1 (MQ=255)

                                    gTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAg                                                        >  1:3733163/1‑91 (MQ=255)

                                      aTCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGcc                                                      <  1:3305223/91‑1 (MQ=255)

                                      aTCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGcc                                                      <  1:671101/91‑1 (MQ=255)

                                        cTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCAt                                                    <  2:1510555/91‑1 (MQ=255)

                                        cTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCAt                                                    >  2:3832609/1‑91 (MQ=255)

                                        cTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCAt                                                    <  1:1386446/91‑1 (MQ=255)

                                         tAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCAtt                                                   >  1:1946659/1‑91 (MQ=255)

                                            ttCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCGAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAg                                                <  2:26357/91‑1 (MQ=255)

                                            ttCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAg                                                >  2:1410418/1‑91 (MQ=255)

                                             tCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGc                                               >  1:77543/1‑91 (MQ=255)

                                             tCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGc                                               >  2:3220187/1‑91 (MQ=255)

                                              cTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCt                                              <  1:967596/91‑1 (MQ=255)

                                              cTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCt                                              <  2:2360289/91‑1 (MQ=255)

                                              cTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCt                                              <  2:2549082/91‑1 (MQ=255)

                                              cTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCt                                              <  2:1849679/91‑1 (MQ=255)

                                              cTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCt                                              <  2:2520782/91‑1 (MQ=255)

                                              cTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCt                                              <  2:2489624/91‑1 (MQ=255)

                                              cTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCt                                              >  2:68977/1‑91 (MQ=255)

                                              cTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAACACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCt                                              <  1:3206044/91‑1 (MQ=255)

                                                 ggAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCa                                           >  1:1276425/1‑91 (MQ=255)

                                                   aaCGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTGAAACTTTCTTTAAAGCCATTGAGCTTCACt                                         >  1:2715343/1‑91 (MQ=255)

                                                   aaCGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACt                                         >  1:1535145/1‑91 (MQ=255)

                                                      gTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACTGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCt                                      >  1:3135243/1‑91 (MQ=255)

                                                      gTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCt                                      >  1:3135264/1‑91 (MQ=255)

                                                        tttCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTAc                                    >  1:641830/1‑91 (MQ=255)

                                                        tttCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTAc                                    >  1:3469860/1‑91 (MQ=255)

                                                           cAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACata                                 >  1:2703380/1‑91 (MQ=255)

                                                            aaaTCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataa                                >  2:3275760/1‑91 (MQ=255)

                                                              aTCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataata                              <  1:1040711/91‑1 (MQ=255)

                                                              aTCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataata                              <  1:1992717/91‑1 (MQ=255)

                                                               tCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataa                             <  1:2208119/91‑1 (MQ=255)

                                                                 tGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataata                           >  2:2207587/1‑91 (MQ=255)

                                                                   cTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATAtt                         >  2:722799/1‑91 (MQ=255)

                                                                    ttGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTg                        <  1:3105343/91‑1 (MQ=255)

                                                                    ttGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTg                        <  1:222228/91‑1 (MQ=255)

                                                                    ttGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTg                        <  1:3569205/91‑1 (MQ=255)

                                                                      gCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTa                      <  1:1418993/91‑1 (MQ=255)

                                                                      gCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTa                      >  1:288818/1‑91 (MQ=255)

                                                                      gCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTa                      >  2:3700545/1‑91 (MQ=255)

                                                                           ccTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtcat                 >  1:241403/1‑91 (MQ=255)

                                                                           ccTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtcat                 >  2:302985/1‑91 (MQ=255)

                                                                           ccTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtcat                 >  2:393535/1‑91 (MQ=255)

                                                                             tGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtcatca               <  1:3278863/91‑1 (MQ=255)

                                                                              gcgcCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTcatcat              >  1:1662556/1‑91 (MQ=255)

                                                                               cgcCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCCCTCCTACATAATAATATTGTATTCatcata             >  1:2464241/1‑91 (MQ=255)

                                                                               cgcCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcata             <  1:452981/91‑1 (MQ=255)

                                                                                 ccAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTATATAATAATATTGTATTCatcatatc           <  1:1492768/91‑1 (MQ=255)

                                                                                 ccAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatatc           <  1:3085735/91‑1 (MQ=255)

                                                                                 ccAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatatc           <  1:2502484/91‑1 (MQ=255)

                                                                                     cACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCAttt       >  2:583200/1‑91 (MQ=255)

                                                                                     cACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCAttt       >  1:1164580/1‑91 (MQ=255)

                                                                                      aCTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCAtttt      >  2:1196280/1‑91 (MQ=255)

                                                                                      aCTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCAtttt      <  2:107074/91‑1 (MQ=255)

                                                                                       cTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCAttttt     >  1:2927636/1‑91 (MQ=255)

                                                                                       cTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCAttttt     <  1:486903/91‑1 (MQ=255)

                                                                                         tgtgTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTaa   >  1:2397629/1‑91 (MQ=255)

                                                                                         tgtgTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTaa   >  2:2146129/1‑91 (MQ=255)

                                                                                         tgtgTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTaa   >  2:2595312/1‑91 (MQ=255)

                                                                                          gtgtCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAAc  >  2:3805090/1‑91 (MQ=255)

ATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAAC  >  USA300TCH1516_ALE/638524‑638704

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