BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,674,892	C→G	6.3%	D42H (GAT→CAT)	dnaG ←	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,674,892	0	C	G	6.3%	774.3 / 32.2	221	D42H (GAT→CAT)	dnaG	DNA primase
Reads supporting (aligned to +/- strand): ref base C (118/89); new base G (8/6); total (126/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

Only 100 of 253 total aligned reads displayed.

TAAATTGAAAAACATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCGGTTTTATCTTT  >  USA300TCH1516_ALE/1674802‑1674982

 aaaTTGAAAAACATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTtcatca                                                                                           <  1:3205359/91‑1 (MQ=255)

   aTTGAAAAACATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCatgatg                                                                                         <  1:462231/91‑1 (MQ=255)

   aTTGAAAAACATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATg                                                                                         >  1:1343532/1‑91 (MQ=255)

     tGAAAAACATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGTCCTTTTCAtgatga                                                                                        <  1:2604749/91‑1 (MQ=255)

     tGAAAAACATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGaa                                                                                       <  2:3791829/91‑1 (MQ=255)

      gAAAATCATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGaaa                                                                                      <  2:3574717/91‑1 (MQ=255)

       aaaaaCATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGaaaa                                                                                     >  1:2044479/1‑91 (MQ=255)

        aaaaCATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAg                                                                                    <  1:3372857/91‑1 (MQ=255)

        aaaaCATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAg                                                                                    <  1:2830935/91‑1 (MQ=255)

        aaaaCATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAg                                                                                    >  2:3125322/1‑91 (MQ=255)

         aaaCATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAgg                                                                                   >  2:2878111/1‑91 (MQ=255)

         aaaCATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAgg                                                                                   <  1:3586141/91‑1 (MQ=255)

         aaaCATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAgg                                                                                   >  1:553732/1‑91 (MQ=255)

         aaaCATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAgg                                                                                   <  2:2316238/91‑1 (MQ=255)

          aaCATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGTCCTTTTCATGATGAAAAga                                                                                   <  1:3662497/91‑2 (MQ=255)

          aaCATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGa                                                                                  <  1:2760980/91‑1 (MQ=255)

          aaCATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGa                                                                                  <  1:2619191/91‑1 (MQ=255)

           aCATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGac                                                                                 <  1:1156177/91‑1 (MQ=255)

             aTTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGTCCTTTTCATGATGAAAAgacac                                                                                <  2:820001/91‑5 (MQ=38)

             aTTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGacac                                                                               <  2:322033/91‑1 (MQ=255)

               tGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACaa                                                                             >  1:90449/1‑91 (MQ=255)

                gCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACaaa                                                                            >  1:3284605/1‑91 (MQ=255)

                gCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACaaa                                                                            >  2:1485544/1‑91 (MQ=255)

                  cACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAAcc                                                                          >  1:1800710/1‑91 (MQ=255)

                  cACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAAcc                                                                          >  1:3625133/1‑91 (MQ=255)

                   aCCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCt                                                                         >  1:1531600/1‑91 (MQ=255)

                    ccTTTTTTACACCCAAAACAATGGCAAATTTGTTTATCTTCAGAAACAGTGAATGAAGGGGT‑CTTTTCATCATGAAAAGGACACAAAccga                                                                        <  1:2793913/91‑3 (MQ=255)

                    ccTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCta                                                                        >  1:3784903/1‑91 (MQ=255)

                     cTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCtat                                                                       >  1:243898/1‑91 (MQ=255)

                       tttttACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCtatat                                                                     <  2:1465228/91‑1 (MQ=255)

                          ttACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAAt                                                                  >  1:1481533/1‑91 (MQ=255)

                           tACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAAtt                                                                 >  2:3724959/1‑91 (MQ=255)

                            aCAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTg                                                                >  1:906799/1‑91 (MQ=255)

                              aaCCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCg                                                              >  1:1667231/1‑91 (MQ=255)

                               aCCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGTTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGt                                                             >  1:3882516/1‑91 (MQ=255)

                                 cAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTcc                                                           <  2:424408/91‑1 (MQ=255)

                                  aaaaCAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCct                                                          >  2:1065275/1‑91 (MQ=255)

                                     aCAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCctctc                                                       <  1:1013667/91‑1 (MQ=255)

                                     aCAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCctctc                                                       <  1:1857596/91‑1 (MQ=255)

                                     aCAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCatctc                                                       >  2:1230872/1‑91 (MQ=255)

                                       aaTGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCtt                                                     <  1:160560/91‑1 (MQ=255)

                                       aaTGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCtt                                                     <  2:1003821/91‑1 (MQ=255)

                                        aTGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCttt                                                    <  1:2816384/91‑1 (MQ=255)

                                        aTGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCttt                                                    >  2:67379/1‑91 (MQ=255)

                                         tGACAAATTTGTTTATCTTCAGAAACTGTTAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCtttt                                                   >  2:3383763/1‑91 (MQ=255)

                                         tGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTAATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCtttt                                                   >  2:3483714/1‑91 (MQ=255)

                                           aCAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCt                                                 >  2:1262241/1‑91 (MQ=255)

                                               aTTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAAtt                                             >  2:2335314/1‑91 (MQ=255)

                                                 ttGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAAtttt                                           <  2:2785442/91‑1 (MQ=255)

                                                 ttGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAAtttt                                           <  2:1922181/91‑1 (MQ=255)

                                                  tGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTa                                          <  1:3803468/91‑1 (MQ=255)

                                                  tGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTa                                          <  2:804445/91‑1 (MQ=255)

                                                   gTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTAc                                         >  1:1032229/1‑91 (MQ=255)

                                                   gTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTAc                                         >  1:2863801/1‑91 (MQ=255)

                                                      tATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACata                                      >  1:2482351/1‑91 (MQ=255)

                                                      tATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACata                                      <  1:1818402/91‑1 (MQ=255)

                                                      tATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACata                                      <  2:1413585/91‑1 (MQ=255)

                                                       aTCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACatat                                     <  1:1958317/91‑1 (MQ=255)

                                                        tCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATAtt                                    >  2:2172714/1‑91 (MQ=255)

                                                         cTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTc                                   >  1:596250/1‑91 (MQ=255)

                                                          ttCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCa                                  >  1:321490/1‑91 (MQ=255)

                                                           tCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCAc                                 <  1:165910/91‑1 (MQ=255)

                                                            cAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACt                                <  2:2919399/91‑1 (MQ=255)

                                                            cAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACt                                >  1:440121/1‑91 (MQ=255)

                                                            cAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACt                                <  2:1912882/91‑1 (MQ=255)

                                                               aaaCTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCg                                                              <  2:2618598/58‑1 (MQ=255)

                                                               aaaCTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTAc                             >  2:2869862/1‑91 (MQ=255)

                                                                aaCTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTCCATATTCACTTAcc                            >  1:1263909/1‑91 (MQ=255)

                                                                 aCTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAGCCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCa                           >  1:3369408/1‑91 (MQ=255)

                                                                  cTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCaa                          <  1:2321909/91‑1 (MQ=255)

                                                                   tGTAAATGACGGTGC‑CTTTTCATCATCAAACGGACACAAACCTATATAATCGCGTCCTCTCCCTTCTAATTTTACATATCCACTTACCAAg                         >  2:1149427/1‑91 (MQ=255)

                                                                   tGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAg                         <  2:2730420/91‑1 (MQ=255)

                                                                    gTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGt                        >  2:2013479/1‑91 (MQ=255)

                                                                     tAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTc                       <  2:1633709/91‑1 (MQ=255)

                                                                       aaTGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTa                     >  2:1894695/1‑91 (MQ=255)

                                                                       aaTGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTa                     >  2:866436/1‑91 (MQ=255)

                                                                          gAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTaaaa                  <  2:1526509/91‑1 (MQ=255)

                                                                          gAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGCCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTaaaa                  <  2:1636728/91‑1 (MQ=255)

                                                                           aaGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAAt                 <  1:2332812/91‑1 (MQ=255)

                                                                             ggTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAAGTTTACATATTCACTTACCAAGTCTAAAATGt               <  2:3350578/91‑1 (MQ=255)

                                                                               tgt‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACACATTCACTTACCAAGTCTAAAATGTCg             <  2:2211546/91‑1 (MQ=255)

                                                                                gt‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTTTTTTTATATATTCACTTACCAAGTCTAAAATGTCgg            >  2:3670167/1‑91 (MQ=255)

                                                                                gt‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCgg            <  1:1367211/91‑1 (MQ=255)

                                                                                gt‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCgg            <  1:3249878/91‑1 (MQ=255)

                                                                                 tccTTTTCATGATGAAAA‑GACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCGGt           >  2:3159169/3‑91 (MQ=255)

                                                                                 t‑CTTTTCATCATGAAAGGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCGGt           >  2:2587991/1‑91 (MQ=255)

                                                                                   cTTTTCATGATGAAAA‑GACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCGGttt         >  1:2698366/1‑91 (MQ=255)

                                                                                    ttttCATGATGAAAA‑GACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCGGtttt        >  1:2211208/1‑91 (MQ=255)

                                                                                    ttttCATGATGAAAA‑AACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACTTATTCACTTACCAAGTCTAAAATGTCGGtttt        >  1:3305055/1‑91 (MQ=255)

                                                                                     tttCATGATGAAAA‑GACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCGGTTTTa       >  1:3772653/1‑91 (MQ=255)

                                                                                     tttCATCATGATGAAGATTCTAATCTTTATGATTGCGTCCTCTCTTTTCTTATTTTACATATTCACTTACCAAGTCTAAAATGTCGGtttt        <  2:294067/91‑1 (MQ=255)

                                                                                     tttCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCGGtttt        >  2:1083450/1‑91 (MQ=255)

                                                                                        catcatGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCGGTTTTATc     >  2:1950026/1‑91 (MQ=255)

                                                                                         atcatGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCGGTTTTATCt    <  2:899560/91‑1 (MQ=255)

                                                                                         atcatGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCGGTTTTATCt    <  1:191869/91‑1 (MQ=255)

                                                                                         atcatGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCGGTTTTATCt    <  1:1073109/91‑1 (MQ=255)

                                                                                          tcatGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCGGTTTTATCtt   <  2:532270/91‑1 (MQ=255)

                                                                                          tcatGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCGGTTTTATCtt   <  2:3521296/91‑1 (MQ=255)

                                                                                           catGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCGGTTTTATCttt  >  1:769151/1‑91 (MQ=255)

                                                                                           catGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCGGTTTTATCttt  <  1:348461/91‑1 (MQ=255)

TAAATTGAAAAACATTGCCACCTTTTTTACAACCAAAACAATGACAAATTTGTTTATCTTCAGAAACTGTAAATGAAGGTGT‑CTTTTCATCATGAAAAGGACACAAACCTATATAATTGCGTCCTCTCTTTTCTAATTTTACATATTCACTTACCAAGTCTAAAATGTCGGTTTTATCTTT  >  USA300TCH1516_ALE/1674802‑1674982

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