BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,654,131	C→G	8.1%	L132V (CTG→GTG)	USA300HOU_RS13650 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,654,131	0	C	G	8.1%	627.0 / 35.3	184	L132V (CTG→GTG)	USA300HOU_RS13650	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (97/72); new base G (9/6); total (106/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 192 total aligned reads displayed.

ACTGCAGGTGGTGCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTCTACTATT  >  USA300TCH1516_ALE/2654041‑2654221

 cTGCAGGTGGTGCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACt                                                                                           >  1:1751212/1‑91 (MQ=255)

    cAGGTGGTGCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGtt                                                                                        <  1:3862286/91‑1 (MQ=255)

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

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

       gtgatgCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCt                                                                                     <  1:3123569/91‑1 (MQ=255)

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

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

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

         ggtgCGTCAACGATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGGACTCATATTTGTCGCCTTTATTTTAGCAACACTGTTTCTaa                                                                                   <  2:2956398/91‑1 (MQ=255)

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

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

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

         ggtgCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTaa                                                                                   <  1:2041207/91‑1 (MQ=255)

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

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

           tgCGTCAACCATTCTTCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTaaaa                                                                                 <  1:3500855/91‑1 (MQ=255)

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

            gCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAAt                                                                                <  2:3671349/91‑1 (MQ=255)

              gTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTc                                                                              >  2:2027972/1‑91 (MQ=255)

               tCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCGGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCg                                                                             >  1:3269343/1‑91 (MQ=255)

                tcaacATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGa                                                                            <  2:3794423/87‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                        catcaAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTaa                                                                    >  2:3081178/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)

                                tttCAACTTACCGTATTGGCTGAGCGCACTCACATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCtgtg                                                            <  2:3600002/91‑1 (MQ=255)

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

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

                                   cAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCtt                                                         <  2:3876906/91‑1 (MQ=255)

                                    aaCTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTg                                                        >  1:310200/1‑91 (MQ=255)

                                     aCTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTgt                                                       >  2:3847564/1‑90 (MQ=255)

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

                                      cTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGc                                                      >  2:106587/1‑91 (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)

                                             tATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTa                                               <  2:1589623/91‑1 (MQ=255)

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

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

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

                                                      aGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCAtttt                                      <  1:798432/91‑1 (MQ=255)

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

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

                                                       gcgcACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCAttttt                                     <  1:631011/91‑1 (MQ=255)

                                                        cgcACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTa                                    >  2:1574770/1‑91 (MQ=255)

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

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

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

                                                           aCTCATATAAGTCGCCTTTATTTTAGCTACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAAtt                                 <  2:2052067/91‑1 (MQ=255)

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

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

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

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

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

                                                               atatTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGa                             <  1:2053435/91‑1 (MQ=255)

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

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

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

                                                                tatTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGAt                            <  1:3373304/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                             tATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGAtt               >  1:3212773/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                 ttAGAAACAGTGTTGCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCgg           >  1:2639644/1‑91 (MQ=255)

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

                                                                                    gCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTct        >  1:3150493/1‑91 (MQ=255)

                                                                                    gCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTct        >  1:2434588/1‑91 (MQ=255)

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

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

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

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

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

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

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

ACTGCAGGTGGTGCGTCAACCATTCATCAAAGTTTCAACTTACCGTATTGGCTGAGCGCACTCATATTAGTCGCCTTTATTTTAGCAACACTGTTTCTAAAATTCGATCGTTTAATTGCTGTGCTTGGCGGTGTTACCCCATTTTTAATTGCGATTGTCATTATGATTGCGGTCTACTATT  >  USA300TCH1516_ALE/2654041‑2654221

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