BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,479,041	A→T	9.3%	K487N (AAA→AAT)	USA300HOU_RS12765 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,479,041	0	A	T	9.3%	656.1 / 42.8	202	K487N (AAA→AAT)	USA300HOU_RS12765	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (89/94); new base T (10/9); total (99/103)
Fisher's exact test for biased strand distribution p-value = 8.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 213 total aligned reads displayed.

TAGATATGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCTAA  >  USA300TCH1516_ALE/2478951‑2479131

tagataTGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCaaa                                                                                            <  2:1865410/91‑1 (MQ=255)

 agataTGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAaat                                                                                           >  2:1424565/1‑91 (MQ=255)

  gataTGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAaata                                                                                          >  1:1697955/1‑91 (MQ=255)

     atGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaatt                                                                                       >  2:3759653/1‑91 (MQ=255)

     atGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaatt                                                                                       <  1:2872334/91‑1 (MQ=255)

      tGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaatta                                                                                      >  2:1197126/1‑91 (MQ=255)

      tGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaatta                                                                                      >  2:468294/1‑91 (MQ=255)

       ggTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaa                                                                                     <  1:1803343/91‑1 (MQ=255)

        gTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                    <  1:1068952/91‑1 (MQ=255)

        gTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATAATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                    >  1:3095890/1‑91 (MQ=255)

         ttttAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATAAAAAAAAAAAAAGCTGATTTCTATCAAATAATTAATa                                                                                   <  2:164996/91‑1 (MQ=21)

           ttAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGa                                                                                 >  1:3040526/1‑91 (MQ=255)

            tAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGaa                                                                                >  1:2553208/1‑91 (MQ=255)

             aaaCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGaaa                                                                               <  2:3637300/91‑1 (MQ=255)

              aaCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                    <  2:842488/85‑1 (MQ=255)

              aaCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                    >  1:842488/1‑85 (MQ=255)

              aaCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAAt                                                                              >  1:1060834/1‑91 (MQ=255)

              aaCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAAt                                                                              >  2:136110/1‑91 (MQ=255)

              aaCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAAt                                                                              >  2:2711356/1‑91 (MQ=255)

               aCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAATc                                                                             <  2:1010674/91‑1 (MQ=18)

               aCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATc                                                                             <  2:122926/91‑1 (MQ=255)

                cGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAATCa                                                                            <  2:1705964/91‑1 (MQ=18)

                cGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCa                                                                            >  2:424681/1‑91 (MQ=255)

                 gTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAg                                                                           >  2:1774074/1‑91 (MQ=255)

                 gTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAg                                                                           <  1:706841/91‑1 (MQ=255)

                 gTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAg                                                                           >  1:1853848/1‑91 (MQ=255)

                  tGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGc                                                                          >  1:1483230/1‑91 (MQ=255)

                  tGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGc                                                                          >  2:767315/1‑91 (MQ=255)

                     cTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttt                                                                       <  1:3760956/91‑1 (MQ=255)

                       tCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttttt                                                                     <  1:2711127/91‑1 (MQ=255)

                       tCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttttt                                                                     >  2:1093483/1‑91 (MQ=255)

                       tCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttttt                                                                     <  1:136110/91‑1 (MQ=255)

                        cAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtttttt                                                                    >  2:3667571/1‑91 (MQ=255)

                        cAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtttttt                                                                    <  2:3720869/91‑1 (MQ=255)

                        cAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtttttt                                                                    <  2:1304049/91‑1 (MQ=255)

                        cAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtttt                                                                      >  2:2258639/1‑89 (MQ=255)

                              atatGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTg                                                              <  2:1527770/91‑1 (MQ=40)

                              atatGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTg                                                              <  1:3880213/91‑1 (MQ=40)

                               tatGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGc                                                             >  1:1926439/1‑91 (MQ=39)

                               tatGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGc                                                             >  2:1934909/1‑91 (MQ=39)

                               tatGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGc                                                             <  1:562312/91‑1 (MQ=39)

                               tatGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGc                                                             >  1:2462682/1‑91 (MQ=39)

                                  gAATGCGACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTa                                                          >  2:1844795/1‑91 (MQ=18)

                                   aaTGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTaa                                                         <  2:2829875/91‑1 (MQ=34)

                                     tGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGa                                                       <  1:3222661/91‑1 (MQ=32)

                                                    atgataAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTGTACATTGCCTAAGAACTTAATGTCCCAGc                                        <  2:292131/88‑1 (MQ=12)

                                                    aTTATAAAATACATAACAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGc                                        <  2:1136193/91‑1 (MQ=17)

                                                     ttATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGcc                                       <  1:116819/91‑1 (MQ=33)

                                                      tataAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCt                                      <  1:1789560/91‑1 (MQ=34)

                                                        taAAATATATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGt                                    <  1:1938756/91‑1 (MQ=37)

                                                         aaaaTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGtt                                   <  2:2346538/91‑1 (MQ=38)

                                                          aaaTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGttt                                  >  1:2526776/1‑91 (MQ=39)

                                                          aaaTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGttt                                  <  1:1885460/91‑1 (MQ=39)

                                                           aaTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGtttt                                 <  2:3864051/91‑1 (MQ=39)

                                                           aaTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGtttt                                 <  2:3198461/91‑1 (MQ=40)

                                                                atatCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTa                            >  2:3009654/1‑91 (MQ=255)

                                                                 tatCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTaa                           <  1:1272158/91‑1 (MQ=255)

                                                                  atCAAAAAAGCTGATTTCTATCAATTAATTAATTGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTaaa                          <  1:3794889/91‑1 (MQ=255)

                                                                  atCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTaaa                          >  1:845779/1‑91 (MQ=255)

                                                                  atCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTaaa                          <  2:1458381/91‑1 (MQ=255)

                                                                   tCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAAc                         >  2:2727259/1‑91 (MQ=255)

                                                                     aaaaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAAccc                       >  2:1871591/1‑91 (MQ=255)

                                                                     aaaaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAAccc                       >  1:1400033/1‑91 (MQ=255)

                                                                      aaaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCa                      <  2:1158034/91‑1 (MQ=255)

                                                                       aaaaGCTGATTTCTATTAATTATTTGATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCaa                     >  1:3353271/1‑91 (MQ=18)

                                                                       aaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCaa                     >  2:2497034/1‑91 (MQ=255)

                                                                        aaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCaaa                    >  1:3260182/1‑91 (MQ=255)

                                                                        aaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCaaa                    <  1:3877172/91‑1 (MQ=255)

                                                                          aGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATc                  >  1:2117872/1‑91 (MQ=255)

                                                                          aGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATc                  <  1:2880261/91‑1 (MQ=255)

                                                                          aGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATc                  <  2:3057986/91‑1 (MQ=255)

                                                                           gCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCg                 <  1:1329139/91‑1 (MQ=255)

                                                                           gCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCg                 >  2:660603/1‑91 (MQ=255)

                                                                           gCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCg                 >  1:2599402/1‑91 (MQ=255)

                                                                           gCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCg                 <  1:1013434/91‑1 (MQ=255)

                                                                             tGATTTCTATTAATTATTTGATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTa               >  2:1041615/1‑91 (MQ=21)

                                                                              gATTTCTATTAATTATTTGATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAt              >  2:3367504/1‑91 (MQ=21)

                                                                              gATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAt              <  2:785046/91‑1 (MQ=255)

                                                                              gATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAt              >  1:2079682/1‑91 (MQ=255)

                                                                              gATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAt              <  1:3667571/91‑1 (MQ=255)

                                                                              gATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAt              >  1:2195508/1‑91 (MQ=255)

                                                                               aTTTCTATTAATTATTTGATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAtt             >  1:3707024/1‑91 (MQ=21)

                                                                               aTTTCTATTAATTATTTGATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAtt             >  2:3247919/1‑91 (MQ=21)

                                                                               aTTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAtt             <  2:1697955/91‑1 (MQ=255)

                                                                                tttCTATTAATTATTTGATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTc            >  1:1798286/1‑91 (MQ=21)

                                                                                tttCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTc            >  2:2559874/1‑91 (MQ=255)

                                                                                tttCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTCAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTc            >  2:918040/1‑91 (MQ=255)

                                                                                  tCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTg          <  2:675889/91‑1 (MQ=255)

                                                                                   cTATCAAATAATTAATAGAAATCAGCTTTTTTGCATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGt         >  2:3816206/1‑91 (MQ=255)

                                                                                    tATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGtt        <  1:2824296/91‑1 (MQ=255)

                                                                                    tATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGtt        <  2:3146522/91‑1 (MQ=255)

                                                                                     aTCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTa       <  1:1873969/91‑1 (MQ=255)

                                                                                     aTCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTa       >  2:702726/1‑91 (MQ=255)

                                                                                      tCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTAt      <  2:1499454/91‑1 (MQ=255)

                                                                                       cAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTc            >  2:2680331/1‑84 (MQ=255)

                                                                                       cAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATc     >  1:2059429/1‑91 (MQ=255)

                                                                                        aaaTAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCt    >  1:2670842/1‑91 (MQ=255)

                                                                                        aaaTAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCt    <  1:1806004/91‑1 (MQ=255)

                                                                                        aaaTAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCt    >  1:1520569/1‑91 (MQ=255)

                                                                                          ataatTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCTaa  <  1:2343885/91‑1 (MQ=255)

TAGATATGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCTAA  >  USA300TCH1516_ALE/2478951‑2479131

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