BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	604,132	Δ1 bp	6.3%	coding (1563/1638 nt)	araB →	Ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	604,132	0	G	.	6.3%	784.5 / 41.2	223	coding (1563/1638 nt)	araB	Ribulokinase
Reads supporting (aligned to +/- strand): ref base G (117/92); new base . (9/5); total (126/97)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 242 total aligned reads displayed.

AGTTGATGCTATGAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATAAAAAATTTAATTTGAC  >  USA300TCH1516_ALE/604042‑604222

 ttcgATGCTATGAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAgg                                                                                           <  2:2565266/88‑1 (MQ=255)

 gTTGATGCTATGAAAGAGCCAATTTTATACGTAATTAATCCAGATGCGGAAAACGTAGAAAGGTATGGAACATTATTTAAAG‑CTTATAAgg                                                                                           <  1:3096143/91‑1 (MQ=255)

 gTTGATGCTATGAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAgg                                                                                           <  1:1336631/91‑1 (MQ=255)

  ttGATGCTATGAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGc                                                                                          >  1:650108/1‑91 (MQ=255)

   tGATGCTATGAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCt                                                                                         >  2:127202/1‑91 (MQ=255)

    gATGCTATGAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAGCCTTATAAgctt                                                                                         <  1:603681/91‑4 (MQ=255)

     aTGCTATGAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAGCCTTATAAgcttt                                                                                        <  1:887030/91‑5 (MQ=255)

     aTGCTATGAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCttt                                                                                       <  2:349797/91‑1 (MQ=255)

      tGCTATGAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAGCCTTATAA‑GCTTTa                                                                                      <  1:1627207/91‑1 (MQ=255)

      tGCTATGAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTa                                                                                      <  1:3328169/91‑1 (MQ=255)

       gCTATGAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAGCCTTATAA‑GCTTTaa                                                                                     <  2:76483/91‑2 (MQ=255)

         tATGAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAACAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACAt                                                                                   <  1:3086852/91‑1 (MQ=255)

         aaTGAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAGCCTTATAA‑GCTTTAAAt                                                                                   <  2:489508/90‑1 (MQ=255)

          aTGAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATg                                                                                  >  2:1571931/1‑91 (MQ=255)

          aTGAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATg                                                                                  >  1:2676427/1‑91 (MQ=255)

            gAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACGTTATTTAAAG‑CTTATAAGGCTTTACATGat                                                                                <  2:1525608/91‑1 (MQ=255)

              aaGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGatat                                                                              >  1:143531/1‑91 (MQ=255)

              aaGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGatat                                                                              <  1:3755668/91‑1 (MQ=255)

               agagCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATc                                                                             <  1:1778398/91‑1 (MQ=255)

               agagCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATc                                                                             <  1:2940341/91‑1 (MQ=255)

                gagCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATcc                                                                            <  2:2970991/91‑1 (MQ=255)

                  gCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCAt                                                                          >  2:3510434/1‑91 (MQ=255)

                   ccAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATg                                                                         >  2:1743372/1‑91 (MQ=255)

                    cAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATgg                                                                        >  2:584992/1‑91 (MQ=255)

                     aaTTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACGTTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGt                                                                       >  1:1726974/1‑91 (MQ=255)

                      aTTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGtt                                                                      >  2:2114031/1‑91 (MQ=255)

                       ttttATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTta                                                                     >  2:1639794/1‑91 (MQ=255)

                       ttttATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTta                                                                     >  1:2792990/1‑91 (MQ=255)

                       ttttATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTta                                                                     >  2:3591305/1‑91 (MQ=255)

                         ttATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTtata                                                                   <  2:255340/91‑1 (MQ=255)

                           ataCGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATaaa                                                                 >  1:2732146/1‑91 (MQ=255)

                           ataCGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATaaa                                                                 <  2:650108/91‑1 (MQ=255)

                            taCGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATaaaa                                                                >  1:190695/1‑91 (MQ=255)

                                aaaTTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGCTATCCATGGTTATAAAAAAGc                                                            >  1:1819173/1‑91 (MQ=255)

                                aaaTTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGc                                                            <  2:2317263/91‑1 (MQ=255)

                                  aTTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTa                                                          >  1:2411900/1‑91 (MQ=255)

                                   ttAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTaa                                                         >  1:3337336/1‑91 (MQ=255)

                                   ttAATACAAAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTaa                                                         <  2:2863241/91‑1 (MQ=255)

                                    tAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAat                                                        <  2:1127385/91‑1 (MQ=255)

                                    tAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAat                                                        <  1:1694051/91‑1 (MQ=255)

                                     aaTCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAata                                                       <  1:2918697/91‑1 (MQ=255)

                                     aaTCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAata                                                       >  1:1603079/1‑91 (MQ=255)

                                       tCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAAtata                                                     >  1:3780471/1‑91 (MQ=255)

                                        ccAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATaa                                                    >  2:158599/1‑91 (MQ=255)

                                        ccAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATaa                                                    >  1:1135924/1‑91 (MQ=255)

                                        ccAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATaa                                                    >  1:3187962/1‑91 (MQ=255)

                                        ccAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATaa                                                    >  2:859028/1‑91 (MQ=255)

                                         cAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAAt                                                   >  1:139309/1‑91 (MQ=255)

                                         cAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAAt                                                   >  1:1499985/1‑91 (MQ=255)

                                         cAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAAt                                                   >  1:1500488/1‑91 (MQ=255)

                                          aGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATg                                                  >  1:2983476/1‑91 (MQ=255)

                                              gCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAg                                              >  1:2259833/1‑91 (MQ=255)

                                              gCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAg                                              >  1:3661458/1‑91 (MQ=255)

                                               cGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGa                                             >  2:3146639/1‑91 (MQ=255)

                                               cGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGa                                             >  2:713544/1‑91 (MQ=255)

                                                   aaaaGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATc                                         <  1:2090053/91‑1 (MQ=255)

                                                   aaaaGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATc                                         >  1:620797/1‑91 (MQ=255)

                                                    aaaGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATcc                                        <  1:169499/91‑1 (MQ=255)

                                                      aGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCag                                      <  1:962089/91‑1 (MQ=255)

                                                      aGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCag                                      <  2:688091/91‑1 (MQ=255)

                                                          cAAAGGTATGAAACTTTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGtt                                  >  2:3241703/1‑91 (MQ=255)

                                                          cAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGtt                                  >  2:222630/1‑91 (MQ=255)

                                                          cAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGtt                                  >  2:196319/1‑91 (MQ=255)

                                                          cAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGtt                                  >  2:3758013/1‑91 (MQ=255)

                                                          cAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGtt                                  >  2:1956516/1‑91 (MQ=255)

                                                           aaaGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGttt                                 <  1:984031/91‑1 (MQ=255)

                                                            aaGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTa                                <  2:78839/91‑1 (MQ=255)

                                                            aaGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTa                                <  1:1329493/91‑1 (MQ=255)

                                                            aaGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTa                                <  1:2494839/91‑1 (MQ=255)

                                                            aaGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTa                                <  1:1328919/91‑1 (MQ=255)

                                                            aaGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTa                                <  2:3467597/91‑1 (MQ=255)

                                                             aGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTaa                               <  1:955173/91‑1 (MQ=255)

                                                             aGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTaa                               >  2:1559805/1‑91 (MQ=255)

                                                               gTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAaga                             <  2:1222134/91‑1 (MQ=255)

                                                                   gAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTCTAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTg                         >  2:3481287/1‑91 (MQ=255)

                                                                    aaaCATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGa                        >  1:512562/1‑91 (MQ=255)

                                                                     aaCATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAg                       >  1:978220/1‑91 (MQ=255)

                                                                         ttattTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGAt                   <  2:1345957/91‑1 (MQ=255)

                                                                           atgtAAAGCCTTATAA‑GCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATaa                 >  1:3417329/4‑91 (MQ=255)

                                                                            tttAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATaaa                >  2:2715534/1‑91 (MQ=255)

                                                                            tttAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATaaa                >  2:1889168/1‑91 (MQ=255)

                                                                             ttAAAG‑CTTATAGGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATaaaa               >  1:2988211/1‑91 (MQ=255)

                                                                             ttAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATaaaa               <  1:228260/91‑1 (MQ=255)

                                                                             ttAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATaaaa               >  2:2654284/1‑91 (MQ=255)

                                                                             ttAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATaaaa               <  1:3510434/91‑1 (MQ=255)

                                                                              tAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATaaaaa              >  2:1060750/1‑91 (MQ=255)

                                                                              tAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATaaaaa              <  1:3441759/91‑1 (MQ=255)

                                                                              tAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATaaaaa              <  1:3461492/91‑1 (MQ=255)

                                                                                 aG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATAAAAaatt           >  2:2710300/1‑91 (MQ=255)

                                                                                 aG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATAAAAaatt           >  2:3302815/1‑91 (MQ=255)

                                                                                  g‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATAAAAaattt          <  2:883617/91‑1 (MQ=255)

                                                                                   ccTTATAA‑GCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATAAAAaattta         >  2:3625886/2‑91 (MQ=255)

                                                                                    cTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATAAAAaattta         >  2:3083269/1‑91 (MQ=255)

                                                                                     ttATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATAAAAaatttaa        <  1:2093948/91‑1 (MQ=255)

                                                                                     ttATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATAAAAaatttaa        <  1:1351681/91‑1 (MQ=255)

                                                                                       ataAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATAAAAaatttaatt      <  1:876481/91‑1 (MQ=255)

                                                                                          aGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGCTAAAAAATTTAATTCGa   >  2:3341261/1‑91 (MQ=255)

                                                                                           ggCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATAAAAAATTTAATTTGAc  <  1:67580/91‑1 (MQ=255)

                                                                                           ggCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATAAAAAATTTAATTTGAc  <  1:1003508/91‑1 (MQ=255)

                                                                                           ggCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATAAAAAATTTAATTTGAc  <  2:1159159/91‑1 (MQ=255)

AGTTGATGCTATGAAAGAGCCAATTTTATACGAAATTAATCCAGAAGCGGAAAAAGTACAAAGGTATGAAACATTATTTAAAG‑CTTATAAGGCTTTACATGATATCCATGGTTATAAAAAAGCTAATATAATGAAAGATATCCAGAGTTTAAGAGTTGAGGGATAAAAAATTTAATTTGAC  >  USA300TCH1516_ALE/604042‑604222

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