BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,089,908	T→G	5.2%	intergenic (+61/‑364)	USA300HOU_RS05380 → / → rlmI	hypothetical protein/Ribosomal RNA large subunit methyltransferase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,089,908	0	T	G	5.2%	1116.7 / 36.3	321	intergenic (+61/‑364)	USA300HOU_RS05380/rlmI	hypothetical protein/Ribosomal RNA large subunit methyltransferase I
Reads supporting (aligned to +/- strand): ref base T (160/144); new base G (12/5); total (172/149)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 342 total aligned reads displayed.

ATTCGTTACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGATGTGAAT  >  USA300TCH1516_ALE/1089819‑1089998

aTTCGTTACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGata                                                                                           <  2:2951711/91‑1 (MQ=255)

 ttCGTTACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGatat                                                                                          >  1:258866/1‑91 (MQ=255)

    gTTACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATAtc                                                                                       <  2:1128954/91‑1 (MQ=255)

     ttACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATAtct                                                                                      <  1:2411708/91‑1 (MQ=255)

     ttACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATAtct                                                                                      <  1:2984917/91‑1 (MQ=255)

     ttACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATAtct                                                                                      <  1:1762235/91‑1 (MQ=255)

       aCAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGAGATATatct                                                                                    <  1:1717377/91‑1 (MQ=255)

        cAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTa                                                                                   <  1:1994114/91‑1 (MQ=255)

         aaTTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTaa                                                                                  <  2:2983545/91‑1 (MQ=255)

          aTTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAac                                                                                 >  2:2174331/1‑91 (MQ=255)

          aTTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAac                                                                                 <  2:1838653/91‑1 (MQ=255)

          aTTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAac                                                                                 <  2:3213501/91‑1 (MQ=255)

            tACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAactt                                                                               <  2:3285957/91‑1 (MQ=255)

              cGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAacttac                                                                             >  2:2654290/1‑91 (MQ=255)

              cGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAacttac                                                                             <  1:3774629/91‑1 (MQ=255)

              cGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAacttac                                                                             >  2:77491/1‑91 (MQ=255)

               gagaGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAacttact                                                                            <  2:47186/91‑1 (MQ=255)

                 gagCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAActtactta                                                                          <  1:2972661/91‑1 (MQ=255)

                  agCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAg                                                                         >  1:1025468/1‑91 (MQ=255)

                  agCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAg                                                                         >  1:926955/1‑91 (MQ=255)

                  agCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAg                                                                         >  1:1341134/1‑91 (MQ=255)

                   gCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGa                                                                        <  2:3383438/91‑1 (MQ=255)

                   gCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGa                                                                        <  2:1408769/91‑1 (MQ=255)

                      aCCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATAc                                                                     >  1:3451561/1‑91 (MQ=255)

                        cATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTGGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAg                                                                   >  2:2499540/1‑91 (MQ=255)

                        cATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATCTTGATGTGTACGATATAAGTTAGATATATATCTACCTTACTCAGATagag                                                                   >  2:2745270/1‑91 (MQ=255)

                          tAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGt                                                                 >  1:1106539/1‑91 (MQ=255)

                              tAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAAt                                                             <  2:2014072/91‑1 (MQ=255)

                              tAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAAt                                                             >  2:3814149/1‑91 (MQ=255)

                              tAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAAt                                                             >  1:218594/1‑91 (MQ=255)

                               aaaaGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATg                                                            <  1:2136343/91‑1 (MQ=255)

                                aaaGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTTAATGa                                                           >  2:2772667/1‑91 (MQ=255)

                                 aaGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGaa                                                          >  2:1055404/1‑91 (MQ=255)

                                  aGCATTATAAAAGTACTATCTATTAAACATTTTGATGGGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGGCAATGAAg                                                         <  2:437474/91‑1 (MQ=255)

                                   gCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGt                                                        <  1:2296004/91‑1 (MQ=255)

                                    cATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGtt                                                       <  2:3905670/91‑1 (MQ=255)

                                    cATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGtt                                                       <  1:3629521/91‑1 (MQ=255)

                                    cATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGtt                                                       >  1:1535185/1‑91 (MQ=255)

                                    cATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGtt                                                       >  1:1885039/1‑91 (MQ=255)

                                     aTTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGttt                                                      <  2:3599223/91‑1 (MQ=255)

                                      ttATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTa                                                     >  1:2050747/1‑91 (MQ=255)

                                      ttATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTa                                                     >  1:2637356/1‑91 (MQ=255)

                                       tataAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTAt                                                    >  1:118522/1‑91 (MQ=255)

                                        ataAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATg                                                   >  1:3485895/1‑91 (MQ=255)

                                        ataAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATg                                                   <  1:1626166/91‑1 (MQ=255)

                                        ataAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATg                                                   >  1:2882058/1‑91 (MQ=255)

                                        ataAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATg                                                   >  1:1252887/1‑91 (MQ=255)

                                        ataAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATg                                                   >  1:2808684/1‑91 (MQ=255)

                                        ataAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATg                                                   >  2:1072467/1‑91 (MQ=255)

                                         taAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATgg                                                  >  2:3170536/1‑91 (MQ=255)

                                             aGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAg                                              >  2:1908320/1‑91 (MQ=255)

                                              gTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGt                                             <  2:2405749/91‑1 (MQ=255)

                                              gTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGt                                             <  1:494328/91‑1 (MQ=255)

                                              gTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGt                                             <  1:2528081/91‑1 (MQ=255)

                                              gTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGt                                             <  1:2756800/91‑1 (MQ=255)

                                               tACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATCTCTCTAACTTACTTAGATACAGGTTAATGAAGTTTATGGATAGtt                                            >  1:1514781/1‑90 (MQ=255)

                                               tACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTa                                            >  2:535364/1‑91 (MQ=255)

                                                aCTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTAc                                           <  1:1426466/91‑1 (MQ=255)

                                                 ctatctatTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACt                                          >  2:3716425/1‑91 (MQ=255)

                                                 ctatctatTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACt                                          >  2:2612883/1‑91 (MQ=255)

                                                  tatctatTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAGGTTTATGGATAGTACtt                                         <  2:1851394/91‑1 (MQ=255)

                                                   atctatTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACttt                                        >  1:742275/1‑91 (MQ=255)

                                                     ctatTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACttttt                                      <  2:3130189/91‑1 (MQ=255)

                                                        ttAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTg                                   <  2:2309254/91‑1 (MQ=255)

                                                          aaaCATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCACTGAAGTTTATGGATAGTACTTTTTTTGTa                                 >  1:3171945/1‑91 (MQ=255)

                                                          aaaCATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTa                                 >  1:2564483/1‑91 (MQ=255)

                                                          aaaCATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTa                                 <  1:3252940/91‑1 (MQ=255)

                                                            aCATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAAc                               <  2:571525/91‑1 (MQ=255)

                                                            aCATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAAc                               >  1:2300055/1‑91 (MQ=255)

                                                            aCATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAAc                               >  1:3862671/1‑91 (MQ=255)

                                                            aCATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAAc                               <  2:82634/91‑1 (MQ=255)

                                                              aTTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTa                             >  1:3750923/1‑91 (MQ=255)

                                                                tttGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTTACTAGa                           >  2:1193553/1‑91 (MQ=255)

                                                                tttGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAgt                           >  2:906232/1‑90 (MQ=255)

                                                                tttGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGa                           <  1:372378/91‑1 (MQ=255)

                                                                tttGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGa                           <  1:3353013/91‑1 (MQ=255)

                                                                     tgtgTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttga                      >  2:931441/1‑91 (MQ=255)

                                                                     tgtgTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttga                      <  1:653891/91‑1 (MQ=255)

                                                                     tgtgTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttga                      <  1:39435/91‑1 (MQ=255)

                                                                      gtgtACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttgat                     >  2:1784040/1‑91 (MQ=255)

                                                                       tgtACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttgatt                    >  1:301143/1‑91 (MQ=255)

                                                                       tgtACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttgatt                    <  2:1663435/91‑1 (MQ=255)

                                                                       tgtACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttgatt                    >  2:1567588/1‑91 (MQ=255)

                                                                         tACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttgattga                  >  1:3105984/1‑91 (MQ=255)

                                                                         tACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttgattga                  >  2:3875236/1‑91 (MQ=255)

                                                                            gCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGGAACTAGATTTGATTGATTc               <  2:3460037/91‑1 (MQ=255)

                                                                                 aaGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTg          >  1:247338/1‑91 (MQ=255)

                                                                                 aaGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTg          <  2:1494616/91‑1 (MQ=255)

                                                                                 aaGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTg          <  2:2862315/91‑1 (MQ=255)

                                                                                  aGTTAGATATATCTCTAACTTACCTAGATACAGGTCAATGACGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGATgc         >  2:3576494/1‑90 (MQ=255)

                                                                                    ttAGATATATCTCTAACCTACTTAGATACAGGTCAATGCAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtg       >  2:659225/1‑91 (MQ=255)

                                                                                    ttAGAGATATATCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtg       >  1:484763/1‑91 (MQ=255)

                                                                                     tAGAGATATATCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGACTCAGGTGAtgt      >  1:991668/1‑91 (MQ=255)

                                                                                     tAGAGATATATCCAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGCtgt      >  1:3058509/1‑91 (MQ=255)

                                                                                      agagATATATCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtgtg     >  1:1034260/1‑91 (MQ=255)

                                                                                      aGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtgtg     >  1:906208/1‑91 (MQ=255)

                                                                                      aGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtgtg     <  1:3828847/91‑1 (MQ=255)

                                                                                       gATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTTTTTTTTGTTTTTATTGATTCAGGTGATGTGa    >  2:2273136/1‑91 (MQ=255)

                                                                                       gATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGATGTGa    >  1:1373095/1‑91 (MQ=255)

                                                                                        atatatCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGATGTGaa   <  2:1698085/91‑1 (MQ=255)

ATTCGTTACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGATGTGAAT  >  USA300TCH1516_ALE/1089819‑1089998

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