BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,654,418	G→C	5.2%	H56Q (CAC→CAG)	USA300HOU_RS08220 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,654,418	0	G	C	5.2%	338.8 / 6.8	104	H56Q (CAC→CAG)	USA300HOU_RS08220	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (51/47); new base C (3/3); total (54/50)
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.94e-01

Only 100 of 109 total aligned reads displayed.

GATTTATAGTAGATATTTGCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTAAGTCTAATTAAGTCTGGAA  >  USA300TCH1516_ALE/1654278‑1654557

gATTTATAGTAGATATTTGCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAACATGATACgg                                                                                                                                             >  3:197195/1‑141 (MQ=255)

 aTTTATAGTAGATATTTGCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGgt                                                                                                                                            >  1:413347/1‑141 (MQ=255)

 aTTTATAGTAGATATTTGCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGgt                                                                                                                                            <  2:114650/141‑1 (MQ=255)

     ataGTAGATATTTGCAGTAAATGCTGTAACATTATTACTCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTg                                                                                                                                        >  1:784251/1‑141 (MQ=255)

      tagtagATATTTGCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGc                                                                                                                                       <  2:532649/141‑1 (MQ=255)

       tttagTTATTTTCAGTAAATGCTGGAACATTATTAATCATTGTTATATTTACTCTGTCATTTACATCTTTAATAATTTTATTATTAATTAAATTAAATTAGATCAATTCTTCACCTTTATGTAAGATGATACGGTGTTGcc                                                                                                                                      <  4:26821/139‑1 (MQ=255)

       agtagATATTTGCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTATGCCGGTACGGTGTTGcc                                                                                                                                      >  1:874699/1‑141 (MQ=255)

       agtagATATTTGCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGcc                                                                                                                                      >  4:22442/1‑141 (MQ=255)

       agtagATATTTGCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATGAATTAATTTGTATTCGATCATTTCTTNACCTTTATGTAANATGATACGGTGTTGcc                                                                                                                                      <  1:980586/141‑1 (MQ=255)

         tagATATTTGCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTa                                                                                                                                    >  1:285089/1‑141 (MQ=255)

           gATATTTGCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGc                                                                                                                                       >  2:1040670/1‑136 (MQ=255)

           gATATTTGCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGc                                                                                                                                       <  1:1040670/136‑1 (MQ=255)

            atatTTGCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAAtt                                                                                                                                 <  2:722208/141‑1 (MQ=255)

                 tGCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAat                                                                                                                            >  1:568831/1‑141 (MQ=255)

                  gCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTTATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGCTGTTGCCTAATTTCAata                                                                                                                           >  2:411800/1‑141 (MQ=255)

                  gCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAata                                                                                                                           >  2:46114/1‑141 (MQ=255)

                  gCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAata                                                                                                                           >  4:355098/1‑141 (MQ=255)

                  gCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAata                                                                                                                           >  1:569218/1‑141 (MQ=255)

                  gCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAata                                                                                                                           >  1:242201/1‑141 (MQ=255)

                      tAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGt                                                                                                                       <  4:256917/141‑1 (MQ=255)

                      tAAATGCTATAACATTATTAATAATTGTTATATTTCCTCTGTTATTTACAACTTTACTAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCCTTATGTACAATGACACGGTGTTGCCTAATTTCAATATCGt                                                                                                                       <  4:109166/141‑1 (MQ=255)

                          tGCTGTAACATTATTTATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGGGTTGCCTAATTTCAATATCGTTtct                                                                                                                   >  2:799402/1‑141 (MQ=255)

                           gCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTtgta                                                                                                                  >  2:472866/1‑141 (MQ=255)

                            cTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTtctat                                                                                                                 <  4:193514/141‑1 (MQ=255)

                             tGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATATGGTGTTGCCTAATTTCAATATCGTTtctatc                                                                                                                >  1:304355/1‑141 (MQ=255)

                             tGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTtctatc                                                                                                                >  2:889293/1‑141 (MQ=255)

                              gTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTtctatct                                                                                                               >  1:51825/1‑141 (MQ=255)

                                aaCATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTTATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGCTGTTGCCTAATTTCAATATCGTTTTTATCTAc                                                                                                             <  4:428781/141‑1 (MQ=255)

                                   attattAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAACATGATCCGCTGTTTCCTAATTTCAATATCGTTTCTATCTACTcc                                                                                                          <  2:48606/141‑1 (MQ=255)

                                    ttattaATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCt                                                                                                         <  2:300847/141‑1 (MQ=255)

                                     tattaATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCtt                                                                                                        <  1:458128/141‑1 (MQ=255)

                                      attaATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCttt                                                                                                       >  1:805283/1‑141 (MQ=255)

                                       ttaATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTa                                                                                                      >  2:430529/1‑141 (MQ=255)

                                         aaTCATTGTTATAGTTCCTCGGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTaaa                                                                                                    >  2:495829/1‑141 (MQ=255)

                                          aTCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTaaaa                                                                                                   >  1:930714/1‑141 (MQ=255)

                                          aTCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTaaaa                                                                                                   >  1:242156/1‑141 (MQ=255)

                                            cATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAatc                                                                                                 >  4:126799/1‑141 (MQ=255)

                                            cATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAatc                                                                                                 >  4:237994/1‑141 (MQ=255)

                                                 ttATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTa                                                                                            >  1:295998/1‑141 (MQ=255)

                                                 ttATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTa                                                                                            >  2:868091/1‑141 (MQ=255)

                                                  tataTTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCTTTCATTTCTTCACTTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAATATCATCTAg                                                                                           <  2:241236/141‑1 (MQ=255)

                                                   atatTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGa                                                                                          <  1:31427/141‑1 (MQ=255)

                                                   atatTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGa                                                                                          <  1:185161/141‑1 (MQ=255)

                                                      tttCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAata                                                                                       >  1:969648/1‑141 (MQ=255)

                                                        tCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATAtc                                                                                     <  3:75190/141‑1 (MQ=255)

                                                         ccTCTGTCATTTACAACTTTAATAATTTTATTATTTATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGCTGTTGCCTAATTTCAATATCGTTTTTATCTACACCTTTAAAATCCTCTAGAATAtct                                                                                    >  2:622410/1‑141 (MQ=255)

                                                            ctGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCtt                                                                                 >  2:704106/1‑141 (MQ=255)

                                                             tGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTg                                                                                >  2:675462/1‑141 (MQ=255)

                                                              gTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTga                                                                               <  2:121628/141‑1 (MQ=255)

                                                               tCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTgag                                                                              >  2:624102/1‑141 (MQ=255)

                                                                 aTTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGaa                                                                            >  1:950873/1‑141 (MQ=255)

                                                                   ttACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGaaaa                                                                          <  3:214644/141‑1 (MQ=255)

                                                                      cAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAAtt                                                                       <  3:267154/141‑1 (MQ=255)

                                                                       aaCTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTc                                                                      <  1:253183/141‑1 (MQ=255)

                                                                       aaCTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTc                                                                      >  2:287850/1‑141 (MQ=255)

                                                                        aCTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTTAGAAAAATTCg                                                                     >  2:956996/1‑141 (MQ=255)

                                                                        aCTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCg                                                                     <  1:1024313/141‑1 (MQ=255)

                                                                           ttAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGaaa                                                                  <  2:359146/141‑1 (MQ=255)

                                                                                aatTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTAc                                                             <  3:146579/141‑1 (MQ=255)

                                                                                    ttattattaATTAATTTGTATTCTATCATTTCTTCACCTTTATGTAAAATGATACTGTTTTGCCTATTTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAgtt                                                         <  1:644668/141‑1 (MQ=255)

                                                                                     tattatttATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGCTGTTGCCTAATTTCAATATCGTTTTTATCTACACCTTTAAAATCCTCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAgttg                                                        >  1:34452/1‑141 (MQ=255)

                                                                                     tattattaATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAgttg                                                        >  1:668674/1‑141 (MQ=255)

                                                                                          ttaattaaTTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTaaa                                                   >  2:1061544/1‑141 (MQ=255)

                                                                                             attaattTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAatc                                                                                                 <  1:1027306/92‑1 (MQ=255)

                                                                                             attaattTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATcc                                                >  3:54162/1‑141 (MQ=255)

                                                                                               taattTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCt                                              <  2:684223/141‑1 (MQ=255)

                                                                                                 attTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTAc                                            >  3:91833/1‑141 (MQ=255)

                                                                                                 attTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTAc                                            <  2:969648/141‑1 (MQ=255)

                                                                                                 attTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTAc                                            <  2:1019790/141‑1 (MQ=255)

                                                                                                 attTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGCTGTTGCCTAATTTCAATATCGTTTTTATCTACACCTTTAAAATCCTCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAATTCCCTAc                                            <  3:56757/141‑1 (MQ=255)

                                                                                                   ttGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACtt                                          >  2:360659/1‑141 (MQ=255)

                                                                                                      tATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCa                                       >  2:76855/1‑141 (MQ=255)

                                                                                                      tATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCa                                       >  4:316803/1‑141 (MQ=255)

                                                                                                      tATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCa                                       <  2:657255/141‑1 (MQ=255)

                                                                                                        ttCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCTTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAAt                                     <  2:296076/141‑1 (MQ=255)

                                                                                                          cGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTa                                   <  2:271611/141‑1 (MQ=255)

                                                                                                          cGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGA‑‑ATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGa                                 <  2:242201/141‑1 (MQ=255)

                                                                                                            aTCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGa                                 >  2:656402/1‑141 (MQ=255)

                                                                                                            aTCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGa                                 <  4:411802/141‑1 (MQ=255)

                                                                                                            aTCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGa                                 <  2:670216/141‑1 (MQ=255)

                                                                                                                tttCTTCACCTTTATGTAAAATGATCCGGTGTTGCCTAATTTCAATCTCGTTTATAACTACACATTCAAAATCATCTTGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAg                             <  1:647370/141‑1 (MQ=255)

                                                                                                                tttCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAg                             >  4:36206/1‑141 (MQ=255)

                                                                                                                   cttcACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGttt                          >  1:46382/1‑141 (MQ=255)

                                                                                                                        ccTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTa                     <  4:25011/141‑1 (MQ=255)

                                                                                                                        ccTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCCACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTa                     <  4:18772/141‑1 (MQ=255)

                                                                                                                         cTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTaa                    <  2:799414/141‑1 (MQ=255)

                                                                                                                             aTGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTATAAAAGTTTTACTAAGTCt                >  1:706860/1‑141 (MQ=255)

                                                                                                                               gTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTAATAAGTCTaa              <  1:954598/141‑1 (MQ=255)

                                                                                                                                tAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTAAGTCTAAt             <  1:430529/141‑1 (MQ=255)

                                                                                                                                 aaaaTGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTAAGTCTAAtt            <  2:506208/141‑1 (MQ=255)

                                                                                                                                 aaaaTGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTAAGTCTAAtt            >  4:195363/1‑141 (MQ=255)

                                                                                                                                  aaaTGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTAAGTCTAATTa           >  2:478696/1‑141 (MQ=255)

                                                                                                                                  aaaTGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTAAGTCTAATTa           >  1:644933/1‑141 (MQ=255)

                                                                                                                                  aaaTGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTAAGTCTAATTa           >  3:375649/1‑141 (MQ=255)

                                                                                                                                   aaTGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTAAGTCTAATTaa          >  2:262047/1‑141 (MQ=255)

                                                                                                                                   aaTGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTAAGTCTAATTaa          <  4:417928/141‑1 (MQ=255)

                                                                                                                                        tACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTAAGTCTAATTAAGTCTg     <  3:448574/141‑1 (MQ=255)

                                                                                                                                        tACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTAAGTCTAATTAAGTCTg     <  2:693862/141‑1 (MQ=255)

                                                                                                                                           ggTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTAAGTCTAATTAAGTCTGGaa  >  1:912185/1‑141 (MQ=255)

                                                                                                                                           ggTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTAAGTCTAATTAAGTCTGGaa  >  3:83231/1‑141 (MQ=255)

GATTTATAGTAGATATTTGCAGTAAATGCTGTAACATTATTAATCATTGTTATATTTCCTCTGTCATTTACAACTTTAATAATTTTATTATTAATTAATTTGTATTCGATCATTTCTTCACCTTTATGTAAAATGATACGGTGTTGCCTAATTTCAATATCGTTTCTATCTACTCCTTTAAAATCATCTAGAATATCTCTTGAGAAAAATTCGAAATCTACAGTTGTTAAATCCCTACTTTCAATTAGAAAAGTTTTACTAAGTCTAATTAAGTCTGGAA  >  USA300TCH1516_ALE/1654278‑1654557

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