BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,654,136	T→G	7.4%	F133L (TTT→TTG)	USA300HOU_RS13650 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,654,136	0	T	G	7.4%	615.2 / 35.8	189	F133L (TTT→TTG)	USA300HOU_RS13650	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (96/79); new base G (8/6); total (104/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 195 total aligned reads displayed.

GGTGGTGCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTCTACTATTTCAC  >  USA300TCH1516_ALE/2654047‑2654225

 gtggtgCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTctct                                                                                         >  1:3127809/1‑91 (MQ=255)

 gtggtgCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCt                                                                                         >  1:1867701/1‑91 (MQ=255)

 gtggtgCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCt                                                                                         <  1:1012052/91‑1 (MQ=255)

  tggtgCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTa                                                                                        <  2:3327430/91‑1 (MQ=255)

  tggtgCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTa                                                                                        >  1:1587855/1‑91 (MQ=255)

  tggtgCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGAAACAGTGTTGCTa                                                                                        <  1:3184671/91‑1 (MQ=255)

  tggtgCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGAAACAGTGTTGCTa                                                                                        <  2:3880460/91‑1 (MQ=255)

  tggtgCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGAAACAGTGTTGCTa                                                                                        <  1:623582/91‑1 (MQ=255)

   ggtgCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTaa                                                                                       >  2:435233/1‑91 (MQ=255)

   ggtgCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTaa                                                                                       >  2:2147998/1‑91 (MQ=255)

   ggtgCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTaa                                                                                       >  1:1711144/1‑91 (MQ=255)

    gtgCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGAAACAGTGTTGCTaaa                                                                                      <  1:102199/91‑1 (MQ=255)

      gCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAAt                                                                                    >  1:1536054/1‑91 (MQ=255)

      gCGTCAACCATTCATCAAAGTTACAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAAt                                                                                    <  2:3132025/91‑1 (MQ=255)

           aaCCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGAt                                                                               <  2:1194509/91‑1 (MQ=255)

           aaCCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCACACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGAt                                                                               <  2:2599528/91‑1 (MQ=255)

             ccATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCg                                                                             <  2:2960630/91‑1 (MQ=255)

              cATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGt                                                                            <  2:1485428/91‑1 (MQ=255)

              cATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGt                                                                            >  2:1352788/1‑91 (MQ=255)

              cATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGt                                                                            >  2:3507169/1‑91 (MQ=255)

               aTTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGtt                                                                           >  2:3241587/1‑91 (MQ=255)

               aTTCATCAAAGTTTCAACGTACCGTATTGGGTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGtt                                                                           <  2:1029982/91‑1 (MQ=255)

                 tcatcaAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTa                                                                         >  2:753793/1‑91 (MQ=255)

                  catcaAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTaa                                                                        <  1:483278/91‑1 (MQ=255)

                  catcaAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTaa                                                                        <  1:1338432/91‑1 (MQ=255)

                     caAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTg                                                                     >  2:2077076/1‑91 (MQ=255)

                      aaaGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGc                                                                    <  2:929571/91‑1 (MQ=255)

                       aaGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCt                                                                   <  1:2402198/91‑1 (MQ=255)

                        aGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCtg                                                                  >  2:3169585/1‑91 (MQ=255)

                          tttCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCtgtg                                                                >  1:493046/1‑91 (MQ=255)

                           ttCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTGAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGc                                                               <  2:2951242/91‑1 (MQ=255)

                            tCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCATCACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCt                                                              <  1:3278578/91‑1 (MQ=255)

                            tCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCt                                                              <  2:1587855/91‑1 (MQ=255)

                            tCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCt                                                              <  2:3642314/91‑1 (MQ=255)

                            tCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCt                                                              <  2:2309300/91‑1 (MQ=255)

                               aCTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTgg                                                           <  1:1172622/91‑1 (MQ=255)

                                   aCCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGgt                                                       <  1:1620244/91‑1 (MQ=255)

                                      gTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGtt                                                    <  1:2747734/91‑1 (MQ=255)

                                       tATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGTCGGTGTTa                                                   >  2:3363625/1‑91 (MQ=255)

                                       tATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTa                                                   >  2:1421523/1‑91 (MQ=255)

                                         ttGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTAcc                                                 <  2:1805575/91‑1 (MQ=255)

                                           ggCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTAcccc                                               <  2:2899486/91‑1 (MQ=255)

                                            gCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCa                                              <  1:2476862/91‑1 (MQ=255)

                                             cTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCAt                                             >  2:1918925/1‑91 (MQ=255)

                                               gAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCAttt                                           <  1:3528536/91‑1 (MQ=255)

                                                aGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCAtttt                                          >  2:1111767/1‑91 (MQ=255)

                                                   gcACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTaa                                       <  2:38595/91‑1 (MQ=255)

                                                    cACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAAt                                      <  1:269972/91‑1 (MQ=255)

                                                    cACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAAt                                      >  2:1987825/1‑91 (MQ=255)

                                                    cACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAAt                                      <  1:800405/91‑1 (MQ=255)

                                                     aCTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAAtt                                     >  2:2408155/1‑91 (MQ=255)

                                                     aCTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAAtt                                     >  2:1988950/1‑91 (MQ=255)

                                                     aCTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAAtt                                     >  1:224373/1‑91 (MQ=255)

                                                      cTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTg                                    >  1:218869/1‑91 (MQ=255)

                                                       tCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGc                                   >  2:837680/1‑91 (MQ=255)

                                                       tCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGc                                   <  1:1131573/91‑1 (MQ=255)

                                                        cATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCg                                  >  1:1457637/1‑91 (MQ=255)

                                                          tatTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCAttt                                           >  2:2871571/1‑80 (MQ=255)

                                                          tatTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCAttt                                           <  1:2871571/80‑1 (MQ=255)

                                                          tatTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGAGGTGTTACCCCATTTTTAATTGCGAt                                >  2:1205157/1‑91 (MQ=255)

                                                           atTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGAtt                               >  1:1988559/1‑91 (MQ=255)

                                                           atTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGAtt                               >  2:764587/1‑91 (MQ=255)

                                                            ttAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTg                              >  1:2044619/1‑91 (MQ=255)

                                                              aGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGCTTGTc                            >  1:1243531/1‑91 (MQ=255)

                                                              aGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTc                            >  2:2181642/1‑91 (MQ=255)

                                                               gTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCa                           >  2:1087813/1‑91 (MQ=255)

                                                               gTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCa                           >  1:763989/1‑91 (MQ=255)

                                                                  gCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTa                        >  1:3673386/1‑91 (MQ=255)

                                                                  gCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTa                        <  1:525500/91‑1 (MQ=255)

                                                                      ttATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGAt                    <  1:753793/91‑1 (MQ=255)

                                                                       tATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGAtt                   <  2:2044619/91‑1 (MQ=255)

                                                                        aTTTTAGCAACAGTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTg                  >  1:2971570/1‑91 (MQ=255)

                                                                        aTTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTg                  >  1:3632037/1‑91 (MQ=255)

                                                                        aTTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTg                  >  2:1618946/1‑91 (MQ=255)

                                                                        aTTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTg                  >  1:2012536/1‑91 (MQ=255)

                                                                        aTTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTg                  >  2:521354/1‑91 (MQ=255)

                                                                         ttttAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGc                 >  2:3272484/1‑91 (MQ=255)

                                                                         ttttAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGc                 >  2:1802134/1‑91 (MQ=255)

                                                                         ttttAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGc                 >  2:285696/1‑91 (MQ=255)

                                                                         ttttAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGc                 >  2:3776895/1‑91 (MQ=255)

                                                                          tttAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCg                <  1:2414929/91‑1 (MQ=255)

                                                                          tttAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCg                <  1:1992740/91‑1 (MQ=255)

                                                                           ttAGAAACAGTGTTGCTAAAATTCGATCGTTTAATTGCTGTGCTTGGTGGTGTTACCCCATTTCTAATTGCGATTGTCATTATGATTGCgg               >  2:1885313/1‑91 (MQ=255)

                                                                            tagaaaCAGTGTTGCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGt              >  2:3032459/5‑91 (MQ=255)

                                                                            tAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGt              <  1:3393620/91‑1 (MQ=255)

                                                                             agaaaCAGTGTTGCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTc             >  1:3119241/4‑91 (MQ=255)

                                                                             agaaaCAGTGTTGCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTc             >  1:3256240/4‑91 (MQ=255)

                                                                             agaaaCAGTGTTGCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACACCATTTTTAATTGCGATTGTCATTATGATTGCGGTc             >  2:3837177/4‑91 (MQ=255)

                                                                              gCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTct            <  2:359202/91‑1 (MQ=255)

                                                                               cAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTcta           <  1:1421523/91‑1 (MQ=255)

                                                                                aaCACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTctac          >  1:705127/1‑91 (MQ=255)

                                                                                 acacTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTctact         <  1:1480310/91‑1 (MQ=255)

                                                                                  cacTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTACGATTGTCATTATGATTGCGGTctacta        <  1:3869024/91‑1 (MQ=255)

                                                                                   acTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTCTACTAt       >  2:327866/1‑91 (MQ=255)

                                                                                    cTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTCTACTAtt      <  1:3363625/91‑1 (MQ=255)

                                                                                      gtgtCTAAAATTCGCTCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTCTACTATTTc    >  2:2396708/4‑91 (MQ=255)

                                                                                      gTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTCTACTATTTc    <  2:1160879/91‑1 (MQ=255)

                                                                                      gTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTCTACTATTTc    <  1:304456/91‑1 (MQ=255)

                                                                                       tttCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTCTACTATTTca   >  1:3524673/1‑91 (MQ=255)

                                                                                        ttCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTCTACTATTTcac  >  2:1665788/1‑91 (MQ=255)

GGTGGTGCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTCTACTATTTCAC  >  USA300TCH1516_ALE/2654047‑2654225

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