BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	756,705	T→A	6.4%	intergenic (‑193/‑20)	uppP ← / → cydD	Undecaprenyl‑diphosphatase/ATP‑binding/permease protein CydD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	756,705	0	T	A	6.4%	482.9 / 15.2	141	intergenic (‑193/‑20)	uppP/cydD	Undecaprenyl‑diphosphatase/ATP‑binding/permease protein CydD
Reads supporting (aligned to +/- strand): ref base T (70/61); new base A (7/2); total (78/63)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 152 total aligned reads displayed.

TAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCTTGGTCAGTACTGG  >  USA300TCH1516_ALE/756615‑756795

tAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTataat                                                                                            >  1:470043/1‑91 (MQ=255)

 aaaGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTataata                                                                                           <  1:1945470/91‑1 (MQ=255)

    gATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGt                                                                                        >  1:3417946/1‑91 (MQ=255)

    gATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGt                                                                                        >  1:3417911/1‑91 (MQ=255)

    gATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGt                                                                                        >  1:1504644/1‑91 (MQ=255)

     aTGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGtt                                                                                       <  1:1169445/91‑1 (MQ=255)

      tGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTATTATAAAAGGttt                                                                                      <  1:1823287/91‑1 (MQ=255)

           atatTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTg                                                                                 <  1:699089/91‑1 (MQ=255)

           atatTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTg                                                                                 <  1:698967/91‑1 (MQ=255)

            tatTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGa                                                                                <  1:3745890/91‑1 (MQ=255)

            tatTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGa                                                                                <  1:3745922/91‑1 (MQ=255)

            tatTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGa                                                                                <  1:3183197/91‑1 (MQ=255)

                cATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGAtata                                                                            >  1:2800657/1‑91 (MQ=255)

                cATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGAtata                                                                            <  1:1379183/91‑1 (MQ=255)

                 aTTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATaa                                                                           >  1:930734/1‑91 (MQ=255)

                 aTTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATaa                                                                           <  1:177046/91‑1 (MQ=255)

                 aTTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATaa                                                                           >  1:887932/1‑91 (MQ=255)

                  ttttCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAAt                                                                          <  1:576636/91‑1 (MQ=255)

                  ttttCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAAt                                                                          <  1:1226919/91‑1 (MQ=255)

                  ttttCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAAt                                                                          >  1:182992/1‑91 (MQ=255)

                   tttCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAAtt                                                                         >  1:3044813/1‑91 (MQ=255)

                    ttCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATtg                                                                        <  1:1829899/91‑1 (MQ=255)

                      cTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATtgtg                                                                      >  1:1418007/1‑91 (MQ=255)

                       ttttACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGGTTTTAAGGTAAATCTTTTATAATAGGTTTAGGTGATATAATTGTGa                                                                     <  1:3315906/91‑1 (MQ=255)

                       ttttACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGa                                                                     <  1:3915675/91‑1 (MQ=255)

                        tttACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTCGGTGATATAATTGTGaa                                                                    >  1:3550518/1‑91 (MQ=255)

                          tACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGaaaa                                                                  <  1:622941/91‑1 (MQ=255)

                          tACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGaaaa                                                                  <  1:2712876/91‑1 (MQ=255)

                          tACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGaaaa                                                                  <  1:3276095/91‑1 (MQ=255)

                           aCATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGaaaaa                                                                 >  1:2060750/1‑91 (MQ=255)

                           aCATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGaaaaa                                                                 <  1:3995433/91‑1 (MQ=255)

                           aCATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGaaaaa                                                                 >  1:329881/1‑91 (MQ=255)

                            cATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGaaaaaa                                                                >  1:2941277/1‑91 (MQ=255)

                            cATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGaaaaaa                                                                >  1:2941451/1‑91 (MQ=255)

                             aTGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAAt                                                               <  1:1148026/91‑1 (MQ=255)

                               gTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTa                                                             >  1:347250/1‑91 (MQ=255)

                               gTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTa                                                             <  1:2354232/91‑1 (MQ=255)

                               gTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTa                                                             >  1:4126079/1‑91 (MQ=255)

                                 aaaaCTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTaac                                                           >  1:3517386/1‑91 (MQ=255)

                                  aaaCTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTaaca                                                          >  1:1748644/1‑91 (MQ=255)

                                  aaaCTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTaaca                                                          >  1:2558882/1‑91 (MQ=255)

                                     cTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAacaaca                                                       >  1:3992004/1‑91 (MQ=255)

                                      tACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAAcaacaa                                                      <  1:400769/91‑1 (MQ=255)

                                      tACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAAcaacaa                                                      >  1:3069038/1‑91 (MQ=255)

                                      tACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAAcaacaa                                                      >  1:34427/1‑91 (MQ=255)

                                       aCAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAAt                                                     >  1:1051978/1‑91 (MQ=255)

                                         aaTATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATAc                                                   >  1:342430/1‑91 (MQ=255)

                                         aaTATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATAc                                                   >  1:342330/1‑91 (MQ=255)

                                          atatTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACt                                                  >  1:1671199/1‑91 (MQ=255)

                                            atTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGt                                                >  1:3427212/1‑91 (MQ=255)

                                             ttttATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGtt                                               <  1:2722944/91‑1 (MQ=255)

                                             ttttATTGAGCTATCTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATCTAATTGTGAAAAAATTAACAACAATACTGtt                                               <  1:3741670/91‑1 (MQ=255)

                                                    gAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAatat                                        >  1:276935/1‑91 (MQ=255)

                                                    gAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAatat                                        <  1:2446162/91‑1 (MQ=255)

                                                     aGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAAtata                                       <  1:3698421/91‑1 (MQ=255)

                                                     aGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAAtata                                       <  1:3482057/91‑1 (MQ=255)

                                                     aGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAAtata                                       <  1:3749959/91‑1 (MQ=255)

                                                     aGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAAtata                                       <  1:3749898/91‑1 (MQ=255)

                                                      gCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATaa                                      <  1:1038463/91‑1 (MQ=255)

                                                      gCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATaa                                      >  1:1310600/1‑91 (MQ=255)

                                                             aaTTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAAttttt                               <  1:3908468/91‑1 (MQ=255)

                                                              aTTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTc                              >  1:482098/1‑91 (MQ=255)

                                                               tttGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTcc                             >  1:1668774/1‑91 (MQ=255)

                                                               tttGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTcc                             <  1:3183943/91‑1 (MQ=255)

                                                               tttGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTcc                             >  1:2860013/1‑91 (MQ=255)

                                                                ttGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCg                            <  1:466856/91‑1 (MQ=255)

                                                                ttGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCg                            <  1:3021263/91‑1 (MQ=255)

                                                                ttGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCg                            >  1:2089602/1‑91 (MQ=255)

                                                                 tGATTTTAAGGAAAACCTTTTATACTAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCgg                           <  1:1105378/91‑1 (MQ=255)

                                                                   aTTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTa                         >  1:3654047/1‑91 (MQ=255)

                                                                    ttttAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTAc                        >  1:2900047/1‑91 (MQ=255)

                                                                     ttttAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGATACt                       <  1:3269616/91‑1 (MQ=255)

                                                                     tttAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACt                       <  1:847112/91‑1 (MQ=255)

                                                                     tttAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACt                       <  1:1442520/91‑1 (MQ=255)

                                                                      ttAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTc                      >  1:1622094/1‑91 (MQ=255)

                                                                      ttAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTc                      <  1:1249452/91‑1 (MQ=255)

                                                                      ttAAGGAAAACCTTTTATAACAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTc                      >  1:1069556/1‑91 (MQ=255)

                                                                         aGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATg                   >  1:1762660/1‑91 (MQ=255)

                                                                         aGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATg                   <  1:2678738/91‑1 (MQ=255)

                                                                          ggAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGt                  <  1:1815775/91‑1 (MQ=255)

                                                                          ggAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGt                  <  1:1552913/91‑1 (MQ=255)

                                                                             aaaCCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCt               >  1:1417728/1‑91 (MQ=255)

                                                                             aaaCCTATTATAAAAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCt               >  1:3565902/1‑91 (MQ=255)

                                                                              aaCCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCtt              >  1:3133625/1‑91 (MQ=255)

                                                                              aaCCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCtt              <  1:566682/91‑1 (MQ=255)

                                                                              aaCCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCtt              >  1:3139255/1‑91 (MQ=255)

                                                                               aCCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCTTg             >  1:1337232/1‑91 (MQ=255)

                                                                                ccTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCTTgg            <  1:2200556/91‑1 (MQ=255)

                                                                                ccTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCTTgg            >  1:64681/1‑91 (MQ=255)

                                                                                ccTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCTTgg            >  1:3518262/1‑91 (MQ=255)

                                                                                 cTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCTTGGt           >  1:3617553/1‑91 (MQ=255)

                                                                                 cTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCTTGGt           >  1:3617332/1‑91 (MQ=255)

                                                                                   attATAAAAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCTTGGTCa         >  1:2719604/2‑91 (MQ=255)

                                                                                   attATAAAAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCTTGGTCa         >  1:2645827/2‑91 (MQ=255)

                                                                                    ttATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCTTGGTCAg        >  1:1420773/1‑91 (MQ=255)

                                                                                     tataATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCTTGGTCAGt       <  1:3677709/91‑1 (MQ=255)

                                                                                     tataATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCTTGGTCAGt       >  1:4066195/1‑91 (MQ=255)

                                                                                     tataATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCTTGGTCAGt       >  1:3714998/1‑91 (MQ=255)

                                                                                      ataataGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCTTGGTCAGTa      <  1:2475217/91‑1 (MQ=255)

                                                                                       taataGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCTTGGTCAGTAc     >  1:2507107/1‑91 (MQ=255)

TAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAACAACAATACTGTTTCAATATAAAATTTTTCCGGTACTCATGTTCTTGGTCAGTACTGG  >  USA300TCH1516_ALE/756615‑756795

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