BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,138,307	A→T	5.2%	F104I (TTT→ATT) ‡	USA300HOU_RS10905 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,138,307	0	A	T	5.2%	867.6 / 26.5	251	F104I (TTT→ATT) ‡	USA300HOU_RS10905	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (123/115); new base T (10/3); total (133/118)
Fisher's exact test for biased strand distribution p-value = 9.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 265 total aligned reads displayed.

TATTCTTTTGAAAAATATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTAAATGGTAA  >  USA300TCH1516_ALE/2138217‑2138397

      tttGAAAAATATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAAtattat                                                                                      <  2:1634317/91‑1 (MQ=255)

        tGAAAAATATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTataa                                                                                    >  1:257111/1‑91 (MQ=255)

         gAAAAATATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTataat                                                                                   <  2:3880293/91‑1 (MQ=255)

         gAAAAATATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTataat                                                                                   <  2:1873337/91‑1 (MQ=255)

          aaaaaTATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTataata                                                                                  >  1:1820039/1‑91 (MQ=255)

           aaaaTATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATATTTATTataatac                                                                                 <  1:3167819/91‑2 (MQ=255)

           aaaaTATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAg                                                                                 <  2:192010/91‑1 (MQ=255)

           aaaaTATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAg                                                                                 >  1:2294386/1‑91 (MQ=255)

           aaaaTATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAAATAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAg                                                                                 <  2:21895/91‑1 (MQ=255)

            aaaTATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATATTTATTataatact                                                                                <  1:1138650/91‑3 (MQ=39)

            aaaTATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGc                                                                                <  2:2194436/91‑1 (MQ=255)

            aaaTATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGc                                                                                <  2:1253761/91‑1 (MQ=255)

             aaTATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCt                                                                               >  1:3738787/1‑91 (MQ=255)

               tatCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTAt                                                                             >  2:1548887/1‑91 (MQ=255)

               tatCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTAt                                                                             <  1:2107782/91‑1 (MQ=255)

               tatCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTAt                                                                             <  1:1245142/91‑1 (MQ=255)

                atCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATc                                                                            <  2:348583/91‑1 (MQ=255)

                   tGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTa                                                                         >  2:3625897/1‑91 (MQ=255)

                   tGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTa                                                                         >  2:1659769/1‑91 (MQ=255)

                   tGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTa                                                                         >  1:984452/1‑91 (MQ=255)

                   tGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTa                                                                         <  1:853109/91‑1 (MQ=255)

                    gCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTAt                                                                        >  1:2869583/1‑91 (MQ=255)

                     cATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTAtt                                                                       <  2:1491429/91‑1 (MQ=255)

                     cATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTAtt                                                                       <  1:767554/91‑1 (MQ=255)

                          cTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATa                                                                  <  1:1130506/91‑1 (MQ=255)

                            tCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGt                                                                <  2:580594/91‑1 (MQ=255)

                            tCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGt                                                                >  2:234856/1‑91 (MQ=255)

                             cacaGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGtt                                                               >  2:2924704/1‑91 (MQ=255)

                             cacaGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGtt                                                               >  2:2494526/1‑91 (MQ=255)

                             cacaGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGtt                                                               >  2:2459664/1‑91 (MQ=255)

                             cacaGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGtt                                                               <  2:1787352/91‑1 (MQ=255)

                               caGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTa                                                             >  2:954824/1‑91 (MQ=255)

                               caGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTa                                                             >  2:2885013/1‑91 (MQ=255)

                                aGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAg                                                            <  1:3393645/91‑1 (MQ=255)

                                   aaTTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAat                                                         <  1:2718824/91‑1 (MQ=255)

                                     ttttATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAatat                                                       <  2:3726880/91‑1 (MQ=255)

                                     ttttATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAatat                                                       >  1:3469210/1‑91 (MQ=255)

                                     ttttATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAatat                                                       >  2:3502203/1‑91 (MQ=255)

                                     ttttATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAatat                                                       <  2:100708/91‑1 (MQ=255)

                                      tttATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTa                                                             >  1:1557003/1‑84 (MQ=255)

                                      tttATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATg                                                      <  2:1203496/91‑1 (MQ=255)

                                       ttATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGc                                                     >  2:1204962/1‑91 (MQ=255)

                                         ataACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCaa                                                   >  2:3797168/1‑91 (MQ=255)

                                          taACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCaaa                                                  <  2:3027117/91‑1 (MQ=255)

                                            aCAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACt                                                <  1:1291432/91‑1 (MQ=255)

                                             cAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACt                                                >  2:599988/1‑90 (MQ=255)

                                             cAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTg                                               <  1:2004318/91‑1 (MQ=255)

                                             cAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTg                                               <  1:289369/91‑1 (MQ=255)

                                             cAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTg                                               <  1:791245/91‑1 (MQ=255)

                                              aGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGa                                              <  1:1621925/91‑1 (MQ=255)

                                                 ttAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaa                                           >  1:2138577/1‑91 (MQ=255)

                                                   aataaAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaat                                         >  2:2024168/1‑91 (MQ=255)

                                                    ataaAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaata                                        >  2:3409299/1‑91 (MQ=255)

                                                    ataaAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCAATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaata                                        >  2:2872844/1‑91 (MQ=255)

                                                    ataaAATTAAAAAAACATAACAACAAGTATTATAATAAATACTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaata                                        >  2:2952364/1‑91 (MQ=255)

                                                      aaaaTTAGAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaataaa                                      >  2:2730025/1‑91 (MQ=255)

                                                      aaaaTTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaataaa                                      >  1:1686435/1‑91 (MQ=255)

                                                       aaaTTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaataaaa                                     <  2:1701326/91‑1 (MQ=255)

                                                       aaaTTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaataaaa                                     >  1:1734914/1‑91 (MQ=255)

                                                       aaaTTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaataaaa                                     >  1:387490/1‑91 (MQ=255)

                                                         aTTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaataaaaat                                   >  1:3557443/1‑91 (MQ=255)

                                                         aTTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaataaaaat                                   <  1:3298932/91‑1 (MQ=255)

                                                          ttAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATc                                  <  1:597999/91‑1 (MQ=255)

                                                           tAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATca                                 <  1:610374/91‑1 (MQ=255)

                                                           tAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATca                                 <  2:1795111/91‑1 (MQ=255)

                                                             aaaaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAACTcaca                               >  2:3581646/1‑91 (MQ=255)

                                                             aaaaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATcaca                               >  1:115771/1‑91 (MQ=255)

                                                             aaaaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATcaca                               >  2:873398/1‑91 (MQ=255)

                                                              aaaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACaa                              <  1:2810452/91‑1 (MQ=255)

                                                               aaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAAt                             >  2:722912/1‑91 (MQ=255)

                                                               aaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAAt                             >  1:1896634/1‑91 (MQ=255)

                                                               aaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAAt                             >  2:3778852/1‑91 (MQ=255)

                                                                  aCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTaa                          <  1:527960/91‑1 (MQ=255)

                                                                     tAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAAt                       >  2:3459276/1‑91 (MQ=255)

                                                                     tAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAAt                       >  1:404382/1‑91 (MQ=255)

                                                                      aacaacAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATa                      >  1:300938/1‑91 (MQ=255)

                                                                       acaacaAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATAc                     <  2:2971165/91‑1 (MQ=255)

                                                                       acaacaAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATAc                     >  2:3258296/1‑91 (MQ=255)

                                                                         aacaaGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTa                   <  1:234856/91‑1 (MQ=255)

                                                                            aaGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGt                >  2:2744383/1‑91 (MQ=255)

                                                                            aaGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGt                >  1:100661/1‑91 (MQ=255)

                                                                            aaGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGt                <  1:1643371/91‑1 (MQ=255)

                                                                             aGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTa               >  1:642967/1‑91 (MQ=255)

                                                                              gTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTaa              <  2:522670/91‑1 (MQ=255)

                                                                                attatAATATTTATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATg            >  1:334450/1‑91 (MQ=255)

                                                                                attatAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATg            >  2:2647639/1‑91 (MQ=255)

                                                                                attatAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATg            <  2:629789/91‑1 (MQ=255)

                                                                                 ttatAATATTTATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGt           >  2:2792576/1‑91 (MQ=255)

                                                                                 ttatAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAAt                       >  1:445524/1‑79 (MQ=255)

                                                                                 ttatAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAAt                       <  2:445627/79‑1 (MQ=255)

                                                                                  tataATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTa          >  2:1131016/1‑91 (MQ=255)

                                                                                   ataataAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTaa         >  2:459714/1‑91 (MQ=255)

                                                                                    taataaATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTaaa        <  2:684482/91‑1 (MQ=255)

                                                                                     aataaataTTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTAAAt       <  2:791215/91‑1 (MQ=255)

                                                                                      ataaataTTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTAAATg      <  1:3502153/91‑1 (MQ=255)

                                                                                        aaataTTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTAAATGGt    >  1:470619/1‑91 (MQ=255)

                                                                                        aaataTTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTAAATGGt    <  2:1363134/91‑1 (MQ=255)

                                                                                          atatTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTAAATGGTaa  <  1:1370918/91‑1 (MQ=255)

                                                                                          atatTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTAAATGGTaa  <  2:3557443/91‑1 (MQ=255)

                                                                                          atatTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTAAATGGTaa  >  1:2641664/1‑91 (MQ=255)

TATTCTTTTGAAAAATATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTAAATGGTAA  >  USA300TCH1516_ALE/2138217‑2138397

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