BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,089,913	C→A	5.7%	intergenic (+66/‑359)	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,913	0	C	A	5.7%	737.3 / 26.3	214	intergenic (+66/‑359)	USA300HOU_RS05380/rlmI	hypothetical protein/Ribosomal RNA large subunit methyltransferase I
Reads supporting (aligned to +/- strand): ref base C (109/91); new base A (8/4); total (117/97)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 237 total aligned reads displayed.

GTTACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGATGTGAATTAAG  >  USA300TCH1516_ALE/1089823‑1090002

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

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

 ttACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATAtct                                                                                          >  1:2417130/1‑91 (MQ=255)

   aCAATTACGAGAGCGATCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATctct                                                                                        <  2:2971368/91‑1 (MQ=255)

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

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

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

         aCGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAactta                                                                                  >  2:3159493/1‑91 (MQ=255)

         aCGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAactta                                                                                  >  2:2346194/1‑91 (MQ=255)

         aCGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAactta                                                                                  <  2:843593/91‑1 (MQ=255)

         aCGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCAATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAactta                                                                                  >  1:3358499/1‑91 (MQ=255)

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

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

           gagaGCGACCATAATTAAAAGCATTATAAAAGTGCTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAacttact                                                                                <  1:3444772/91‑1 (MQ=255)

           gagaGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAacttact                                                                                >  1:1308903/1‑91 (MQ=255)

            agagCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGAGATATATCTAacttatag                                                                               <  1:759792/91‑4 (MQ=38)

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

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

              agCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAg                                                                             <  2:2565613/91‑1 (MQ=255)

               gCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGa                                                                            >  2:2106119/1‑91 (MQ=255)

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

                 gACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATa                                                                          <  2:2426856/91‑1 (MQ=255)

                 gACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATa                                                                          <  2:910148/91‑1 (MQ=255)

                 gACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGAAATATCTCTAACTTACTTAGATa                                                                          >  2:3123364/1‑91 (MQ=255)

                   ccATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACa                                                                        >  1:99876/1‑91 (MQ=255)

                   ccATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACa                                                                        >  2:2627264/1‑91 (MQ=255)

                   ccATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACa                                                                        >  1:25157/1‑91 (MQ=255)

                     aTAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAgg                                                                      <  2:1582302/91‑1 (MQ=255)

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

                      tAATTAAAAGCATTAGAAAAGTACTGGCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGt                                                                     <  2:3361926/91‑1 (MQ=255)

                          tAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATTTCTAACTTACTTAGATACAGGTCAAt                                                                 <  2:2212609/91‑1 (MQ=255)

                           aaaaGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATg                                                                >  2:75676/1‑91 (MQ=255)

                           aaaaGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATg                                                                >  1:410012/1‑91 (MQ=255)

                            aaaGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGa                                                               >  2:1143847/1‑91 (MQ=255)

                            aaaGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGa                                                               >  1:3808653/1‑91 (MQ=255)

                             aaGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGaa                                                              <  2:621884/91‑1 (MQ=255)

                             aaGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGaa                                                              <  1:574073/91‑1 (MQ=255)

                              aGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAg                                                             >  2:195022/1‑91 (MQ=255)

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

                                cATTATAAAAGTACTATCTATTAAAGATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTGCTTAGATACAGGTCAATGAAGtt                                                           <  1:3597427/91‑1 (MQ=255)

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

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

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

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

                                      aaaaGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGa                                                     <  2:1631877/91‑1 (MQ=255)

                                      aaaaGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGa                                                     <  2:1994006/91‑1 (MQ=255)

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

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

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

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

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

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

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

                                                   atTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACttttttt                                        >  1:1624810/1‑91 (MQ=255)

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

                                                         cATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACt                                  >  2:2204239/1‑91 (MQ=255)

                                                         cATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACt                                  >  1:1184410/1‑91 (MQ=255)

                                                         cATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACt                                  <  1:2280322/91‑1 (MQ=255)

                                                         cATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACt                                  <  2:2479626/91‑1 (MQ=255)

                                                         cATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACt                                  <  1:3710276/91‑1 (MQ=255)

                                                              tGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGAtt                             >  1:2937883/1‑91 (MQ=255)

                                                               gATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGAttt                            <  2:1594554/91‑1 (MQ=255)

                                                                aTGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttg                           >  2:1814579/1‑91 (MQ=255)

                                                                aTGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttg                           >  1:229709/1‑91 (MQ=255)

                                                                aTGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttg                           >  1:906079/1‑91 (MQ=255)

                                                                aTGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttg                           <  2:2703334/91‑1 (MQ=255)

                                                                aTGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttg                           >  1:1990121/1‑91 (MQ=255)

                                                                aTGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttg                           >  2:418512/1‑91 (MQ=255)

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

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

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

                                                                  gtgtACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttgat                         <  2:2243047/91‑1 (MQ=255)

                                                                  gtgtACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttgat                         <  1:2723053/91‑1 (MQ=255)

                                                                    gtACGCTATAAGTTAGATATATCTCTAACTTACTTAGCTACAGGTCAATGACGTTTATGGATAGTGCTTTTTTTGTAACTAGATtcgattg                       >  1:3612559/1‑91 (MQ=255)

                                                                        gCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTc                   <  1:826622/91‑1 (MQ=255)

                                                                         cTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCa                  <  2:2810818/91‑1 (MQ=255)

                                                                         cTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCa                  >  2:351657/1‑91 (MQ=255)

                                                                          tataAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAg                 <  2:3457864/91‑1 (MQ=255)

                                                                          agtaAGTTAGAGATATATCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAg                 >  1:1282028/3‑91 (MQ=39)

                                                                          agtaAGTTAGAGATATATCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAg                 >  1:1043447/3‑91 (MQ=39)

                                                                           gtaAGTTAGAGATATATCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAgg                >  1:1135415/2‑91 (MQ=255)

                                                                           ataAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAgg                >  1:1861337/1‑91 (MQ=255)

                                                                           ataAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAgg                <  2:1531274/91‑1 (MQ=255)

                                                                           ataAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAgg                <  1:61909/91‑1 (MQ=255)

                                                                           ataAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAgg                <  1:61200/91‑1 (MQ=255)

                                                                            taAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGt               <  1:519386/91‑1 (MQ=255)

                                                                                ttAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtg           >  2:1993054/1‑91 (MQ=255)

                                                                                ttAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtg           >  2:387643/1‑91 (MQ=255)

                                                                                ttAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtg           >  1:3855853/1‑91 (MQ=255)

                                                                                 tAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtgt          <  1:967870/91‑1 (MQ=255)

                                                                                 tAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtgt          <  1:2148096/91‑1 (MQ=255)

                                                                                 tAGAGATATATCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtgt          >  1:1508563/1‑91 (MQ=255)

                                                                                  agagatatatCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtgtg         >  1:3335462/5‑91 (MQ=255)

                                                                                  agagatatatCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtgtg         >  1:2870307/5‑91 (MQ=255)

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

                                                                                      atatCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGATGTGAAtt     >  2:706287/1‑91 (MQ=255)

                                                                                      atatCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGATGTGAAtt     <  1:1843798/91‑1 (MQ=255)

                                                                                       tatCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGATGTGAATTa    >  2:2327272/1‑91 (MQ=255)

                                                                                        atCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTTTTTCTTTTTTTGATTGATTCAGGTGATGTGAATTaa   >  2:3757054/1‑91 (MQ=255)

                                                                                         tctcTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGATGTGAATTAAg  <  1:212148/91‑1 (MQ=255)

GTTACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGATGTGAATTAAG  >  USA300TCH1516_ALE/1089823‑1090002

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