BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,668,300	G→T	5.1%	intergenic (‑109/‑442)	USA300HOU_RS13735 ← / → USA300TCH1516_02539	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,668,300	0	G	T	5.1%	819.8 / 24.7	233	intergenic (‑109/‑442)	USA300HOU_RS13735/USA300TCH1516_02539	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (108/113); new base T (8/4); total (116/117)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 241 total aligned reads displayed.

TATGTTTATTTTATTTTCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTAACTTTTTACT  >  USA300TCH1516_ALE/2668210‑2668390

tATGTTTATTTTATTTTCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTg                                                                                            <  2:204931/91‑1 (MQ=255)

tATGTTTATTTTATTTTCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTg                                                                                            <  2:2279252/91‑1 (MQ=255)

 aTGTTTATTTTATTTTCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGa                                                                                           >  1:1267009/1‑91 (MQ=255)

  tGTTTATTTTATTTTCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGaa                                                                                          >  2:1677666/1‑91 (MQ=255)

      tattttatttTCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAatat                                                                                      <  2:130686/91‑1 (MQ=255)

        ttttattttCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGGAatattt                                                                                    >  2:1697219/1‑91 (MQ=255)

        ttttattttCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAatattt                                                                                    <  2:320278/91‑1 (MQ=255)

         tttattttCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAatatttg                                                                                   >  1:13555/1‑91 (MQ=255)

           tattttCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAatatttgat                                                                                 >  1:1523420/1‑91 (MQ=255)

            attttCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAatatttgata                                                                                >  1:545617/1‑91 (MQ=255)

            attttCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAatatttgata                                                                                <  2:894224/91‑1 (MQ=255)

             ttttCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAatatttgatat                                                                               >  2:565239/1‑91 (MQ=255)

             ttttCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAatatttgatat                                                                               <  2:77276/91‑1 (MQ=255)

              tttCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAatatttgatatt                                                                              >  2:349791/1‑91 (MQ=255)

               ttCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAatatttgatattt                                                                             >  2:1415754/1‑91 (MQ=255)

                tCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAatatttgatatttg                                                                            >  1:2799568/1‑91 (MQ=255)

                tCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAatatttgatatttg                                                                            >  2:2903628/1‑91 (MQ=255)

                tCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAatatttgatatttg                                                                            <  1:257386/91‑1 (MQ=255)

                  aGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTtgtg                                                                          <  2:1019540/91‑1 (MQ=255)

                  aGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTtgtg                                                                          <  2:2482766/91‑1 (MQ=255)

                  aGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTtgtg                                                                          <  1:2682122/91‑1 (MQ=255)

                   gACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGa                                                                         >  1:813806/1‑91 (MQ=255)

                   gACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGa                                                                         <  1:2373628/91‑1 (MQ=255)

                   gACCATAAACATCGTACACCCTCAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGa                                                                         <  2:3119260/91‑1 (MQ=255)

                     ccATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGaaa                                                                       >  1:1323569/1‑91 (MQ=255)

                       aTAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAACTATTTGATATTTGTGAAAGa                                                                     >  1:1001151/1‑91 (MQ=255)

                        tAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGAt                                                                    >  1:1685616/1‑91 (MQ=255)

                        tAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGAt                                                                    <  2:577530/91‑1 (MQ=255)

                           aCATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGAtttt                                                                 <  1:2971470/91‑1 (MQ=255)

                           aCATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGAtttt                                                                 <  1:1132236/91‑1 (MQ=255)

                           aCATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGAtttt                                                                 <  1:776339/91‑1 (MQ=255)

                             aTCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTa                                                               >  1:674994/1‑91 (MQ=255)

                               cGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTaaa                                                             >  1:1932328/1‑91 (MQ=255)

                                gcacacCCTTAAGAAAACGTTTTCGACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTaaaa                                                            <  1:3720163/89‑1 (MQ=255)

                                gTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTaaaa                                                            <  1:1115386/91‑1 (MQ=255)

                                  acacCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATg                                                          <  2:1954668/91‑1 (MQ=255)

                                     cccTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGAt                                                       <  1:3274127/91‑1 (MQ=255)

                                     cccTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGAt                                                       <  1:2232151/91‑1 (MQ=255)

                                      ccTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGAtt                                                      >  2:1319721/1‑91 (MQ=255)

                                        ttAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGt                                                    <  1:2095661/91‑1 (MQ=255)

                                        ttAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGt                                                    <  1:1964255/91‑1 (MQ=255)

                                        ttAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGt                                                    >  1:3761517/1‑91 (MQ=255)

                                          aaGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGttt                                                  >  1:403883/1‑91 (MQ=255)

                                          aaGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGttt                                                  <  1:897917/91‑1 (MQ=255)

                                           aGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTc                                                 >  1:3465407/1‑91 (MQ=255)

                                            gAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGtt                                                   >  2:1142521/1‑88 (MQ=255)

                                             aaaaCGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCaa                                               <  1:1570207/91‑1 (MQ=255)

                                               aaCGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCaaaa                                             <  1:3248603/91‑1 (MQ=255)

                                                aCGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCaaaaa                                            >  1:216497/1‑91 (MQ=255)

                                                    tttCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAactac                                        >  2:1568060/1‑91 (MQ=255)

                                                     ttCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGttt                                                  >  1:3085703/1‑80 (MQ=255)

                                                     ttCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAactact                                       >  2:1934280/1‑91 (MQ=255)

                                                     ttCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAactact                                       >  1:1357516/1‑91 (MQ=255)

                                                     ttCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAactact                                       >  2:3746342/1‑91 (MQ=255)

                                                      tCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACtt                                      >  1:844491/1‑91 (MQ=255)

                                                       cAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTa                                     <  2:1060473/91‑1 (MQ=255)

                                                       cAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTa                                     <  2:1857996/91‑1 (MQ=255)

                                                        aaCTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTAAAAAACTACTTAt                                    <  2:3164612/91‑1 (MQ=255)

                                                           ttttATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTAtgtt                                 <  1:514196/91‑1 (MQ=255)

                                                           ttttATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTAtgtt                                 >  1:22811/1‑91 (MQ=255)

                                                                tCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGt                            >  1:920470/1‑91 (MQ=255)

                                                                tCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGt                            >  1:2537107/1‑91 (MQ=255)

                                                                 cTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGtt                           >  2:2740273/1‑91 (MQ=255)

                                                                 cTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGtt                           <  1:3581349/91‑1 (MQ=255)

                                                                 cTGTTATCAAATCAAATATTTAAGTGAAATACTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTACAAAAACTACTCATGTTGTCGtt                           >  2:2959024/1‑91 (MQ=255)

                                                                  tGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTa                          <  2:680872/91‑1 (MQ=255)

                                                                   gTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTaa                         <  1:2063309/91‑1 (MQ=255)

                                                                    ttATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAat                        <  1:1889976/91‑1 (MQ=255)

                                                                    ttATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAat                        >  1:2559176/1‑91 (MQ=255)

                                                                     tATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAata                       <  2:3059086/91‑1 (MQ=255)

                                                                      atcaaatcaaATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAatat                      <  2:1526194/91‑1 (MQ=255)

                                                                        aatatcaaatATTTCACTTAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATAttt                    >  1:3321534/4‑91 (MQ=37)

                                                                          aatcaaatATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTAc                  <  2:3860852/91‑1 (MQ=255)

                                                                          aatcaaatATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTAc                  >  1:384925/1‑91 (MQ=255)

                                                                            tcaaatATTTCACTTAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTa                >  2:3806910/1‑91 (MQ=38)

                                                                            tcaaatATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTa                >  1:1314744/1‑91 (MQ=255)

                                                                            tcaaatATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTa                >  2:2414247/1‑91 (MQ=255)

                                                                            tcaaatATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTa                >  2:2153470/1‑91 (MQ=255)

                                                                            tcaaatATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTa                >  2:1722463/1‑91 (MQ=255)

                                                                             caaatATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTaa               <  2:3248795/91‑1 (MQ=255)

                                                                             caaatATTGAAGTGACATATTTGATCTTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTaa               <  2:3367209/91‑1 (MQ=255)

                                                                              aaatATTTCACTTAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAAt              >  2:2563634/1‑91 (MQ=39)

                                                                              aaatATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAAt              >  1:570921/1‑91 (MQ=255)

                                                                                atatTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTa            >  1:599950/1‑91 (MQ=255)

                                                                                 tatTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTaa           >  2:1857043/1‑91 (MQ=255)

                                                                                 tatTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTaa           <  1:1447129/91‑1 (MQ=255)

                                                                                  atTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTAAc          >  2:1378040/1‑91 (MQ=255)

                                                                                  atTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTAAc          <  1:3749244/91‑1 (MQ=255)

                                                                                  atTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTAAc          >  2:71925/1‑91 (MQ=255)

                                                                                  atTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTAAc          >  1:3654623/1‑91 (MQ=255)

                                                                                   tttAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTAACt         >  2:3879323/1‑91 (MQ=255)

                                                                                   tttAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTAACt         >  1:1397495/1‑91 (MQ=255)

                                                                                    ttAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTAACtt        <  2:2891219/91‑1 (MQ=255)

                                                                                     tAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTAACttt       <  1:3031560/91‑1 (MQ=255)

                                                                                     tAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTAACttt       <  1:3538625/91‑1 (MQ=255)

                                                                                      aaGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTAACtttt      <  2:3645736/91‑1 (MQ=255)

                                                                                        gTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTAACTTTTTa    <  2:815539/91‑1 (MQ=255)

                                                                                          gAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTTATTAACTTTTTAct  <  2:3292665/91‑1 (MQ=255)

                                                                                          gAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTAACTTTTTAct  <  2:3549528/91‑1 (MQ=255)

                                                                                          gAAATATTTGATATTTGTGAAAGATTTTTAAAAAAAAAAATAAAAAAAACTACTTATGTTGTCGTTAATATTTACTAATTAACTTTTTAct  <  2:2664370/91‑1 (MQ=255)

TATGTTTATTTTATTTTCAGACCATAAACATCGTACACCCTTAAGAAAACGTTTTCAACTTTTATCTGTTATCAAATCAAATATTTAAGTGAAATATTTGATATTTGTGAAAGATTTTTAAAATGGATTGTTTCAAAAACTACTTATGTTGTCGTTAATATTTACTAATTAACTTTTTACT  >  USA300TCH1516_ALE/2668210‑2668390

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: `‑`